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INTRODUCTION 

Breast  cancer  remains  one  of  the  most  frequently  diagnosed  cancers  today.  One  in  eight  women  is  expected 
to  present  with  breast  cancer  within  her  lifetime  in  developed  countries.  An  estimated  1,000,000  cases  are 
detected  each  year  worldwide  and  in  Canada  alone,  an  estimated  20,700  women  will  be  diagnosed  in  2002 
with  breast  cancer,  and  5,400  women  will  be  lost  to  this  disease  (1).  For  women  with  recurrent  disease,  the 
median  time  of  survival  is  about  two  years  (2).  Despite  such  striking  statistics,  breast  cancer  related 
mortality  is  slowly  decreasing  as  continuing  research  has  led  to  earlier  detection,  more  treatment  options  for 
breast  cancer  patients  and  an  improved  chance  of  long-term  survival.  Improving  the  diagnosis  and  clinical 
management  of  breast  cancer  requires  access  to  and  characterization  of  biomarkers  that  are  able  to  reflect  the 
molecular  phenotype  of  breast  tissue.  Indeed,  It  is  believed  that  the  knowledge  of  the  genetic  changes  that 
manifest  themselves  as  alteration  of  gene  expression  and  that  occur  during  breast  tumorigenesis  and  breast 
tumor  progression,  will  allow  the  identification  of  targets  for  new  preventive  and  curative  strategies. 

During  the  first  year  of  this  award,  we  focused  our  research  on  two  breast  specific  genes  (mammaglobin  and 
hSBEM)  and  on  genes  involved  in  estrogen  signaling  pathway  {ER-beta  isoforms  and  co-regulators). 

BODY 

Mammaglobin 

Mammaglobin  A  was  first  identified  in  1996,  using  differential  display  analysis,  as  a  breast  specific  member 
of  the  uteroglobin  gene  family  overexpressed  in  some  breast  tumors  (3, 4).  Today,  a  search  for  breast  specific 
ESTs  performed  using  the  Differential  Gene  Expression  Displayer  Tool  at  the  CGAP  website 
(http://cgap.nci.nih.gov/Tissues/GXSl  shows  that  mammaglobin  A  related  ESTs  have  been  found  in  9 
different  breast  cDNA  libraries  but  only  two  non  breast  libraries,  further  confirming  the  relatively  restrained 
breast  specificity  of  mammaglobin  A  expression.  Using  a  subtractive  hybridization  approach,  I  previously 
identified  mammaglobin  A  mRNA  as  overexpressed  in  the  in  situ  compared  to  the  invasive  element  within  an 
individual  breast  tumor  (5, 6).  Further  in  situ  hybridization  analysis  (6),  performed  in  breast  tumors  selected 
to  include  normal,  in  situ  and  invasive  primary  tumor  elements  and  in  most  cases  axillary  lymph  node 
metastases,  revealed  that  mammaglobin  A  expression  can  be  detected  in  all  elements  (normal,  DCIS,  invasive 
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and  nodal  metastasis).  Mammaglobin  A  expression  was  found  to  be  restricted  to  epithelial  cells  and  to  be 
significantly  increased  in  tumor  cells  compared  to  normal  cells  (6).  This  suggests  that  mammaglobin  A 
expression  is  differentially  regulated  in  normal  and  breast  tumor  tissue,  and  could  be  involved  in  the 
mechanisms  underlying  breast  tumorigenesis.  Analysis  of  mammaglobin  A  expression  within  independent 
primary  breast  tumors  and  their  corresponding  axillary  lymph  nodes  further  revealed  that  all  lymph  node 
positive,  but  none  of  the  nodes  negative  for  metastatic  breast  carcinoma  by  histology,  were  mammaglobin  A 
positive  by  reverse  transcription-PCR  (6).  Moreover,  it  has  recently  been  demonstrated  that  circulating 
mammary  carcinoma  cells  can  also  be  detected  in  the  blood  of  breast  cancer  patients  via  PCR  detection  of 
mammaglobin  A  mRNA  (7-10).  Mammaglobin  A  is  now  therefore  considered  as  being  a  specific  marker  of 
axillary  lymph  node  breast  metastases  as  well  as  of  occult  breast  cancer. 

As  shown  in  APPENDIX  1  (Submitted  to  Brit.  J.  Cancer),  we  evaluated  MGBl  expression  at  the  mRNA 
and  at  the  protein  level  by  reverse-transcription  polymerase  chain  reaction  and  immunocytochemistry  in  52 
and  32  breast  tumors,  respectively.  Both  MGBl  mRNA  and  protein  expression  were  significantly  higher  in 
estrogen  receptor  positive  compared  to  estrogen  negative  tumors  (Mann- Whitney  rank  sum  test,  p  =  0.04; 
Chi-square  test,  p  =  0.01;  respectively).  In  contrast,  MGBl  expression  did  not  correlate  with  progesterone 
receptor  levels  or  Notthingham  grade.  As  estrogen  and  antiestrogen  treatment  of  estrogen  positive  breast 
cancer  cell  lines  does  not  modify  MGBl  expression  we  suggest  that  MGBl  could  be  a  new  independent 
breast  cancer  prognostic  marker. 

hSBEM 

There  is  an  urgent  need  to  identify  novel  genes  whose  expression  is  restricted  to  the  mammary  epithelium,  as 
these  genes  have  the  greatest  potential  to  enhance  the  detection  of  micrometastatic  disease  and  the  potential 
to  report  on  proliferative  changes  in  the  breast. 

As  detailed  in  APPENDIX  2  (1 1),  we  have  identified  a  novel  putative  breast-specific  gene  (hSBEM,  Human 
Small  Breast  Epithelial  Mucin),  that  belongs  to  a  recently  regrouped  cluster  (UniGene  identifier  Hs.348419), 
and  which  represents  an  attractive  candidate  for  a  new  breast  tumor  marker  with  obvious  potential  for  cancer 
diagnostics.  We  showed  that  hSBEM  gene  is  expressed  in  most  breast  cancer  cell  lines  and  widely  expressed 
in  epithelial  cells  from  both  normal  and  tumor  breast  tissues.  In  normal  tissues  hSBEM  gene  is  only 
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expressed  in  breast  and  salivary  gland.  Moreover,  hSBEM  mRNA  can  be  detected  in  axillary  lymph  node 
metastases  but  not  in  normal  lymph  nodes,  suggesting  that  hSBEM  mRNA  could  be  used  to  detect  micro- 
metastases. 

Among  the  primary  breast  tumors  examined  in  this  study  (representing  mostly  invasive  ductal  carcinoma), 
SBEM  mRNA  was  observed  by  RT-PCR  analysis  in  the  majority  (>90%)  of  cases.  Despite  a  significant 
overall  correlation  between  the  expression  of  SBEM  and  mammaglobin  A  mRNA,  a  significantly  higher 
SBEM:mammaglobin-l  ratio  was  observed  in  primary  tumors  associated  with  positive  axillary  lymph  nodes 
as  compared  with  node-negative  tumors.  This  was  mostly  attributable  to  a  trend  toward  higher  SBEM 
expression  in  node-positive  tumors.  Although  further  analysis  of  a  larger  number  of  tumors  will  be  required 
to  confirm  these  observations,  this  may  suggest  differences  in  the  biology  of  these  tumors  and  also  a  possible 
role  of  SBEM  and  mammaglobin  A  in  the  mechanisms  involved  in  tumor  metastasis.  Our  findings  indicate, 
however,  that  SBEM  expression  is  a  common  feature  of  breast  cancer  and  can  furthermore  serve  as  a  useful 
marker  for  breast  nodal  metastasis,  both  for  detection  of  micrometastatic  cells  within  lymph  nodes  as  well  as 
in  the  differential  diagnosis  of  the  primary  origin  of  an  unknown  metastasis.  This  potential  is  enhanced  by 
the  conserved  SBEM  expression  in  high  grade  and  ER/PR-negative  tumors  that  are  most  likely  to 
metastasize. 

£R-beta  isoforms 

Estrogens  regulate  the  growth  and  development  of  normal  human  mammary  tissue  and  are  also  involved  in 
breast  tumor  progression  (12).  Indeed,  estrogens  are  thought  to  promote  the  growth  of  breast  tumors  through 
their  mitogenic  effects  on  breast  cancer  cells.  The  ability  of  antiestrogens  such  as  tamoxifen  or  raloxifene  to 
inhibit  estrogenic  action  provides  the  basic  rationale  for  the  use  of  endocrine  therapies.  Estrogen  action  is 
believed  to  be  mediated  mainly  through  two  estrogen  receptors  (ERs):  ER-alpha  (13)  and  ER-61  (14,  15). 
These  two  receptors,  which  are  encoded  by  two  different  mRNAs  containing  eight  exons  each  (16,  17), 
belong  to  the  steroid/thyroid/retinoic  acid  receptor  superfamily  (18).  Several  variant  forms  of  ER-Bl  mRNAs 
have  been  identified  (for  reviews  see  Refs.  19-22).  ER-62  variant,  deleted  of  regions  encoded  by  ER-131  exon 
8  sequences,  has  been  shown  to  heterodimerize  with  both  ER-fil  and  ER-alpha  and  to  inhibit  ER-alpha  DNA 
binding  capability  (23,  24).  The  ability  of  ER  variants  to  potentially  interfere  with  ER-alpha  and  ER-beta 
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signaling  pathways  raised  the  question  of  their  possible  involvement  in  mechanisms  underlying  breast 
tumorigenesis  and  tumor  progression. 

In  collaboration  with  Dr.  Murphy,  we  investigated  the  putative  functional  characteristics  of  human  receptor 
beta  isoforms  (see  APPENDIX  3,  submitted  to  Endocrinology).  We  showed  that  only  ER-  61  was  able  to 
bind  ligand  whereas  all  ER-beta  isoforms  bind  to  DNA  even  though  their  binding  abilities  differ.  ER-beta 
isoforms  inhibition  of  ER-alpha  and  ER-  61  transcriptional  activity  is  promoter  specific.  Overall,  our  data 
suggest  that  ER-beta  isoforms  may  have  a  differentially  modulating  estrogen  action.  We  hypothesize  that  the 
modification  of  the  ratio  between  the  different  ERs  (25, 26),  together  with  these  different  modulating  actions 
play  a  role  in  the  mechanisms  underlying  breast  tumorigenesis  and  tumor  progression. 

Co-regulators 

Recent  observations  using  laboratory  models  (27-30)  have  demonstrated  that  altered  levels  of  ER  isoforms 
and/or  alteration  of  expression  of  co-activators  and  co-repressors  can  result  in  altered  estrogen  and 
antiestrogen  activity  in  target  cells,  suggesting  the  hypothesis  that  altered  levels  of  ER  isoforms  and/or  co¬ 
regulators  in  vivo  could  be  a  mechanism  of  tamoxifen  resistance.  As  mentioned  earlier  in  the  text,  we  have 
previously  demonstrated  that  the  relative  expression  of  ERo/ERp  as  well  as  the  relative  expression  of  some 
ER  coactivators  to  corepressors  is  significantly  altered  during  breast  tumorigenesis  in  vivo  (26,  31). 
Furthermore,  since  these  alterations  parallel  the  marked  changes  in  estrogen  action  that  accompany  breast 
tumorigenesis,  they  may  have  a  role  in  this  process. 

*  In  collaboration  with  Dr.  Murphy  we  explored  the  hypothesis  that  such  changes  could  imderlie  de  novo 
tamoxifen  resistance  in  vivo.  The  expression  of  ER  isoforms,  two  known  co-activators  (steroid  receptor 
RNA  activator,  SRA,  (28)  and  amplified  in  breast  cancer-1,  AIBl  (32))  and  one  co-repressor  (repressor  of 
estrogen  receptor  activity,  REA  (29)  of  ER  activity  have  been  investigated  in  primary  breast  tumors  from 
node  negative  patients  whose  tumors  were  ER  positive  and  that  subsequently  responded  or  had  disease 
progression  while  on  adjuvant  tamoxifen  therapy.  As  shown  in  APPENDIX  4  (In  press,  Brit  J.  Cancer)  we 
found  no  evidence  for  altered  coregulators  expression  in  breast  tumors  that  are  de  novo  resistant.  However, 
our  data  provide  preliminary  evidence  that  the  expression  of  ER-beta  protein  isoforms  may  differ  in  primary 
tumors  of  breast  cancer  patients  who  prove  to  have  differential  sensitivity  to  tamoxifen  therapy. 
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*  When  characterizing  SRA,  which  was  found  to  function  as  an  RNA  transcript  and  to  be  specific  for 
steroid  receptors,  Lanz  et  al.  reported  that  none  of  the  different  SRA  transcript  isoforms  isolated  encoded 
detectable  levels  of  SRA  protein  following  in  vitro  translation  experiments  (28).  We  have  identified  three  new 
SRA  isoforms  encoding  in  vitro  stable  SRA  proteins.  One  of  these  isoforms  contains  a  point  mutation 
followed  by  a  codon  insertion  that  is  found  in  both  normal  and  tumor  breast  tissue  as  well  as  in  some  breast 
tumor  cell  lines.  This  insertion  is  present  at  the  DNA  level  and  is  likely  to  represent  a  genetic 
polymorphism.  We  hypothesize  that  wild-type  and  inserted  SRA  mRNA  and  protein  could  play  a  role  in 
the  mechanisms  underlying  breast  tumorigenesis  and  tumor  progression  (APPENDIX  5,  submitted  BBRC). 

KEY  RESEARCH  ACCOMPLISHMENTfi 

Mammaglobin  expression  correlates  with  ER  levels  but  not  with  progesterone  receptor  levels  or 
Notthingham  grade  in  primaiy  breast  tumors. 

Estrogen/antiestrogen  treatment  does  not  modify  Mammaglobin  expression. 

*  We  identified  a  novel  putative  breast-specific  gene  (hSBEM,  Human  Small  Breast  Epithelial  Mucin). 

We  investigated  the  putative  fimctional  characteristics  of  human  receptor  beta  isoforms. 

We  found  that  expression  of  ER-beta  protein  isoforms  differ  in  primary  tumors  of  breast  cancer  patients 
who  prove  to  have  differential  sensitivity  to  tamoxifen  therapy. 

*  We  have  identified  new  SRA  isoforms  encoding  stable  SRA  proteins. 

REPORTABLE  OUTCOMES 

*  Two  articles  published  or  in  press  (APPENDICES  2  and  4) 

*  Three  articles  submitted  (APPENDICES  1,  3  and  5) 

CONCLUSION 

As  only  two  projects  are  currently  fimded  (SBEM  and  SRA),  the  initial  statement  of  work  will  be  changed 
and  set  up  in  coordination  with  the  Grants  Officer. 
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Abstract 

Mammaglobin  (MGBl)  has  previously  been  found  over-expressed  in  breast  tumors  but 
possible  associations  between  its  expression  and  established  prognostic  tumor 
characteristics  such  as  the  levels  of  estrogen  and  progesterone  receptors  have  not  yet  been 
investigated.  We  evaluated  MGBl  expression  at  the  mRNA  and  at  the  protein  level  by 
reverse-transcription  polymerase  chain  reaction  and  immimocytochemistry  in  52  and  32 
breast  tumors,  respectively.  Both  MGBl  mRNA  and  protein  expression  were  significantly 
higher  in  estrogen  receptor  positive  compared  to  estrogen  negative  tumors  (Mann- Whitney 
rank  sum  test,  p  =  0.04;  Chi-square  test,  p  =  0.01;  respectively).  In  contrast,  MGBl 
expression  did  not  correlate  with  progesterone  receptor  levels  or  Notthingham  grade.  As 
estrogen  and  antiestrogen  treatment  of  estrogen  positive  breast  cancer  cell  lines  does  not 
modify  MGBl  expression  we  suggest  that  MGBl  could  be  a  new  independent  breast 
cancer  prognostic  marker. 

Key  words:  MGBl,  estrogen  receptor,  progesterone  receptor,  grade,  breast  cancer. 
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Introduction 

Mammaglobin  (MGBl)  was  first  identified  in  1996,  using  differential  display  analysis,  as  a 
breast  specific  member  of  the  secretoglobin  (SCGB)  gene  family  overexpressed  in  some 
breast  tumors  (Watson  and  Fleming,  1994;  Watson  and  fleming,  1996).  Today,  a  search  for 
breast  specific  ESTs  performed  using  the  Differential  Gene  Expression  Display er  Tool  at 
the  CGAP  website  (http://cgap.nci.nih.gov/Tissues/GXSl  shows  that  MGBl  related  ESTs 
have  been  found  in  9  different  breast  cDNA  libraries  but  only  two  non  breast  libraries, 
further  confirming  the  relative  breast  specificity  of  MGB 1  expression.  Using  a  subtractive 
hybridization  approach,  we  previously  identified  MGB  1  mRNA  as  overexpressed  in  the  in 
situ  compared  to  the  invasive  element  within  an  individual  breast  tumor  (Leygue  et  al., 
1996a;  Leygue  et  al.,  1999).  Further  in  situ  hybridization  analysis,  performed  in  breast 
tumors  selected  to  include  normal,  in  situ  and  invasive  primaiy  tumor  elements  revealed  that 
MGBl  expression,  restricted  to  epithelial  cells,  could  be  detected  in  all  elements  and  was 
significantly  increased  in  tumor  cells  compared  to  normal  cells  (Leygue  et  al.,  1999).  This 
higher  MGBl  expression  in  malignant  versus  non  malignant  breast  epithelium  has  also  been 
confirmed  at  the  protein  level  by  immunocytochemistry  (Watson  et  al.  1999).  In  this  latter 
study,  Watson  et  al.  concluded  that  MGBl  expression  was  independent  of  tumor  grade  and 
histological  type. 


4 


.  Leygue  et  al.  Appendix  1 

It  has  recently  been  demonstrated  that  circulating  mammary  carcinoma  cells  can  also  be 
detected  in  the  blood  of  breast  cancer  patients  via  PCR  detection  of  MGBl  mRNA  (Zach  et 
aJ.,  1999;  Zach  et  al.  2001;  Gal  et  al.,  2001;  Silva  et  al.,  2001).  Even  though  its  biological 
function  remains  unknown,  MGBl  is  now  considered  as  a  relatively  specific  marker  of 
axillary  lymph  node  breast  metastases  as  well  as  of  occult  breast  cancer  (Kataoka  et  al., 
2000;  Ooka  et  al.,  2000;  Silva  et  al.,  2002;  Marchetti  et  al.,  2001).  Interestingly,  Zach  et  al. 
detected  MGBl  mRNA  expression  by  nested  reverse-transcription  PCR  (RT-PCR)  more 
frequently  in  the  blood  of  patients  with  estrogen  receptor  positive  (ER+)  breast  tumor  than 
in  the  blood  of  estrogen  receptor  negative  (ER-)  breast  cancer  patients,  suggesting  a  possible 
relationship  between  MGBl  and  ER  levels  in  primary  breast  tumors  (Zach  et  al.,  1999).  In 
order  to  investigate  further  possible  associations  between  MGBl  expression  and  estrogen 
and  progesterone  receptors  in  primary  breast  tumors,  we  assessed  MGB 1  expression  at  the 
mRNA  and  at  the  protein  level  in  a  cohort  of  breast  tumors. 

Materials  and  methods 
Human  breast  tissues  and  cell  lines 

All  breast  tumor  cases  used  for  this  study  were  selected  from  the  NCIC-Manitoba  Breast 
Tumor  Bank  (Winnipeg,  Manitoba,  Canada).  As  it  has  been  previously  described  (Hiller  et 
al.,  1996),  tissues  are  accrued  to  the  Bank  from  cases  at  multiple  centers  within  Manitoba, 
rapidly  collected  and  processed  to  create  matched  formalin-fixed-embedded  and  frozen 
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tissue  blocks  for  each  case  with  the  mirror  image  surfaces  oriented  by  colored  inks.  The 
histology  of  every  sample  in  the  Bank  is  uniformly  interpreted  by  a  pathologist  in 
Hematoxylin/Eosin  (H&E)  stained  sections  from  the  face  of  the  paraffin  tissue  block.  This 
information  is  available  in  a  computerized  database  along  with  relevant  pathological  and 
clinical  information  and  was  used  as  a  guide  for  selection  of  specific  paraffin  and  frozen 
blocks.  Fifty  two  tumors  were  selected,  spanning  a  wide  range  of  estrogen  and  progesterone 
receptor  levels,  as  determined  by  ligand  binding  assay.  Within  these  tumors,  9  were  ER- 
/PR-  (ER  <  3  finol/mg  total  protein;  PR  <  10  finol/mg),  10  were  ER-  /PR+  (ER  <  3 
finol/mg;  PR  >  1 0  finol/mg),  1 0  were  ER+  /  PR-  (ER  >  3  finol/mg;  PR  <  1 0  finol/mg)  and  23 
were  ER+/  PR+  (ER  >  3  finol/mg,  PR  >  10  finol/mg).  These  tumors  also  spanned  a  wide 
range  of  Nottingham  grade  for  ER-  (n  =  19,  grade  ranging  from  5  to  9,  median  8)  and  ER+  (n 
=  33,  grade  ranging  from  5  to  9,  median  6)  tumors.  MGB 1  mRNA  expression  was  assessed 
by  RT-PCR  on  total  RNA  extracted  from  frozen  tissue  sections.  ParaflSn  blocks 
corresponding  to  32  out  of  these  52  tumors  were  sectioned  and  processed  for 
immunohistochemical  analysis  of  MGB  1  expression. 

ZR-75  cells,  ER  positive  breast  cancer  cells  known  to  express  MGBl,  were  grown,  and 
treated  with  Estradiol- 17p  lO'^M  in  charcoal-stripped  medium  or  with  the  antiestrogen  ICI 
182,780  (10  ^M)  in  regular  medium  for  6h,  24  h  or  48  h,  as  previously  described  (Courts  et 
al.,  1999).  Total  RNA  was  extracted  from  frozen  tissue  sections  or  cell  lines  using  Tri¬ 
reagent  (MRCI,  Cincinnati,  OH). 
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RT-PCR  analysis 

One  ng  of  total  RNA  was  reverse  transcribed  in  a  final  volume  of  20  pi  and  one  pi  of  the 
reaction  mixture  subsequently  amplified  by  PCR  as  previously  described  (Leygue  et  al., 
1996b;  Leygue  et  al.,  1996c).  Primers  used  corresponded  to:  MGBl  (sense  5'- 
CCGACAGCAGCAGCCTCAC-3',  located  in  MGBl  sequence  between  bases  41-59,  and 
antisense  5'-TCCGTAGTTGGTTTCTCAC-3',  located  between  bases  401-383  (Watson  et 
al.,  1996);  glyceraldehyde-3 -phosphate  dehydrogenase  (GAPDH)  gene  (sense  5’- 
ACCCACTCCTCCACCTTTG-3’  and  antisense  5'-CTCTTGTGCTCTTGCTGGG-3'); 
and  to  psoriasin  gene  (sense  5'-AAGAAAGATGAGCAACAC-3'  and  antisense  5'- 
CCAGCAAGGACAGAAACT-3')  To  amplify  cDNA  corresponding  to  MGBl,  GAPDH 
and  psoriasin,  30  cycles  (30  seconds  at  94°C,  30  seconds  at  55°C  and  30  seconds  at  72°C) 
of  PCR  were  used.  Ten  microliters  of  PCR  products  were  loaded  on  prestained  (15  pg/ml 
ethidium  bromide)  2%  agarose  gels.  Identity  of  fragments  corresponding  to  MGBl, 
GAPDH  and  psoriasin  had  previously  been  confirmed  by  sequencing. 

Three  independent  PCRs  were  performed  using  MGB 1 ,  psoriasin  and  GAPDH  primers  and 
signals,  visualized  with  U.V.  irradiation  on  a  GelDoc2000/ChemiDoc  System  (Biorad),  were 
quantified  by  densitometry  using  the  Quantity  One  software  (Version  4.2,  Biorad).  MGB  1 
and  psoriasin  expression  was  expressed  relative  to  GAPDH  expression  as  previously 
described  (Leygue  et  al,  1998).  Briefly,  three  independent  PCRs  were  performed  using  each 
set  of  primers.  In  order  to  control  for  variations  between  experiments,  a  value  of  1  was 
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arbitrarily  assigned  to  the  signal  of  one  particular  tumor  measured  in  each  set  of  PCR 
experiments  (always  the  same  tumor)  and  all  signals  were  expressed  relative  to  this  signal. 
Levels  of  MGBl  were  then  expressed  relative  to  the  GAPDH  signal  corresponding  to  each 
individual  tumor  sample.  Correlation  between  normalized  MGBl  expression  and  tumor 
characteristics  was  tested  by  calculation  of  the  Spearman  coefficient,  r.  Comparison 
between  tumor  subgroups  was  performed  using  the  Maim- Whitney  rank  sum  test,  two 
sided. 

Inununohistochenucal  analysis  of  MGBl  expression. 

Detection  of  MGBl  protein  was  performed  using  an  antibody  previously  characterized  and 
kindly  provided  by  Dr.  Timothy  Fleming  (see  Watson  et  al.,  1999).  Paraffin-embedded 
breast  tissue  sections  were  processed  using  the  automated  Discovery  Staining  Module, 
Ventana  System  (Tucson,  Arizona)  and  the  Research  IHC  DAB  paraffin  protocol  according 
to  the  manufacturer's  instructions.  All  steps  were  performed  automatically:  briefly, 
following  deparaffination  of  tissue  sections,  slides  were  incubated  60  minutes  at  42°C  in  the 
presence  of  rabbit  anti-MGBl  antibody  (1/1000  final  concentration),  washed,  incubated 
with  Biotinylated  secondary  anti-rabbit  antibody  (14  minutes  42"C),  washed,  incubated  8 
minutes  with  Avidin-HRPO  complex  subsequently  detected  with  DAB-H2O2  solution. 
Counter-staining  was  also  performed  automatically  by  the  Ventana  apparatus 
(hematoxylin/bluing  reagent). 

Levels  of  mammaglobin  expression  were  assessed  by  brightfield  microscopic  examination  at 
low  power  magnification  and  using  a  previously  described  semi-quantitative  approach 
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(Leygue  et  al.,  1998).  Scores  were  obtained  by  estimating  average  signal  intensity  (on  a  scale 
of  0  to  3)  and  the  proportion  of  epithelial  cells  showing  a  positive  signal  (0,  none;  0.1,  less 
than  one  tenth,  0.5,  less  than  one  half;  1.0  greater  than  one  half).  The  intensity  and 
proportion  scores  were  then  multiplied  to  give  an  overall  score.  Cases  with  a  score  lower 
than  or  equal  to  1  were  considered  negative  or  weakly  positive,  whereas  tumors  with  scores 
higher  than  1.0  were  classified  as  positive  for  MGBl  expression.  Statistical  comparisons 
between  tumor  subgroups  have  been  performed  using  the  Chi-square  test. 

Results 

Assessment  of  MGBl  mKNA  expression  in  a  cohort  of  52  human  breast  tumor  samples. 

To  establish  whether  MGBl  mRNA  expression  paralleled  established  known  prognostic 
parameters  such  as  ER  and  PR  levels,  a  cohort  of  52  cases  was  selected  from  the  NCIC- 
Manitoba  Breast  Tumor  Bank.  For  each  case,  clinical  characteristics  of  the  tumor  (i.e  ER 
and  PR  levels,  Nottingham  grade)  were  known  (see  Materials  and  methods  section  for  a 
summary  of  tumor  subgroup  characteristics).  Total  RNA  was  extracted  from  frozen 
primary  tumor  sections,  reverse-transcribed  and  analyzed  by  RT-PCR  using  primers 
recogmzing  specifically  MGBl  cDNA  and  chosen  to  span  intronic  regions.  As  shown 
Figure  1,  MGBl  corresponding  signal  can  be  detected  in  the  majority  of  cases,  even  though 
levels  of  expression  varied  from  one  sample  to  another.  Amplification  of  the  ubiquitously 
expressed  GAPDH  cDNA  in  the  same  cDNA  samples  was  performed  in  parallel  and  for 
each  case,  a  normalized  MGB 1  mRNA  expression  value  was  calculated  (see  Materials  and 
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Methods).  MGBl  expression  was  found  to  strongly  correlate  with  ER  levels  (n  =  52, 
Spearmann  coefficient  r  =  0.282,  p  =  0.042)  but  not  with  PR  levels  or  Grade  (data  not 
shown).  Similarly  (Figure  2),  using  the  established  clinical  cut-off  of  ER  positivity  (ER 
positive  tumors  have  a  binding  higher  than  3  finol/mg  of  total  protein),  MGBl  mRNA 
expression  was  significantly  (Mann- Whitney  rank  sum  test,  two  sided,  p  =  0.040)  higher  in 
ER+  (n  =  33,  median  value  MGBl  =  0.62)  than  in  ER  -  (n  =  19,  median  MGBl  value  = 
0.33).  Interestingly,  even  though  not  statistically  significant,  a  trend  towards  a  higher 
MGBl  expression  in  PR  +  (median  value  MGBl  =  0.56,  n  =33)  than  in  PR-  (median  value 
MGBl  =  0.36,  n  =  19)  was  observed.  However,  as  this  trend  disappears  when  comparing 
PR+  and  PR-  subgroups  within  ER-  and  ER+  cohorts,  we  concluded  that  this  apparent 
trend  resulted  from  sampling  biases  (data  not  shown). 

Similarly,  the  apparent  trend  of  a  higher  MGBl  expression  in  Low  grade  tumors  (median 
values  MGBl :  Low  grade  =  0.56,  n  =29;  High  grade  =  0.36,  n  =  23)  likely  results  from  the 
strong  inverse  correlation  between  grade  and  ER  levels  observed  within  the  cohort. 

Assessment  of  MGBl  protein  expression  in  a  cohort  of  32  human  breast  tumor  san^les. 

In  order  to  determine  whether  MGBl  protein  expression  correlated  with  MGBl  mRNA 
expression  and  whether  a  similar  association  between  ER  status  and  MGBl  expression 
could  be  observed  at  the  protein  level,  paraffin  blocks  corresponding  to  32  out  of  these  52 
tumors  were  sectioned  and  processed  for  immunohistochemical  analysis  of  MGBl 
expression  (see  Materials  and  Methods).  Slides  were  scored  blindly  for  MGBl  protein 
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expression  by  a  pathologist  as  described  in  the  Materials  and  Methods  section.  Some 
sections  showed  no  (Figure  3  A,  MGBl  score  =  0)  or  low  (Figure  3B,  MGBl  score  =  1) 
MGBl  expression  whereas  others  presented  strong  MGBl  protein  signal  (Figure  3C, 
MGBl  score  =  3,  Figure  3-D,  MGBl  score  =  2).  Comparison  of  MGBl  protein  scores  and 
previously  obtained  normalized  MGBl  mRNA  levels  revealed  a  strong  correlation  (n  =  32, 
Spearman  r  coefficient  r  =  0.575,  p  =  0.0006)  between  protein  and  mRNA  levels.  Tumors 
were  classified  as  low  (scores  between  0  and  1)  and  high  (1.5  to  3)  MGBl  protein 
expressors,  and  differences  between  tumor  subgroups  (ER+/ER-,  PR+/PR-,  Low  grade/high 
grade)  were  assessed  using  Chi-square  test.  As  observed  for  MGBl  mRNA,  MGBl  protein 
positivity  was  associated  (Chi-square  test,  p  =  0.017)  with  ER  status  but  not  with  PR 
status  or  grade  (Figure  4). 

Absence  of  estrogen  regulation  of  MGBl  expression. 

These  data  suggested  that  estrogen  might  regulate  MGBl  expression.  In  order  to  address  the 
question  of  a  possible  regulation  of  MGBl  expression  in  breast  cancer  cells,  ZR-75  cells, 
known  to  express  MGBl  (Min  et  al.,  1998)  were  treated  by  estradiol-17P  10*M  or  the 
antiestrogen  ICI- 182,780  lO  '^M  for  6,  24  and  48  hours  as  described  in  the  Materials  and 
Methods  section.  Total  RNA  was  extracted  and  analyzed  by  RT-PCR  using  primers 
recognizing  either  GAPDH,  MGBl  or  psoriasin  cDNAs.  Psoriasin  was  chosen  as  its 
expression  has  previously  been  shown  to  be  regulated  by  estrogen  treatment  (Moog-Lutz  et 
al.,  1995).  MGBl  mRNA  expression  was  not  changed  under  any  treatment  condition  (data 
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not  shown),  whereas,  as  expected,  psoriasin  signal  was  found  to  be  increased  by  estradiol 
and  decreased  by  antiestrogen  treatment  as  soon  as  6  hours  of  treatment,  with  a  maximum 
effect  after  24  and  48  hours  of  treatment  (estradiol  treatment:  1.5,  2.8  and  4.5  fold  control 
and  antiestrogen  treatement  0.90, 0.80,  and  0.70  fold  control,  respectively). 

Discussion 

Assessment  of  MGBl  expression  at  the  mRNA  and  at  the  protein  levels  in  a  cohort  of 
breast  tissue  samples  showed  a  statistically  significant  relationship  between  MGB 1  levels 
and  ER  status.  However,  within  the  same  cohort,  no  association  was  found  between 
MGBl  expression  and  other  known  prognostic  marker  such  as  PR  levels  or  Notthingham 
grade. 

To  the  clinician,  a  factor  is  considered  a  prognostic  factor  when  it  is  associated  wdth  the 
outcome  of  the  disease,  i.e  predicts  how  the  disease  would  evolve  if  not  treated,  whereas  a 
predictive  factor  is  associated  to  the  degree  of  response  to  therapy,  i.e  predicts  the 
likelihood  of  response  to  a  particular  treatment.  A  high  level  of  ER  in  tumor  tissue  has  a 
good  prognostic  value  and  also  predicts  a  good  likelihood  of  responding  to  hormonal 
adjuvant  therapy  such  as  tamoxifen  (Aapro,  2001;  Bundred,  2001).  As  PR  expression  is 
positively  regulated  by  estrogens,  higher  PR  levels  in  ER+  tumors  support  the  hypothesis 
of  an  operational  ER  signaling  pathway  and  is  therefore  also  considered  as  a  good 
prognostic  and  predictive  parameter.  Whereas  the  parallel  between  MGBl  and  ER 
expression  suggested  that  MGB  1  could  be  a  new  ER  target  gene,  the  lack  of  association 
with  a  known  regulated  gene  such  as  PR  suggested  that  MGBl  expression  was 
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independent  of  ER  signaling  pathway.  This  later  hypothesis  was  further  supported  by  the 
absence  of  estrogen  and  antiestrogen  regulation  of  MGBl  expression  in  ZR-75  cells,  even 
though  ER  signaling  pathway  appears  functional,  as  shown  by  the  induction  of  a  known 
ER  regulated  gene,  psoriasin.  It  should  be  noted  that  a  similar  absence  of  regulation  was 
also  observed  in  another  ER  positive  breast  cancer  cell  line  MCF-7  cells  (our  unpublished 
results  and  Watson  et  al.,  1998).  However,  even  though  MGBl  gene  was  not  grossly 
rearranged  in  MCF-7  cells  (Watson  et  al.,  1998),  these  cells  do  not  express  endogenous 
MGBl  (our  unpublished  results  and  Watson  et  al.,  1996).  It  might  therefore  be 
hypothesized  that  MGBl  expression  in  MCF-7  cells  is  negatively  regulated  by  other 
factors,  resulting  in  an  absence  of  estrogen  regulation  in  these  cells. 

Interestingly,  the  general  expression  of  MGBl  as  well  its  association  with  ER  levels 
observed  in  vivo  in  breast  tissue  contrasts  with  in  vitro  observations  made  on  mammary 
epithelial  cancer  cell  lines.  Indeed,  looking  at  a  panel  of  different  breast  cancer  cell  lines, 
Watson  et  al.  reported  the  detection  of  MGBl  transcripts  only  in  few  cell  lines  (MB361, 
MB415,  MB468,  BT474,  MB  175)  with  no  expression  in  MCF7,  MB  134,  MB231  or 
MCFIOA  cells  (Watson  et  al.,  1996).  Similarly,  we  did  not  detect  MGBl  expression  in 
breast  cell  lines  such  as  BT20,  T47D  or  MCFIOATI  even  though  a  strong  signal  was  seen 
in  ZR-75  (our  unpublished  observation).  As  cells  such  as  MB468  and  MB361  are  ER- 
and  cells  such  as  ZR-75  or  BT474  are  ER+,  MGBl  expression  does  not  appear  related  to 
ER  status  in  cells  grown  in  vitro.  Overall,  this  suggests  that  most  of  cell  lines,  through 
selection,  medium  conditions  and/or  dedififerentiation  lost  their  ability  to  express  MGBl 
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Presently,  no  data  are  available  regarding  the  possible  biological  function  of  MGBl.  It  has 
however  recently  been  reported  that  MGBl  existed  in  a  tetrameric  complex  with  BUI 01  ( 
lipophilin  B),  another  member  of  the  secretoglobin  family  expression  of  which  correlated 
with  MGBl  expression  in  breast  tissue  (Colpitts  et  al.,  2001).  The  role  of  this  complex  as 
well  as  the  possible  regulation  of  its  components  remains  to  be  determined. 

In  conclusion,  we  found  that  MGBl  expression  correlated  with  ER  levels  in  breast  tumor 
tissue.  As  ER  is  considered  as  a  good  prognostic  factor  and  as  MGBl  does  not  appear  to 
be  directly  regulated  by  the  ER  signaling  pathway,  we  hypothesize  that  MGBl 
expression  could  be  a  new  independent  prognostic  marker  in  breast  cancer. 
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Figure  legend 
Figure  1 

RT-PCR  analysis  of  MGBl  and  GAPDH  mRNA  expression  in  primary  breast  tumors. 
Total  RNA  was  extracted  from  frozen  tissue  sections  corresponding  to  ER  positive  (ER+) 
and  ER  negative  (ER-)  cases,  reverse-transcribed  and  PCR  amplified  as  described  in  the 
Materials  and  Methods  section  using  MGBl  or  GAPDH  specific  primers.  PCR  products 
were  then  separated  on  2%  agarose  gels  prestained  with  ethidium  bromide.  .  Black  arrow: 
product  corresponding  to  MGBl,  grey  arrow:  product  corresponding  to  GAPDH.  M: 
Molecular  weight  marker  (0x174  RF  DNA/Hae  III  fragments,  Gibco  BRL,  Grand  Island, 
NY). 

Figure  2 

Quantification  of  MGBl  mRNA  expression  in  different  breast  tumor  subgroups. 

Total  RNA  was  extracted  from  frozen  tissue  sections  corresponding  to  52  cases  and 
analyzed  as  described  in  Figure  1.  MGBl  mRNA  expression  was  quantified  relative  to 
GAPDH  mRNA  as  described  in  the  Materials  and  Methods  section.  Tumors  were  grouped 
according  to  their  ER  status  (ER+,  ER-),  their  PR  status  (PR+,  PR-)  or  their  Grade  (Low 
Gr:  Notthingham  scores  between  5  and  7,  High  Gr:  Notthingham  scores  between  8  and  9). 
Difference  between  subgroups  were  tested  using  the  Mann- Whitney  rank  sum  test,  two 


sided. 
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Figure  3 

Detection  of  MGBl  protein  in  breast  tumors  by  immunohistochemistry. 

MGBl  protein  was  detected  on  paraffin-embedded  breast  tumor  tissue  sections  using  a 
rabbit  polyclonal  primary  anti-MGBl  antibody  (Kindly  provided  by  Dr.  T  Fleming)  and 
the  Ventana-Discovery  system  as  described  in  the  Materials  and  Methods  section.  Panel  A 
and  B  :  Two  independent  ER-  cases  showing  no  (A)  or  low  (B)  MGBl  protein  expression 
in  tumor  epithelial  cells.  Panel  C  and  D:  two  independent  ER  +  cases  presenting  a  strong 
MGBl  signal.  Blue  bar:  20  pm. 

Figure  4 

Quantification  of  MGBl  protein  in  breast  tumor  subgroups. 

Paraffin-embedded  tissue  section  corresponding  to  32  cases  were  processed  as  shown 
Figure  3.  Slides  were  independently  reviewed  and  scored  as  described  in  the  Materials  and 
Methods  section.  For  each  tumor  subgroup  (ER-,  ER+,  PR+,  PR-,  Low  Grade  and  High 
Grade)  the  number  of  cases  negative  (White  coliunns)  or  positive  (black  columns)  is  shown. 
Differences  between  subgroups  were  tested  using  the  Chi-square  test. 
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Abstract 

Expression  profiling  using  the  public  expressed  sequence  tag  (EST)  and 
serial  analysis  of  gene  expression  (SAGE)  databases  resulted  in  the  iden¬ 
tification  of  a  putative  breast-specific  mRNA  that  we  have  termed  small 
breast  epithelial  mucin  (SBEM).  Hybridization  analysis  performed  on  43 
normal  human  tissues  revealed  that  the  gene  was  only  expressed  in 

mammary  and  salivary  glands.  Further  reverse-transcription  PCR  anal¬ 
yses  confirmed  SBEM  expression  in  most  of  established  human  breast 
epithelial  cell  lines  analyzed  (7  of  8)  but  not  in  cell  lines  of  non-breast 
origin  (0  of  5!B£d/mRNA  expression  was  detected  in  >90%  of  invasive 
ductal  carcinomas  and  correlated  with  the  expression  of  a  previously 
characterized  breast-specific  gene,  mammaglobin-1  (n  ==  54;  Spearman 
r  =  0J4,  P  =  0.011).  Interestingly,  a  higher  SBEM:mammaglobm-l  ratio 
was  observed  in  primary  tumors  with  axillary  lymph  node  metastasis  fiian 
in  node-negative  tumors  (n  =  46;  Mann-Whitney,  P  -  0.04).  In  a  subset 
of  20  primaty  breast  tumors  and  their  matched  axillary  lymph  nodes,  a 
high  concordance  (Fisher’s  exact  test,  P  <  0,001)  was  seen  between  PCR 
detection  of  SBEM  mRNA  in  lymph  node  tissue  and  their  histopatholog- 
kal  status,  indicating  that  SBEM  mRNA  expression  is  conserved  in  nodal 
metastasis.  The  SBEM  gene  is  predicted  to  code  for  a  putative  low 
molecular  weight,  secreted  sialoglycoprotein,  potentially  useful  for  the 
diagnosis  of  metastatic  breast  cancer. 

Introduction 

Early  detection  remains  a  central  goal  in  breast  cancer  treatment  to 
enable  intervention  at  a  localized  and  potentially  curable  stage  and  to 
maximize  the  opportunity  for  breast  conservation.  The  5-year  survival 
rate  for  women  with  breast  cancer  increases  dramatically  when  it  can 
be  diagnosed  at  an  early  stage,  from  >95%  in  patients  with  a  localized 
tumor  to  ^-75%  with  regional  disease  and  <25%  in  women  with 
disseminated  cancer  (1).  Nevertheless,  only  60%  of  all  breast  cancers 
are  diagnosed  at  a  local  stage,  and  any  inprovement  in  early  detection 
would  have  a  significant  inq)act  on  reducing  overall  breast  cancer 
mortality. 

Improving  the  diagnosis  and  clinical  management  of  breast  cancer 
requires  access  to  a  wider  range  of  biomarkers  able  to  reflect  the 
molecular  phenotype  of  breast  tissue.  A  special  need  exists  to  identify 
novel  genes  whose  expression  is  restricted  to  the  mammaiy  epithe¬ 
lium,  because  these  genes  have  the  greatest  potential  to  enhance 
detection  of  micrometastatic  disease  and  the  potential  to  report  on 
proliferative  changes  in  the  breast,  analogous  to  the  ability  of  elevated 
serum  prostate-specific  antigen  levels  to  indicate  the  presence  of 
hyperplasia  or  cancer  of  the  prostate  gland  (2). 
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The  identification  of  new  tissue-specific  markers  has  benefited 
especially  from  expansion  of  public  and  private  databases  for  ESTs"^ 
(3, 4)  and  by  large-scale  efforts  to  profile  patterns  of  gene  expression 
using  techniques  such  as  serial  analysis  of  gene  expression  (5).  Using 
sequence  analysis  software  and  web-based  tools  developed  for  mo¬ 
lecular  profiling,  we  have  identified  a  novel  putative  breast-specific 
gene,  belonging  to  a  recently  regrouped  cluster  (UniGene  identifier 
Hs.348419),^  which  represents  an  attractive  candidate  for  a  breast 
tumor  marker  with  obvious  potential  for  cancer  diagnostics. 

Materials  and  Methods 

Database  and  Sequence  Analysis.  The  cDNA  xProfilor  tool®  was  used  to 
search  for  novel  breast-specific  ESTs.  Protein  sequence  analysis  used  the 
Signal?  algorithm^  to  search  for  the  presence  of  a  signal  sequence  (6)  and  the 
NetOGlyc  algorithm®  to  predict  sites  of  potential  glycosylation  (7). 

RNA  Hybridization  Analysis  for  Tissue  Specificity.  A  ^^P-labeled 
SBEM  probe,  generated  using  the  cloned  SBEM  PCR  product  (396  bp)  and  tiie 
RadPrime  DNA  labeling  system  (Life  Technologies,  Inc.,  Burlington,  Ontario, 
Canada),  was  hybridized  to  a  commercially  available  RNA  Master  Blot  (Clon- 
tech,  Palo  Alto,  CA),  containing  poly(A)4-  RNA  (100-500  ng)  isolated  from 
a  variety  of  adult  and  fetal  human  tissues,  according  to  die  manufecturer’s 
instructions. 

Cell  Culture  and  RNA  Preparation.  Cell  lines  were  obtained  from  the 
American  Type  Culture  Collection  or  other  sources  and  were  cultured  as 
follows:  DMEM  with  10%  fetal  bovine  serum  (MCF7,  MCFIOATI, 
MCFlOATSc,  SK-UT-IB,  and  HepG2);  DMEM  with  10%  calf  serum  (MDA 
MB-231,  Hec  lA,  and  HeLa);  DMEM:Ham’s  F12  (1:1)  with  10%  fetal  bovine 
serum  (ZR-75-1  and  RL95-2);  RPMI  1640  with  10%  fetal  bovine  serum 
(T-47D  and  LNCaP);  or  MSU-1  medium  (8)  with  5%  fetal  bovine  serum 
(M13SV-1).  All  media  were  supplemented  with  penicillin  (50  units/ml),  strep¬ 
tomycin  (50  ixg/ml),  HEPES  (pH  7.4;  5  mM),  and  glutamine  (2  mM).  MCF7, 
T-47D,  and  ZR-75-1  cells  also  received  bovine  insulin  (10  ^g/ml).  Media  and 
sera  were  obtained  from  Life  Technologies,  Inc.  (Life  Technologies,  Inc., 
Grand  Island,  NY).  RNA  was  extracted  from  cultured  cells  using  guanidinium 
isothiocyanate,  followed  by  centrifrxgation  through  a  5.7  M  cesium  chloride 
cushion  as  described  (9).  RNA  from  cultured  primary  HMECs  obtained  by 
reduction  mammoplasty  was  a  kind  gift  from  P.  Ervin  (Biotherapies,  Inc.,  Ann 
Arbor,  MI). 

Breast  Tumors  and  Axillary  Lymph  Nodes.  Fifty- four  invasive  ductal 
carcinomas  were  selected  from  the  Manitoba  Breast  Tumor  Bank  (Winnipeg, 
Manitoba,  Canada).  Cases  spanned  many  ER  (0-298  ftno!/mg  protein)  and  PR 
(0-1199  frnol/mg  protein)  levels,  as  determined  by  ligand  binding  assay. 
Tumors  also  spaimed  many  grades  (Nottingham  grade  scores  from  5  to  9).  For 


^  The  abbreviatioDs  used  are:  EST,  expressed  sequence  tag;  SBEM,  small  breast 
epithelial  mucin;  HMEC,  human  mammary  epithelial  cell;  ER,  estrogen  receptor-o;  PR, 
progesterone  receptor  RT-PCR,  reverse  tianscription-PCR;  GAPDH,  glyceraldehyde 
phosphate  dehydrogenase;  MUCl,  mucin  1. 

*  Internet  address:  litlp://www.ncbi.olm.ttih.gov/UalGeoe/.  UniGene  is  a  system  for 
automatically  partitioning  GenBank  sequences,  including  ESTs,  into  a  nonredundaot  set 
of  gene-oriented  clusters. 

®  Internet  address:  http://cgap.nci.mh.gov/CGAP/Tissues/xProfiler. 

’  Internet  address:  http://www.cbs.dtu.dk/services/SignalP/. 

*  Internet  address:  ht^://www.cbs.dtu.dk/scrvices/NetOGIyc/. 
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Fig.  1.  Tissue  expression  of  SBEM  mRNA.  An 
RNA  Master  Blot  (Clontech  Laboratories,  Palo 
Alto,  G\)  containing  po]y(A)+  RNAs  from  differ¬ 
ent  human  tissues  (A)  was  screened  with  ^^P-labeled 
SBEM  probe  as  described  in  **Materials  and  Meth¬ 
ods.’*  SBEM  transcript  (B)  was  expressed  in  the 
human  salivary  (DT)  and  mammary  gland  (DS). 


46  tumors,  pathological  axillary  lymph  node  status  (presence  or  absence  of 
metastasis)  was  known. 

In  a  subset  of  cases  (n  =  20),  frozen  primary  human  breast  tumor  samples 
and  their  matched  frozen  lymph  nodes  containing  (n  =  14)  or  not  (n  =  6) 
histologically  detectable  metastatic  cancer  cells  were  available  from  the  Mani¬ 
toba  Breast  Tumor  Bank.  For  the  primary  tumor  samples,  the  ER  levels, 
determined  by  ligand  binding  assays,  ranged  from  2.3  finol/mg  protein  to  298 
finol/mg  protein,  whereas  PR  levels  ranged  from  10.1  finol/mg  protein  to  112 
finol/mg  protein. 

RNA  Analysis  by  RT-PCR.  Total  RNA  was  extracted  from  20-/im  frozen 
tissue  sections  (five  sections/tumor)  and  reverse  transcribed  as  described 
previously  (10).  The  primers  used  for  SBEM  amplification  consisted  of 
SBEM-U  (5'-CTTTGAAGCATTTTTGTCTGTG-3';  sense)  and  SBEM-L  (5'- 
AAGGTAAGTAGTTGGATGAAAT-3';  antisense).  PCR  amplifications  were 
performed,  and  PCR  products  were  analyzed  as  described  previously  (1 1),  with 
minor  modifications.  Briefiy,  aliquots  of  each  reverse  transcription  mixture  (2 
|il  for  Fig.  2  and  0.8  fil  for  Fig.  3)  were  amplified  in  a  final  volume  of  20  /il, 
in  the  presence  of  20  mM  Tris-HCl  (pH  8.4),  50  mM  KCl,  1 .5  mM  Mga2, 200 
jjM  of  each  deoxynuclcotide  tripho^hate,  5  ng/jd  of  each  SBEM  primer,  and 
0.5  unit  of  Taq  DNA  polymerase,  ^ch  PCR  consisted  of  35  cycles  (15  s  at 
94°C,  30  s  at  58°C,  and  60  s  at  74®C). 

Primers  used  for  mammaglobm-1  were:  Mam-1  (5'-CCGACAGCAG- 
CAGCCTCAC-3',  sense  strand)  and  Mam-2  (5'-TCCGTAGTTGGTTTCT- 
CAC-3\  antiseose  strand).  Primers  for  the  ubiquitously  expressed  GAPDH 
gene  were  GAP- 1  (5'-ACCCACTCCTCCACCTTTG-3\  sense  strand)  and 
GAP-2  (5^-CTCTTGTGCTCTTGCTGGG-3',  antis^e  strand).  To  amplify 
cDNA  corresponding  to  mammaglobin-1  and  GAPDH,  30  cycles  of  PCR  were 
used  (30  s  at  94*C,  30  s  at  55°C,  and  30  s  at  72°C).  All  buffers  were  the  same 
as  for  SBEM  PCR,  except  that  2  mM  MgCl2  was  used  when  amplifying 
mammaglobin-1  cDNA.  PCR  products  were  fiien  separated  on  a  1.5%  agarose 
gel.  After  electrophoresis,  the  gels  were  stained  with  ethidium  bromide  (0.5 
ftg/ml). 

Quantification  and  Statistical  Analysis.  Three  independent  PCRs  were 
performed  for  tumor  specimens  using  SBEM,  mammaglobin-1,  and  GAPDH 
primers  and  signals,  visualized  with  UV  irradiation  on  a  GeIDoc2000/Chemi- 
Doc  System  (Bio-Rad),  were  quantified  by  densitometry  using  the  Quantity 
One  software  (version  4.2;  Bio-Rad).  SBEM  and  mammaglobin-1  expression 
was  normalized  to  GAPDH  eTqpression  as  described  previously  (11),  Correla¬ 
tion  between  SBEM  expression  and  tumor  characteristics  or  mammaglobin-l 
expression  was  tested  by  calculation  of  the  Spearman  coefficient  r.  Differences 
between  tumor  subgroups  were  tested  using  the  Mann- Whitney  two-tailed  test 
or  Fisher’s  exact  test. 

In  Situ  Analysis  of  SBEM  mRNA  Expression.  In  situ  hybridization  was 
performed  on  adjacent  paraffin-embedded  breast  tumor  tissue  sections 
corresponding  to  a  case  shown  to  express  high  levels  of  SBEM  mRNA  by 
RT-PCR,  using  SBEM  ^'P-labcled  sense  and  antisense  probes,  as  described 
previously  (12). 


Results 

Identification  of  a  Putative  Novel  Breast-specifk  Gene.  The 
cDNA  xProfiler  tool  (see  “Materials  and  Methods”)  was  used  to 
identify  tissue-restricted  cDNAs  with  preferential  representation  in 
libraries  prepared  from  normal  breast  tissue  and  breast  tumors.  This 
search  identified  a  new  cluster  of  ESTs  now  grouped  under  the 
UniGene  identifier  number  Hs.348419.  Of  the  30  ESTs  found  in  this 
cluster,  15  are  ascribed  to  breast  cDNA  libraries,  9  were  isolated  from 
random  activation  of  gene  expression  or  pooled  tissues,  5  were  iso¬ 
lated  from  fetal  sources  (fetal  heart  and  fetal  skin),  and  1  came  from 
a  head  and  neck  tumor  cDNA  library.  Alignment  of  these  ESTs  led  to 
the  construction  of  a  500-bp  consensus  cDNA  sequence  containing  a 
90-amino  acid  open  reading  frame  in  which  the  initiating  methionine 
is  framed  by  a  nearly  perfect  consensus  motif  for  translation  initiation 
(5'-CCACCATGA-3';  Ref.  13).  Further  database  analysis  showed 
that  this  sequence,  interrupted  by  three  introns,  is  present  on  chromo¬ 
some  12ql3.2.  Primers  were  designed  to  span  the  open  reading  frame, 
and  we  cloned  a  396-bp  fragment  from  both  MCF-7  cells  and  breast 
tissue,  which  we  called  SBEM  (GenBank  accession  number 
AF4 14087).  The  presence  of  a  hydrophobic  signal  peptide  (residues 
1“19;  Ref  6)  within  the  protein  sequence  (GenBank  accession  num¬ 
ber  AAL02119)  suggests  that  SBEM  is  a  secreted  protein  subject  to 
proteolytic  processing.  The  NetOGlyc  glycosylation  algorithm  (7) 
further  predicts  this  protein  to  be  0-gIycosyIated  on  most  of  its  16 
threonine  residues.  The  SBEM  protein  contains  three  tandem  copies 
of  a  neutral  octapeptide  core  repeat  (ThrThrAIaAlaXxxThrThrAla, 
where  Xxx  corresponds  to  Ala,  Pro,  or  Ser).  The  NH2  and  COOH 
termini  of  the  processed  polypeptide  are  otherwise  charged  and  feirly 
polar.  These  features  suggest  strong  similarity  to  many  sialomucins, 
although  this  protein  lacks  a  transmembrane  domain  and  is  substan¬ 
tially  shorter  than  most  other  known  epithelial  mucins  (14,  15). 

Expression  of  SBEM  mRNA  Is  Restricted  to  the  Mammary  and 
Salivary  Glands.  Database  searches  suggested  that  SBEM  expres¬ 
sion  was  mainly  restricted  to  breast  tissue.  To  confirm  this  prediction, 
we  performed  hybridization  analysis  with  an  RNA  MasterBIot  con¬ 
taining  highly  purified  polyadenylated  RNA  from  43  adult  and  7  fetal 
human  tissues  arrayed  on  a  nylon  membrane.  A  SBEM  cDNA  probe 
hybridized  exclusively  to  mRNA  from  the  mammary  and  salivary 
glands  (Fig.  1).  Of  note,  no  expression  was  observed  in  colon,  lung, 
uterus,  ovary,  liver,  pancreas,  kidney,  or  prostate,  all  of  which  repre¬ 
sent  common  primary  tumor  sites.  Additionally,  no  hybridization  to 
any  of  the  fetal  RNAs  was  observed. 
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SBEM  mRNA  Is  Expressed  in  Breast  Cancer  Cell  Lines  but  not 
in  Cell  Lines  of  Non-Breast  Origin.  The  profile  of  SBEM  mRNA 
expression  was  ftirther  assessed  using  RT-PCR.  followed  by  PCR 
amplification,  in  a  panel  of  human  breast  and  non-breast  cell  lines.  A 
SBEM  PCR  product  of  the  expected  size  (396  bp)  was  readily  de¬ 
tected  in  MCF7  and  ZR-75-1  breast  tumor  cells  (data  not  shown). 
Lower  but  reproducible  expression  was  also  observed  in  primary 
HMECs  and  in  several  established  breast  epithelial  cell  lines  including 
T-47D,  M13SV-1  (8),  MCFIOATI,  and  MCF10AT3c  (16).  MDA 
MB-231  breast  tumor  cells  were  negative  for  SBEM  expression,  as 
were  six  tumor  cell  lines  of  non-breast  origin  (uterus:  RL95-2,  SK- 
UT-IB,  Hec  lA;  cervix:  HeLa;  prostate:  LNCaP;  and  liver  HepG2). 
As  controls,  we  also  examined  the  expression  of  a  housekeeping  gene 
{GAPDH)  and  mammaglobin-1,  an  established  mammaiy-specific 
gene  that  is  being  independently  investigated  as  a  promising  marker 
for  breast  tumor  diagnosis  and  nodal  metastasis  (12,  17).  Of  the  cell 
lines  tested,  only  HMEC  and  ZR-75-1  cells  expressed  mammaglo¬ 
bin-1,  consistent  with  published  reports. 

Analysis  of  SBEM  mRNA  in  Human  Breast  Tumors.  Northern 
blot  analyses  performed  on  a  small  series  of  10  cases  revealed  that 
SBEM  mRNA  was  600  bp  long  and  differentially  expressed  from  one 
sair5)le  to  another  (data  not  shown).  To  determine  whether  SBEM 
mRNA  was  widely  expressed  in  human  breast  tumor  tissue,  54  human 
breast  tumors,  spanning  many  ER  and  PR  levels  as  well  as  tumor 
grade  and  nodal  status,  were  selected  from  the  Manitoba  Breast 
Tumor  Bank.  Total  RNA  was  extracted  fi-om  frozen  tissue  sections 
and  reverse  transcribed.  PCR  aiiq)lification  of  GAPDH  (control), 
mammaglobin^l  and  SBEM  cDNA  was  then  performed.  A  PCR 
product,  396-bp  long  was  detected  in  all  but  three  tumors  (data  not 
shown)  when  using  SBEM-specific  primers.  After  cloning  and  se¬ 
quencing,  this  product  was  shown  to  correspond  to  SBEM  cDNA. 
Quantification  of  the  SBEM  signal  relative  to  the  GAPDH  signal  was 
performed  as  described  in  “Materials  and  Methods.”  No  correlation 
was  found  between  SBEM  expression  and  tumor  characteristics  such 
as  ER  (ff  -  54;  Spearman  r  =  -0.01,  P  =  0.89)  and  PR  {n  =  54; 
Spearman  r  =  -0.03,  P  =  0.77)  levels  or  tumor  grade  («  =  44; 
Spearman  r  ==  -0.06,  P  =  0.68).  Interestingly,  however,  the  SBEM 
signal  correlated  positively  with  mammaglobM  expression  («  =  54; 
Spearman  r  =  0.340,  P  =  0.01 1).  Subgroup  comparison  of^PPA/and 
mammaglobin-1  expression  confirmed  our  previous  observation®  that 
mammaglobin-1  expression  is  higher  in  ER-positive  and  low-grade 
tumors  (Table  1).  Interestingly,  although  not  statistically  significant 
(P  ==  0.09),  higher  SBEM  expression  was  found  in  lyn^h  node¬ 
positive  compared  with  node-negative  tumors.  Also  of  interest  is  the 

®  E.  Leygue,  L.  C.  Murphy,  and  P.  H.  Watson,  unpublished  results. 

Table  1  Median  values  (arbitrary  units)  of  SBEM  agfression,  mammaglobin-l  (Mam) 
expression,  and  SBEM:mammaglobm-l  ratio  in  different  tumor  subgroups 


Node- 


Nodc  + 


ER  +  34  2.29  0.66  0.63  0.04  3.39  0.22 

ER~  20  2.15  0.34  4.44 

PR+  34  2.30  0.63  0.59  0.27  3!56  0  50 

PR-  20  1,97  0.38  3.58 

Node  +  36  2,30  0.09  0.45  0.62  3.97  0.04 

Node  “  10  1.96  0.70  2.42 

Grades  5-7  24  2.05  0.84  0.59  0.03  3.33  0.10 

Grades  8-9  20  2.24 _ 0.34  4.02 

ER  +,  >3  fnK)l/ing  of  protein;  ER  £3  fmol/mg  of  protein;  PR  +,  >10  ftnol/mg 
of  protein:  PR  ^sio  fmol/nrg  of  protein,  as  determined  by  ligand  binding  assay.  Node 
+,  confirmed  metastasis  in  axillary  lymph  nodes;  Node  absence  of  metastasis  in 
analyzed  a.xinary  lymph  nodes.  Grade,  Nottingham  grading  system.  Subgroup  comparison 
was  performed  using  the  Maon>Whitney  two-tailed  test 

Bold  face  P  values  correspond  to  statistically  significant  differences  between  tumor 
subgroups  (P  <  0.05). 


Fig.  2.  RT-PCR  analysis  of  mammaghbin-I.  SBEM,  and  GAPDH  mRNA  expression 
m  primary  breast  tumors  (P)  and  their  corre^onding  axillary  \ytrph  nodes  (JV),  histolog- 
i^y  shown  to  contain  (Node  +)  or  not  to  contain  (Node  -)  metastases.  Mammaglobin-1 
(MAM)  PCR  products  were  tun  separately  (A,  dotted  arrow),  whereas  SBEM  md  GAPDH 
PCR  products  were  mb^d  before  separation  on  2%  agarose  gels  prestained  with  ethidium 
bromide  (P).  Gray  arrow,  product  corresponding  to  SBEM;  black  arrow,  product  cone- 
ending  to  GAPDH.  M,  molecular  weight  markers  (d>ptl74  RF  DNA/Hae  m  fragments; 
Life  Technologies,  Inc.,  Grand  Island,  NY).  N.  negative  control,  no  cDNA  added  during 
the  PCR  reaction. 


fact  that  the  SBEM:mammaglobin-l  ratio  is  significantly  {n  =  46; 
Mann-Whitney,  P  =  0.04)  higher  in  these  lymph  node-positive 
tumors. 

SBEM  mRNA  Expression  in  Primaiy  Breast  Tumors  and  Their 
Corresponding  Axillary  Nodes.  We  next  investigated  the  possibility 
that  SBEM  mRNA  could  be  a  tissue  marker  of  axillaiy  lyn^h  node 
metastasis.  Twenty  independent  cases  were  selected,  including  14 
tumors  that  were  axillary  lynqjh  node  positive  and  6  that  were  node 
negative.  Total  RNA  was  extracted  from  frozen  primary  tumor  sec¬ 
tions  and  frozen  node  sections  of  corresponding  axillary  lynqjh  nodes. 
The  histological  status  of  all  tissues  was  confirmed  in  paraffin  sec¬ 
tions  cut  fi-om  adjacent  mirror  image  paraffin  tissue  blocks  that  had 
been  processed  in  parallel  to  the  frozen  blocks.  These  RNAs  were 
reverse-transcribed  and  analyzed  by  RT-PCR  using  5»£A/-specific 
primers.  PCR  was  performed  three  times,  giving  the  same  result.  A 
representative  experiment  is  shown  in  Fig.  2.  A  signal  corresponding 
to  SBEM  was  detected  in  all  lyn^h  nodes  containing  metastatic  cells 
by  histopathological  assessment  (14  of  14  cases).  In  contrast,  no 
signal  was  detectable  in  lymph  nodes  from  cases  without  histologi¬ 
cally  detectable  tumor  cells  (0  of  6  cases).  RT-PCR  detection  of 
SBEM  mRNA  in  axillary  lymph  nodes  is  therefore  strongly  associated 
(Fisher’s  exact  test,  P  <  0.001)  with  the  histopathological  detection  of 
lymph  node  metastases.  The  higher  sensitivity  afforded  by  RT-PCR 
detection  therefore  indicates  that  SBEM,  perhaps  together  with  mam¬ 
maglobin-1,  represents  an  excellent  marker  for  the  detection  or  con¬ 
firmation  of  occult  breast  tumor  metastasis,  where  histopathology  may 
not  be  definitive. 

SBEM  mRNA  Is  Expressed  in  Mammary  Epithelial  Tumor 
Cells  in  Vivo,  To  further  establish  whether  SBEM  was  expressed  by 
mammaiy  epithelial  cells  in  vivo,  paraffin  breast  tumor  tissue  sections 
corresponding  to  a  case  shown  to  strongly  express  SBEMniRHA  by 
RT-PCR  were  studied  by  in  situ  hybridization.  No  signal  was  detect¬ 
able  when  using  a  sense  probe  (Fig.  3A).  In  contrast,  a  signal  was 
observed  in  epithelial  tumor  cells  when  using  an  antisense  probe 
(Fig.  3R).  SBEM  mRNA  was  not  detected  in  stromal  or  inflammatory 
cells  in  any  of  the  sections  studied. 

Discussion 

This  article  reports  the  identification,  cloning,  and  preliminary 
characterization  of  a  cDNA  encoding  a  novel  mucin-like  protein  that 
displays  an  unusually  narrow  pattern  of  expression.  Hybridization 
analysis  revealed  that  SBEM  mRNA  was  only  detectable  in  two 
normal  tissues,  breast  and  salivary  gland.  Interestingly,  the  tissue- 
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Fig.  3.  Expression  of  SBEM  mRNA  in  a  primaiy 
breast  tumor  studied  by  in  situ  hybridization.  These 
plates  illustrate  consecutive  sections  fiom  a  single 
breast  tumor  and  show  H&E>stained  paraffin  sec¬ 
tion  treated  wiUi  a  sense  probe  (A)  and  SBEM 
mRNA  expression  in  epitl^lial  cells  detected  using 
an  ai^sense  probe  (B).  Mauve  and  Uack  arrow¬ 
heads  ^ow  tumor  epithelial  cells  and  regions 
of  stroma  with  inflammation,  respectively.  Bar, 
30  pm. 


specific  expression  that  we  observed  experimentally  directly  reflects 
the  distribution  of  ESTs  within  the  Hs.348419  cluster.  Indeed,  as 
mentioned  above,  only  two  adult  tissues  (breast  and  head/neck  tu¬ 
mors)  have  been  shown  to  express  SBEM-related  ESTs,  The  fact  that 
SBEM  is  also  expressed  in  salivaiy  tissue  does  not  undermine  the 
possible  use  of  SBEM  as  a  marker  of  breast  cancer,  because  tumors  of 
the  salivary  gland  are  less  common  and  can  readily  be  distinguished 
clinically. 

Among  the  primary  breast  tumors  examined  in  this  study  (repre¬ 
senting  mostly  invasive  ductal  carcinoma),  SBEM  mRNA  was  ob¬ 
served  by  RT-PCR  analysis  in  the  majority  (>90%)  of  cases.  Despite 
a  significant  overall  correlation  between  the  expression  of  SBEM  and 
mammaglobin-1  mRNA,  a  significantly  higher  SBEM:mammaglo- 
bin-1  ratio  was  observed  in  primary  tumors  associated  with  positive 
axillaiy  lyn^)h  nodes  as  con^ared  with  node-negative  tumors.  This 
was  mostly  attributable  to  a  trend  toward  higher  SBEM  expression  in 
node-positive  tumors.  Although  further  analysis  of  a  larger  number  of 
tumors  will  be  required  to  confirm  these  observations,  this  may 
suggest  differences  in  the  biology  of  these  tumors  and  also  a  possible 
role  of  SBEM  and  mammaglobin-l  in  the  mechanisms  involved  in 
tumor  metastasis.  Our  findings  indicate,  however,  that  SBEM  expres¬ 
sion  is  a  common  feature  of  breast  cancer  and  can  furthermore  serve 
as  a  useful  marker  for  breast  nodal  metastasis,  both  for  detection  of 
micrometastatic  cells  within  lymph  nodes  as  well  as  in  the  differential 
diagnosis  of  the  primary  origin  of  an  unknown  metastasis.  This 
potential  is  enhanced  by  the  conserved  SBEM  expression  in  high 
grade  and  ER/PR-negative  tumors  that  are  most  likely  to  metastasize. 

The  potential  diagnostic  relevance  of  SBEM  is  also  increased  by  its 
predict^  biochemical  structure.  The  SBEM  cDNA  sequence  codes  for 
a  90-amino  acid  polypeptide  that  contains  a  distinctive  tandem  repeat, 
rich  in  alanine  and  threonine  residues,  that  represents  a  probable  target 
for  D-glycosylation.  Consistent  with  such  posttranslational  modifica¬ 
tion  is  the  presence  of  a  well-defined  signal  peptide,  leading  us  to 
predict  that  SBEM  is  likely  to  be  processed  at  the  apical  surface  of 
luminal  epitiielial  cells  and  to  be  secreted  into  the  alveolar  or  ductal 
lumen.  Further  study  is  needed  to  ascertain  whether  higher  SBEM 
expression  occurs  in  association  with  tumors. 

Secreted  (or  transmembrane)  proteins  that  contain  internally  re¬ 
peated,  densely  glycosylated  neutral  core  motifs  such  as  this  are 
characteristic  of  mucins,  which  are  typically  expressed  by  the  surface 
q)ithelium  of  secretory  mucosae  and  by  exocrine  glands  (14, 15).  The 
role  of  mucins  is  primarily  one  of  hydrating  and  lubricating  epithelial 
linings,  although  several  mucins  have  been  implicated  in  modulating 
both  cell  adhesion  and  growth  factor  signaling  (18, 19).  Furthermore, 
mucins  have  a  well-established  link  to  cancer,  best  illustrated  by  the 
product  of  the  MUCl  gene.  MUCl  is  overexpressed  in  a  variety  of 


epithelial  tumors  including  breast  cancer  and  gives  rise  to  several 
well-characterized  tumor  antigens  including  CA15.3  and  CA27.29 
(20,  21).  Combined  with  the  loss  of  cell  polarity  and  changes  in 
glycosylation  patterns  observed  in  transformed  epithelial  cells,  over¬ 
expression  of  MUCl  results  in  the  appearance  of  mucin-derived  tumor 
antigens  in  the  sera  of  cancer  patients  that  are  not  seen  in  normal 
controls  (14, 22).  We  hypothesize  that  a  similar  situation  may  hold  for 
SBEM  expression  in  the  human  mammary  epithelium  and  in  human 
breast  tumors.  However,  MUCl  displays  relatively  broad  expression 
among  epithelial  tissues  including  the  colon,  breast,  pancreas,  ovary, 
prostate,  tracheobronchial  tree,  stomach,  and  uterus.  For  this  reason, 
A/l/Ci-derived  tumor  antigens  have  relatively  poor  specificity  for 
individual  tumor  types,  and  their  clinical  utility  is  limited  to  monitor¬ 
ing  the  efficacy  of  cancer  therapy  and  warning  of  tumor  relapse  or 
malignant  spread  (21-23). 

Parallels  between  SBEM  and  known  epithelial  mucins  such  as 
MUCl,  together  with  its  more  narrowly  restricted  pattern  of  expres¬ 
sion,  suggest  that  this  novel  gene  represents  an  attractive  candidate  for 
a  breast  biomarker  with  potential  for  cancer  diagnostics,  as  well  as 
being  a  possible  future  target  for  the  development  of  a  breast  tumor 
vaccine.  Moreover,  the  absence  of  SBEM  expression  in  normal  lyn^h 
node  tissue  suggests  that  this  gene  could  also  be  used  to  detect  breast 
micrometastases  in  axillary  lyn^h  nodes. 
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Abstract. 

Estrogen  receptors  (ERa  and  ERP)  are  clearly  multifaceted  in  terms  of  structure 
and  function.  Several  relatively  abundant  ERP  isoforms  have  been  identified,  which 
can  be  differentially  expressed  in  various  tissues.  In  order  to  provide  insight  into  the 
possible  role  of  the  ERP  family  in  breast  tissue  a  study  of  the  putative  functions  of  the 
hERpi,  hERP2  and  hERpS  isoforms  was  undertaken.  Only  hERpi  was  found  to  bind 
ligand,  which  induced  conformational  changes  as  determined  by  protease  digestion 
assays.  All  ERP  isoforms  could  bind  to  and  bend  DNA  although  the  relative 
efficiency  with  which  they  bound  DNA  differed  with  hERa  > 
hERpi>hERP2  »hERP5  .  All  ERP  isoforms  inhibited  ERa  transcriptional  activity 
on  an  ERE-reporter  gene.  The  relative  activities  were  hERpi>hERP2> 
hERpS,  however  only  hERpi  had  transcripitonal  activity  alone.  Both  LY117018- 
hERa  and  LY117018-hERpi  complexes  alone  could  activate  transcription  on  a  TGF- 
PS-  reporter  gene.  Although  hERP2  and  hERps  had  no  activity  alone,  they  inhibited 
ERa  but  not  hERpi  transcriptional  activity  of  TGF-PS-CAT.  In  marked  contrast  to 
activity  on  an  ERE-CAT  reporter  gene,  hERpi  did  not  modulate  ERa  transcriptional 
activity  on  a  TGF-P3-CAT  reporter  gene.  These  data  support  promoter-specific 
differential  activities  of  ERP  isoforms  with  respect  to  models  of  ERa  regulated  gene 
expression,  and  suggest  that  they  may  have  a  role  in  differentially  modulating 
estrogen  action. 


Introduction. 

The  estrogen  receptor  family  of  steroid  hormone  receptors  is  clearly  multifaceted 
(1)  and  more  complex  than  originally  thought.  There  are  two  genes  which  encode 
estrogen  receptors  (ER),  ERaand  ERp.Both  are  ligand  regulated  transcription 
factors  which  classically  modulate  target  gene  transcription  by  binding  as  homo- 
and/or  hetero-dimers  to  estrogen  responsive  sequences  in  target  gene  promoters  (2). 
These  receptors  likely  have  distinct  roles  in  estrogen  action,  independent  of  each 
other  (3)  when  they  are  expressed  separately  but  can  also  have  direct  interactions  due 
to  heterodimerization  when  the  receptors  are  expressed  together  in  the  same  target 
cell  (4).  In  addition  both  ERs  may  encode  variant  isoforms  generated  by  alternative 
splicing  mechanisms  (5,  6).  In  particular  there  are  data  to  support  variant  isoforms  of 
ERP  at  the  protein  level  (7,  8)  Furthermore,  we  have  shown  in  human  breast  tissues 
that  variant  forms  of  ERP  are  more  abundant  than  the  wild-type  at  least  at  the  RNA 
level  (9). 

Human  ER32  (also  called  hERpbx  (10))  and  ERPS  variant  mRNAs  are  missing 
the  wild  type  exon  8  sequences  and  contain  extra  sequences  which  are  distinct  from 
each  other,  followed  by  sequences  that  are  then  identical  with  each  other  (see  Fig  1). 
They  are  predicted  to  encode  C-terminally  truncated  ERpiike  proteins  identical  to 
wild-type  until  amino  acid  residue  468  (by  reference  to  the  long  form  of 
hERpl  (11).  After  amino  acid  468  hERp2is  predicted  to  encode  28  novel  amino 
acids,  with  the  full-length  protein  having  a  predicted  molecular  mass  of  55.5  kDa.  In 
contrast,  after  amino  acid  468  hERP5  is  predicted  to  encode  only  5  novel  amino  acids 
with  the  full-length  protein  having  a  predicted  molecular  mass  of  53  kDa. 

Although  total  ERP  expression  appears  to  decrease  between  normal  breast  and 
ER-i-  breast  tumors  (12,  13),  the  relative  expression  of  the  variant  ERP  isoforms  to 
the  wild-type  ERP  can  also  change  during  breast  tumorigenesis,  at  least  at  the  RNA 
level  (9).  This  suggests  that  the  expression  and/or  the  activity  of  the  ERP  family  of 
receptors  changes  during  breast  tumorigenesis  and  may  have  a  role  in  this  process  as 
well  as  have  a  role  in  altered  estrogen  action  that  occurs  during  breast  tumorigenesis. 
In  order  to  provide  insight  into  the  possible  role  of  the  ERP  family  in  breast  tissue  we 
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have  undertaken  a  study  of  the  putative  functions  of  the  hERpi,hERP2and 
hERpS  isoforms. 

MATERIALS  AND  METHODS. 

Materials.  17P-estradiol  (E2),  4-hydroxytamoxifen  (4-OH-TAM)  and  CAPS  (3- 
cyclohexylamino-l-propanesulfonic  acid)  were  from  Sigma  Chemical  Co  (St  Louis, 
MO).  ICI  182,780  was  a  gift  from  Dr  AE  Wakeling  (AstraZeneca  Pharmaceuticals, 
England).  LYl  17018  was  a  gift  from  Eli  Lilly  Co.  (Indianapolis,  IN).  (^H)-17P- 
estradiol  (^H-Ej),  ('‘’C)-chloramphenicol,  and  (”S)-methionine  were  from  New 
England  Nuclear  (Boston,  MA).  (a-'^P)-dCTP,  (y-''P)-ATP,  and  (''S)-ATP  were 
from  ICN  Pharmaceuticals  (Irvine,  CA).  All  cell  culture  reagents  were  obtained  from 
GIBCO/BRL  (Burlington,  Ontario). 

In  vitro  transcription  and  translation.  In  vitro  transcription/translation  reactions 
were  performed  using  a  coupled  transcription/translation  system  (TnT  coupled 
Reticulocyte  Lysate  System,  Promega,  Madison,  WI).  Reactions  were  performed 
according  to  the  manufacturer's  instructions. 

Scatchard  Analysis.  Human  ERa  (pcDNA3.1/wild-type  human  ERa  from  HEGO, 
(14)),  human  ERpl  (pcDNA3.1hERpl,  long  form  of  530  amino  acids  (11,  15)), 
human  ERp2  (pcDNA3.1  hERP2,  long  form  (10))  and  human  ERP5  (pcDNA3.1 
hERP5,  long  form)  proteins  were  synthesized  by  in  vitro  transcription-translation  as 
described  above.  Ligand  binding  studies  were  conducted  as  previously  described 
(16).  In  vitro  generated  receptor  was  diluted  10-fold  in  buffer  (10  mM  Tris-Cl,  pH 
7.5,  1.5  mM  EDTA,  10  mg/ml  BSA,  10%  glycerol)  and  kept  on  ice  until  use.  One 
hundred  pi  of  the  diluted  protein  were  used  in  each  binding  reaction,  that  contained 
varying  concentrations  of  (^Hl-Ej  (0.01  -  100  nM),  followed  by  overnight  incubation 
at  4°C.  Non-specific  binding  was  determined  by  parallel  incubations  containing  a 
200-fold  excess  of  unlabelled  E2.  Unbound  steroid  was  removed  by  addition  of  500  pi 
of  0.5%  charcoal-0.05%  dextran  in  the  above  dilution  buffer  for  30  mins  at  4°C 
followed  by  centrifugation  at  10000  x  g  for  10  mins  at  4“C.  Radioactivity  was 
determined  in  an  aliquot  of  the  supernatant  and  in  aliquots  of  total  (^H)-E2  solutions 
using  a  scintillation  counter.  The  ratio  of  specifically  bound/unbound  steroid  and  the 
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concentration  of  specifically  bound  steroid  were  used  for  Scatchard  analysis,  from 
which  was  determined  the  equilibrium  dissociation  constant,  K^. 

Limited  proteolytic  digestion  analysis.  Conformational  studies  were  performed  as 
described  previously  (17)./«  vitro  synthesized  ERs  were  incubated  with  agonist  (E2, 
DES)  and  antagonists  (TAM,  LY)  for  overnight  at  4®  C.  The  liganded  receptors  were 
then  diluted  1:  10  (v/v)  in  TE  buffer,  then  20  ul  of  these  ER  solution  was  treated  with 
increased  concentration  of  trypsin  (0.2  to  5  ug)  for  20  min  at  room  temperature  and 
stopped  loading  buffer  was  added.  The  samples  were  boiled  and  were  directly 
analyzed  by  SDS -polyacrylamide  gel  electrophoresis  (10%  w/v).  The  gel  was  dried 
and  the  digested  bands  were  visualized  by  autoradiography. 

Electrophoretic  Mobility  Gel-Shift  Assay  (EMSA).  In  vitro  synthesized  human 
ERs  were  used  for  EMSA.  Typically  1  |xl  of  programmed  lysates  containing  equal 
amounts  of  each  receptor  as  determined  by  polyacrylamide  gel  electrophoresis  of  ^^S- 
methionine  labelled  protein  generated  in  parallel  in  vitro  transcription-translation 
assays  was  assayed  in  EMSA.  One  |xl  of  lysate  was  incubated  in  a  final  volume  of  20 
pi,  and  the  reaction  solution  was  5  mM  Tris-Cl,  pH  7.5,  100  mM  KCl,  2  mM  DTT, 
5%  v/v  glycerol  and  contained  2  mg  polyd(I-C).  The  binding  reaction  was  initiated 
by  adding  1  pi  (approx.  10  fmol)  of  5'-(^^P)-end-labeled,  double  stranded  ERE 
oligonucleotide  (35mer,  5'-AACTTTGATCAGGTCACTGTGACCTGACTTTGGAC 
-3’  containing  the  vitellogenin  A2  ERE  sequence),  and  the  mix  was  incubated  at  20‘’C 
for  30  min.  DNA-bound  complexes  were  electrophoretically  separated  on  a  4.5% 
loosely  cross-linked  acrylamide  gel  (1:29  bisracrylamide)  at  150V  for  90  min  at  20®C 
temperature  in  0.5xTBE  buffer.  Gels  were  then  vacuum  dried  and  autoradiographed. 
To  identify  immunoreactive  ER  within  retarded  DNA-bound  complexes,  parallel 
incubations  containing  Ipg  of  ER  antibody  (usually  1  pi  of  H222  for  ERa,  or  1  pi  of 
PAI-310  for  ERpi  and  ERP2)  were  run  to  determine  the  presence  of  super-shifted 
antibody-bound  ER-ERE  complexes  (data  not  shown).  Reticulocyte  lysates 
containing  in  vitro  translated  ER  proteins  were  incubated  with  or  without  saturating 
concentrations  of  ligand  (estrogen  or  antiestrogen)  at  4®C  overnight  to  allow  receptors 
to  bind  ligand,  followed  by  EMSA. 
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DNA  Bending  Assay.  The  DNA  bending  vector  ERE  Bend  I  (kindly  provided  by  Dr 
A  Nardulli,  University  of  Illinois,  Urbana,  IL)  (18)  was  digested  with  EcoRI  and 
EcoRV  to  produce  a  430  bp  DNA  fragment  with  a  single  consensus  ERE  either  at  the 
end  (EcoRI  fragment)  or  in  the  middle  (EcoRV  fragment).  The  fragments  were  then 
gel  purified,  labeled  by  incubation  with  polynucleotide  kinase  in  the  presence  of  (y- 
^  P)-ATP,  and  purified  on  a  G50  Sephadex  column.  Gel  mobility  shift  assays  were 
earned  out  essentially  as  described  above.  Aliquots  of  the  binding  reactions  were  run 
on  8%  non-denaturing  acrylamide  gels,  dried  and  exposed  to  X-ray  film.  The  degree 
of  DNA-bending  was  determined  using  the  method  of  Thompson  and  Landy  (19). 
Cells,  Cell  Culture  and  Transient  Transfection.  For  transient  transfection  analysis 
tagged  ER  expression  vectors  were  generated.  Human  ERa,  ER(3l,  ERp2  and  ERp5 
were  tagged  at  their  N-terminus  with  a  polyhistidine  and  an  Xpress'^“  epitope  tag 
using  the  pcDNA4/HisMax(A)  plasmid  (Invitrogen  Canada  Inc,  Burlington,  ON). 
Cos-1  and  Cos-7  cells  were  obtained  from  the  ATCC  (Manassas, V A).  The  cells  were 
routinely  cultured  in  DMEM  containing  5%  v/v  fetal  calf  serum  (FBS),  1%  w/v 
glucose,  glutamine  and  penicillin-streptomycin  (5%CM).  To  obtain  estrogen  depleted 
cells,  the  culture  medium  of  stock  cells  was  changed  to  phenol  red-free  DMEM 
supplemented  with  5%  charcoal  dextran-treated  FBS,  1%  w/v  glucose,  glutamine  and 
penicillin-streptomycin  (5%CS)  and  replaced  every  2  days.  Six  days  later  the  medium 
was  replaced  by  10%CS  until  required  for  experiments.  For  transient  transfection 
experiments,  the  cells  were  transfected  using  the  Effectene  transfection  reagent 
according  to  the  manufacturer's  instruction  (QIAGEN,  Mississauga,  Ontario).  Briefly, 
the  day  before  transfection,  the  estrogen  depleted  cells  were  seeded  in  six-well  plates 
at  2.5  X  10^  cells  per  well  in  2  ml  of  5%CS  and  left  overnight.  The  plates  were  70- 
80%  confluent  on  the  day  of  transfection.  The  transfection  mixture  was  prepared 
according  to  the  manufacturer's  protocol,  then  fresh  medium  (5%CS)  was  added  to 
the  transfection  mixture  and  0.6  ml  per  well  of  the  above  mixture  either  ERE-II-TCO- 
CAT  (gift  from  P  Webb,  (20))  or  TGF-|33-CAT-reporter  plasmid  DNA  (21).  ER 
expression  plasmid  (50-450  ng)  or  empty  vector  and  100  ng  P-gal  pCHllO  plasmid 
DNA  (Pharmacia  canda,  Mississauga,  ON)  was  added  drop-wise  into  the  medium 
(CS)  and  the  plates  were  gently  swirled  to  ensure  uniform  distribution  of  the  DNA- 


Effectene  complexes.  Vehicle  (ethanol),  estradiol- 1 7 13  or  LYl  17018  was  then  added 
20  to  30  minutes  later.  The  cells  were  left  for  48  hr  and  then  harvested.  Cell  extracts 
were  prepared  by  freeze/thawing  and  used  to  determine  CAT  and  P-galactosidase 
activity  as  previously  described  (22) 

Western  Blot  Analysis.  For  Western  blot  analysis,  2.5  x  10^  Cos-1  cells  were  set  up 
in  6  well  plates,  then  transiently  transfected  with  plasmids  and  treated  with  estrogen 
or  antiesotrgen  under  the  same  conditions  as  for  the  CAT  assay  described  above. 
Cells  were  harvested  48  hour  after  transfection,  washed  once  with  Isoton  II  and  then 
the  washed  cell  pellets  were  resuspended  in  200  pi  Isoton  II.  Aliquots  of  cell 
suspension  (150  pi)  were  extracted  and  used  for  Western  blots  and  the  remainder  was 
used  for  determination  of  P-galactosidase  activity.  For  Western  blotting,  the  cells 
were  pelleted  and  then  extracted  using  40  pi  of  hot  (95°C)  extraction  buffer  J  with 
shaking  for  20  min  at  95°C  as  previously  described  by  Joel  et  al,  (23).  The  entire 
extract  was  subjected  to  10%  SDS-polyacrylamide  gel  electrophoresis  as  previously 
described  (24).  The  separated  proteins  were  transferred  to  nitrocellulose  membranes 
and  processed  as  previously  described  (24).  Detection  of  the  tagged  estrogen  receptor 
proteins  was  by  incubation  of  blots  with  anti-Xpress  antibody  (1/5000  in  TEST, 
Cat#R9 10-25,  Invitrogen  Canada  Inc,  Burlington,  ON)  overnight  at  4°C,  followed  by 
washing  and  incubation  with  secondary  antibody  (horseradish  peroxidase  conjugated 
goat  anti-mouse  antibody,  1/5000  in  TEST,  Jackson  Immuno  Research  Labs  Inc,  PA) 
at  room  temperature  for  2  hours.  Visualization  was  carried  out  using  the  SuperSignal 
West  Dura  Extended  Duration  Substrate  kit  (Pierce,  IL)  according  to  the 
manufacturer’s  instructions. 

Statistical  Analysis.  Differences  in  variance  were  tested  using  ANOVA  where 
appropriate.  Differences  between  individual  mean  values  were  then  determined  using 
students  t-tests.  All  tests  were  performed  using  GraphPad  Prism  statistical  analysis 
software  (GraphPad  Software  Inc). 

Results. 

Identiflcation  and  Organization  of  hER^l,  hERP2  and  hER^S  cDNA  Sequences 
at  the  hER^  Locus  on  Chromosome  14. 
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The  estrogen  beta  gene  has  been  localized  to  human  chromosome  14q22-24  and  the 
genomic  structure  of  the  8  exons  comprising  hERpi  has  been  previously  published 
(4).  However,  the  previously  described  variant  hERp2  (also  called  hER|3cx)  and 
hERps  mRNA  only  contain  sequences  corresponding  to  exons  1  to  7  of  hERpi  and 
then  diverge  (see  Figure  lA).  They  do  not  contain  exon  8  sequences  of  hERpi,  but 
contain  sequences  termed  exon  9,  which  are  located  downstream  of  exon  8  on 
chromosome  14  (Figure  IB),  identified  using  database  sequences  of  chromosome  14 
(Accession  numbers  CNSOIRHJ  and  AF215937)  and  the  Human  Genome  Working 
Draft.  It  should  be  noted  that  neither  of  these  genomic  sequences  contain  an  extra  A 
5  of  the  start  site  of  translation  for  tiERP,  that  would  place  another  upstream  ATG  in 
frame  with  the  known  coding  region  and  introduce  18  amino  acids  to  the  N-terminal 
of  the  known  coding  region  as  recently  described  (25).  Interestingly,  hERp5  mRNA 
also  contains  sequences  between  exon  7  and  part  of  exon  9  which  are  not  present  in 
either  hERPi  or  hERP2  mRNA.  These  hERp5  mRNA  specific  sequences  can  be 
found  immediately  following  exon  7  sequences  in  intron  7  of  the  human  ERp  gene 
(Figure  IB),  suggesting  that  the  normal  splice  donor  site  is  not  recognized  and  a 
cryptic  splice  donor  site  is  present  in  intron  7.  Furthermore  the  exon  9  sequences 
present  in  hERP5  cDNA  start  28  nucleotides  downstream  of  those  present  in  hERP2, 
suggesting  a  cryptic  splice  acceptor  site  is  present  within  exon  9.  There  are  also 
multiple  non-coding  exons  5’  to  exon  1  as  previously  identified  (4),  since  several 
hERP  cDNAs  contain  sequences  in  their  5’  UTR  which  are  found  further  upstream  of 
exon  1  on  chromosome  14  (Figure  IB,  IH-IC  seen  in  AB006589,  IB  seen  in 
NM_001437,  AX234658,  AF05428,  AF060555,  AB006589  references).  The 
sequences  of  a  recently  characterized  promoter  region  of  hERP  (26)  are  found 
immediately  upstream  and  overlapping  with  exon  1®.  However,  the  presence  of  hERP 
mRNAs  whose  5’UTR  contain  exonic  sequences  found  upstream  of  this  documented 
promoter  suggest  that  there  are  alternative  promoters  for  the  hERP  gene.  This  is 
similar  to  the  hERa  gene  and  suggests  that  regulation  of  expression  of  these  genes  is 
complex  (27). 

The  predicted  open  reading  frames  for  hERpi,  hERP2  and  hERps  are  shown  in 
Figure  2.  hERp2  contains  amino  acids  1-468  which  are  identical  to  hERpi  and  then 
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diverges  with  another  28  novel  amino  acids  encoded  in  the  open  reading  frame. 
hERp5  was  isolated  as  a  partial  cDNA  but  is  likely  to  also  be  identical  to  hERpl  from 
amino  acids  1-468  and  then  diverges  containing  another  5  novel  amino  acids.  Both 
these  variant  hER|3  proteins  would  be  truncated  at  the  C-terminus,  disrupted  in  helix 
1 1  and  missing  helix  12  and  therefore  unlikely  to  bind  ligand  or  have  AF2  mediated 
transcriptional  activity.  Lack  of  ligand  binding  has  been  confirmed  using  in  vitro 
generation  of  these  proteins  as  outlined  below. 

Ligand  Binding  Activity  of  hER^l  and  variant  isoforms  hERp2  and  hERps 
Proteins. 

Human  ERpi  has  been  shown  previously  (4)  to  bind  E2  with  high  affinity  and 
specificity,  and  our  data  confirm  these  findings.  Figure  3A  shows  specific  saturable 
binding  of  (^H)E2  to  in  vitro  translated  hERpi  with  a  calculated  =  0.11  nM. 
However,  the  open  reading  frames  of  hERP2  and  hERPS  cDNA  predict  for  proteins 
C-terminally  truncated  compared  to  hERpi  and  predicted  not  to  bind  ligand.  As 
shown  in  Figure  3  B  and  C  no  saturable  binding  of  (^H)E2  to  in  vitro  translated 
hERP2  or  5  was  observed. 

Human  Estrogen  Receptor  Isoform  Conformational  Status  and  Ligand  Induced 
Changes. 

To  determine  the  possible  conformational  status  of  variant  hERP  isoforms,  a 
previously  used  limited  trypsin  digestion  assay  (17, 28)  was  employed  to  compare  the 
proteolytic  digestion  patterns  of  variant  S^Lmethionine  labelled  hERa,  hERpi, 
hERP2  and  hERpS  in  the  presence  and  absence  of  estrogens  and  antiestrogens  (4-OH- 
tamoxifen  and  LYl  17018).  The  results  are  shown  in  Figure  4.  In  the  absence  of  any 
ligand  all  ER  isoforms  were  sensitive  to  proteolysis,  with  hERa  and  hERpi  being 
relatively  more  resistant  than  hERP2  and  hERP5.  In  the  presence  of  estradiol  both 
hERa  and  hERpi  become  more  resistant  to  digestion,  and  a  32.5  kDa  resistant  band 
(shown  by  *  in  Fig  4)  was  observed.  Addition  of  antiestrogens  4-OH-tam  and 
LYl  17018  increased  slightly  the  sensitivity  of  hERa  and  hERpi  to  trypsin  digestion 
compared  to  the  receptors  in  the  absence  of  ligand.  The  sensitivity  of  the  variant 
isoforms  hERP2  and  hERps  was  not  effected  by  ligand,  consistent  with  their  inability 
to  bind  ligand.  Furthermore,  the  variant  isoforms  were  more  sensitive  to  trypsin 


digestion  than  either  the  unliganded  hERa  or  hER|3l,  suggesting  that  the  variant 
isoforms  are  unlikely  to  be  in  an  activated  conformation. 

DNA  Binding  and  Bending  Activity  of  hER^l  and  variant  isoforms  hERP2  and 
hERps  Proteins. 

As  previously  demonstrated  hERpiand  hER32  can  bind  to  an  ERE  in  a  gel  mobility 
shift  assay  (Figure  5),  although  the  efficiency  of  hER32  DNA  binding  was  less  than 
hER3l(6, 10).  hER35  also  has  the  ability  to  bind  an  ERE  in  gel  mobility  shift  assays 
(Figure  5),  but  was  less  efficient  than  hER32.  The  specificity  of  the  binding  was 
determined  by  competition  with  excess  unlabelled  ERE  whereas  no  competition  was 
seen  with  an  excess  of  unlabelled  nonspecific  33  mer  oligonucleotide. 

DNA  bending  assays  demonstrated  that  hERa,  hERpi  and  hER32  were  all  able  to 
bend  DNA  as  demonstrated  by  the  reduced  mobility  of  complexes  when  the  ERE  is  in 
the  middle  of  the  DNA  fragment  (M)  in  comparison  to  the  mobility  of  complexes 
when  the  ERE  is  at  the  end  of  the  DNA  fragment  (16,  29).  The  calculated  bending 
angle  for  hERa  was  64.8  ±  1  (mean  ±  SEM,  n  =  3),  for  hERpl  was  53.6  ±  0.5  and 
for  hER32  was  54.6  +  0.7.  A  lower  overall  signal  of  the  retarded  complexes  were 
seen  with  hERpi  and  hER32  compared  to  hERawhich  likely  reflects  the  lower 
efficiency  of  the  hERpi  isoforms  compared  to  the  hERa  of  binding  to  an  ERE. 
Furthermore,  the  DNA  binding  ability  of  hER35  was  too  low  to  obtain  accurate  data 
for  DNA  bending  calculations.  No  effect  of  ligand  was  observed  (data  not  shown). 
Transcriptional  Activity  of  hERpi,  liER32  and  hERpS. 

The  ability  of  tagged  ERs  to  activate  transcription  was  initially  investigated  using 
Cos-1  cells  and  an  ERE  containing  reporter  gene,  ERE452-delta-TCO-CAT,  which  has 
two  vitellogenin  A2  (-333/-288)  EREs  upstream  of  a  CAT  reporter  (20).  Epitope 
tagged  receptors  were  used  so  that  relative  expression  of  all  the  ERs  could  be 
measured  using  antibodies  to  the  epitope  tag,  and  preliminary  experiments 
demonstrated  that  the  tagged  ERa  and  ERP 1  were  similar  to  their  untagged 
counterparts  in  activating  transcription  with  and  without  ligand  (data  not  shown). 
Preliminary  studies  showed  that  transfection  of  50  ng  ERa  expression  vector  gave 
maximal  estradiol  induced  transactivation  of  this  reporter  gene.  hERa  and  hERpi 
activated  transcription  in  a  ligand  inducible  manner  (Figure  6),  but  hERpi  was 


Peng  etal 


-  II  - 


overall  less  active  than  hERa(P  <  0.0001,  n  =  5),  and  increased  expression  of 
hERpi  did  not  alter  this  relationship.  These  data  are  consistent  with  previous 
findings.  As  shown  in  Figure  6,  low  doses  of  estradiol  (0.1  nM)  which  significantly 
activated  hERot  did  not  activate  hERpi  (P  <  0.0001,  n  =  5),  and  the  apparent  ligand 
independent  activity  (zero  ligand  added)  of  hERa  was  significantly  higher  than  that 
of  hERpi  (P  =  0.024,  n  =  5).  This  ligand-independent  activity  was  inhibited  by  0.1 
and  100  nM  of  the  antiestrogen  LYl  17018  (a  raloxifene  analog)  as  well  as  by  hER^l 
(50  and  450  ng)  and  hER32  (450  ng)(data  not  shown).  These  results  were  not  due  to 
different  levels  of  expression  of  hERaand  hER^l  since  Western  blot  analysis 
showed  similar  levels  of  expression  of  the  two  receptors  when  similar  amounts  of 
plasmid  were  transfected  (Figure  7).  The  expression  of  the  variant  isoforms  hERp2 
and  hERP5  alone  demonstrated  little  if  any  transcriptional  activity  under  these 
conditions  (Figure  6). 

ERa  and  hERP  isoforms  can  heterodimerize  (2),  which  may  underlie  the  functional 
interactions  between  ER  isoforms.  All  hERP  isoforms  tested  inhibited  the 
transcriptional  activity  of  hERa  on  an  ERE  containing  promoter  (Figure  8  A  and  8B) 
but  the  various  hERP  isoforms  had  different  efficiencies  with  hERpi  >  hERP2  > 
hERP5.  Ligand  activation  of  hERP  1  did  not  effect  its  ability  to  decrease  the  activity 
of  hERa,  since  under  conditions  when  it  was  not  activated  (0  or  0. 1  nM  estradiol,  see 
Figure  6)  hERpi  activity  was  similar  to  that  under  conditions  where  it  was  activated 
(100  nM  estradiol).  Variant  isoforms  of  hERP  had  little  if  any  effect  on  hERpi 
activity  on  ERE-containing  promoters  (data  not  shown). 

The  transcriptional  activity  of  ER  isoforms  was  next  examined  on  the  non-ERE- 
containing  promoter,  TGFP-3-CAT  where  the  DNA  binding  domain  of  ERa  is  not 
required  for  activity  (21).  This  promoter  was  shown  to  be  preferentially  activated  by 
the  raloxifene  bound  hERa  compared  to  estradiol  in  cultured  cells  (21),  and  we  have 
previously  shown  differential  abilities  of  murine  ERP  isoforms  to  effect  this  promoter 
compared  to  ERE-containing  promoters  (16).  Therefore,  the  activity  of  hERP 
isoforms  on  TGFP-3-CAT  was  examined  (Figure  9).  Optimal  activity  for  hERa  was 
obtained  with  transfection  of  50  ng  of  expression  plasmid  (data  not  shown).  A 
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significant  increase  in  transcription  was  obtained  with  0.1  nM  of  LY 11 70 18  that  was 
not  further  increased  with  lOOnM  LY  11 70 18  treatment  (P  =  0.0061,  n  =  3). 
LYl  17018  significantly  increased  the  transcriptional  activity  of  hERpi  on  the  TGPP- 
3-CAT  reporter  gene  at  the  lower  levels  of  hERpi  expression  (50ng,  P  =  0.008;  150 
ng,  P  =  0.02,  n  =  3)  but  at  high  levels  of  hERpi  expression  (450  ng),  a  significant 
increase  in  ligand-independent  activity  was  seen,  and  no  further  increase  was  seen 
due  to  ligand.  Overall,  hERpi  was  significantly  less  active  than  hERa  in  inducing 
TGFP-3-CAT  (P<0.0001,  n  =  3).  Although  there  was  a  trend  towards  inhibition  of 
TGFp-3-CAT  with  increasing  expression  of  hERP2  or  hERP5  (data  not  shown),  this 
was  not  statistically  significant. 

When  the  ability  of  hERP  isoforms  to  effect  hERa  activity  was  investigated  at  the 
TGF-P3  promoter,  differences  between  the  wild-type  and  variant  isoforms  were 
observed.  The  wild-type  hERpi  did  not  significantly  effect  hERa  transcriptional 
activity  at  any  level  of  expression  tested  (Figure  lOA  and  lOB).  However,  under  the 
same  conditions  hERP2significantly  inhibited  hERa  transcriptional  activity  on 
TGFP-3-CAT  (P=  0.0002,  n  =  3),  and  as  expected  the  effect  was  not  influenced  by 
LYl  17018,  since  hERp2does  not  bind  ligand.  However,  hERP2  inhibits  both  the 
ligand  activated  and  the  non-ligand  activated  (data  not  shown)  hERa  (P  =  0.017,  n  = 
3)  at  the  TGFP-3-CAT  promoter.  hERP5  also  inhibited  hERa  transcriptional  activity 
on  TGFP-3-CAT  but  only  at  the  highest  expression  of  hERP5  (Figure  lOB,  P  = 
0.038,  n  =  3).  Similar  to  their  action  at  an  ERE  containing  promoter,  the  truncated 
ERP  variants  2  and  5  do  not  modulate  wild-type  hERpi  transcriptional  activity  on 
TGFP-3-CAT  (data  not  shown). 

Discussion. 

There  is  a  growing  body  of  evidence  that  ERa  and  ERP  can  be  expressed  together  in 
some  cell  types  and  independently  expressed  in  others  (30-32).  If  expressed  together 
they  form  heterodimers,  which  under  experimental  conditions  are  preferred  over 
homodimerization  (2).  Further,  transient  coexpression  of  ERa  and  ERP  in  cell  lines, 
results  in  ERpi  induced  reduction  of  ERa  activity  at  low  ligand  concentrations,  as 
measured  using  ERE-regulated  reporters  (33).  A  conclusion  from  these  data  is  that 
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ERP  can  directly  modulate  ERa  activity.  This  has  significance  since  many  reports 
exist  of  differential  expression  of  the  two  receptors  under  conditions  of  altered 
estrogen  sensitivity.  For  example,  ERP  expression  is  significantly  downregulated  and 
ERa  expression  upregulated  during  human  breast  tumorigenesis,  suggesting  that 
ERP’s  abilityto  modulate  ERa  is  significantly  altered  during  breast  tumorigenesis 
(12,  13).  In  addition,  current  data  show  that  in  normal  and  neoplastic  breast  tissues, 
the  level  of  expression  of  the  C-terminally  truncated  ER|3  variants,  ERP2and 
ER35,  is  markedly  higher  than  the  ligand  binding  ER3I .  These  data  suggest  that 
the  variant  ER3  isoforms  may  also  have  a  role  in  modulating  estrogen  and  possibly 
antiestrogen  action  in  human  breast  cells.  The  experiments  described  in  this 
manuscript  were  undertaken  to  gain  insight  into  the  possible  role  of  the  truncated 
ER3  variants. 

Our  data  show  that  only  hER3l  is  able  to  bind  ligand.  Steroid  hormone  receptors 
are  known  to  undergo  conformational  changes  during  the  process  of  activation 
especially  due  to  ligand  binding  and  differences  are  seen  between  agonist  and 
antagonist  binding  (17,  28).  Recent  structural  analyses  of  the  Ligand  Binding  Domain 
(LBD)  of  several  nuclear  receptors  suggest  that  the  LBD  contains  common  structural 
motifs  that  generate  a  conserved  ligand  binding  pocket,  and  that  agonist  and 
antagonists  bind  to  the  same  site  but  induce  different  conformational  changes  that  are 
now  known  to  affect  transcriptional  function,  providing  structural  evidence  for 
antagonism  (34).  The  variant  hER3  isoforms,  while  not  binding  ligand,  may  exist  in 
an  activated  state  in  the  absence  and  presence  of  ligand,  however,  our  data  suggest 
that  hER32  and  hER35  are  unlikely  to  be  in  an  activated  conformation,  and  this  is 
consistent  with  their  inability  to  activate  transcription  of  either  a  ‘classical’  or  a’ non- 
classical’  estrogen  receptor  regulated  reporter  gene. 

All  ER3  isoforms  examined  (ER31, 2,  and  5)  inhibit  the  transcriptional  activity  of 
ERa  on  ERE  containing  promoters,  while  only  ER31  has  any  activity  alone.  This 
confirms  and  extends  previous  data  and  demonstrates  the  relative  inhibitory  activity 
of  the  ER3  isoforms  is  ER3I  ^R32  >ER35.  This  correlates  with  the  relative 
efficiencies  with  which  ER3  homodimers  bind  to  DNA  and  may  suggest  a 
competition  of  the  beta  isoform  homodimers  with  ERa  homodimers  for  DNA 
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binding.  However,  since  heterodimers  are  preferred  under  these  conditions,  it  is  likely 
that  these  predominant  under  our  experimental  conditions  and  the  intrinsically  lower 
transcriptional  activity  of  the  heterodimers  predominants.  Cowley  et  al,  (2) 
demonstrated  that  when  hERa  and  hERpi  are  expressed  at  both  a  1:1  and  1:2  ratio 
the  ERot/ERpi  heterodimer  was  predominant.  This  heterodimer  had  a  DNA  binding 
affinity  similar  to  that  of  the  ERa  homodimer,  and  was  capable  of  recruiting  SRC-1. 
However,  the  heterodimer  was  less  transcriptional  activity  than  the 
ERoc/ERa  homodimer,  suggesting  that  it  may  be  less  efficient  in  recruiting 
coactivators  than  the  ERahomodimer.  In  contrast,  the  C-terminally  truncated 
hERP2  has  markedly  reduced  ability  to  bind  to  DNA  and  likely  the  ERa/ERp2 
heterodimer  also  binds  less  well  than  ERcc/ERa  homodimers  to  an  ERE  (6,  10).  But 
in  contrast  to  hERpi,  hERp2  does  not  recruit  coactivators  (10).  Our  data  show  that 
hERpS  is  less  efficient  than  hERP2  in  binding  to  DNA,  and  is  also  unlikely  to  recruit 
coactivators.  However,  at  an  ERE  the  wild-type  hERpi  is  more  potent  than  either  of 
the  two  variants  in  inhibiting  the  ability  of  ERa  to  activate  transcription.  So  it  appears 
that  the  inability  to  recruit  coactivators  is  not  correlated  with  the  ability  of 
ERP  isoforms  to  inhibit  ERa  activity.  Since  DNA  binding  activity  is  also  a  reflection 
of  efficiency  of  dimerization,  it  is  speculated  that  the  truncated  ERP  isoforms  have 
reduced  ability  to  dimerize  with  ERa  and  form  stable  heterodimers  than  the  wild-type 
ERP  1.  Together  with  our  Western  blot  data  it  seems  that  significant  inhibition  of 
ERa  transcription  occurs  at  levels  of  ERpi  expression  that  are  less  than  or  equivalent 
to  ERa.  Therefore  our  data  would  be  consistent  with  the  mechanism  of  inhibition 
being  related  to  a  high  efficiency  of  dimerization  and  reduced  efficiency  in  recruiting 
coactivators,  but  not  the  inability  to  recruit  coactivators. 

Interestingly,  marked  differences  in  the  ability  of  the  ERP  isoforms  to  affect 
ERa  activity  are  seen  at  an  estrogen  receptor  responsive  site  where  the  mechanism  of 
transcriptional  regulation  is  quite  distinct  from  that  operating  at  a  classical  ERE,  e.g. 
the  so-called  raloxifene  responsive  element  in  the  TGF-PS  promoter  (21).  This  is  in 
marked  contrast  to  the  results  seen  at  an  ERE  regulated  reporter  gene.  The  ER 
responsive  site  in  the  TGF-p3  promoter  is  poorly  activated  by  the  estradiol-ERa 
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complex,  but  strongly  activated  by  the  raloxifene-ERa  complex.  In  addition  the  DNA 
binding  domain  of  the  ER  is  not  required  for  this  activation.  It  is  assumed  that 
proteinrprotein  interactions  between  ERa  and  other  transcription  factors  bound  to  this 
promoter  are  involved  in  regulation.  However,  the  identity  of  these  ‘other’ 
transcription  factors  is  unknown.  Using  an  analogue  of  raloxifene,  LYl  17018  (16), 
we  have  confirmed  that  this  promoter  is  poorly  activated  by  the  estradiol-ERa 
complex  (and  this  was  not  altered  in  our  hands  by  treatment  of  the  transfected  Cos-1 
cells  with  EOF  (35),  data  not  shown)  but  was  significantly  activated  by  the 
LYl  17018-ERcc  complex.  Similarly,  the  LYl  17018-  hERpi  complex  was  found  to 
activate  transcription  from  the  TGF-P3  promoter,  but  in  contrast  to  the  murine  ERpi 
(16),  is  less  active  than  the  LY117018-hERa  complex.  Human  ERP2  and  hERp5 
alone  could  not  activate  this  promoter.  This  is  in  contrast  to  the  murine  ERP2  variant 
(16)  which  is  structurally  quite  different  to  the  hER|32  Furthermore,  no  murine 
equivalent  to  either  hERP2orhERP5  isoforms  that  are  frequently  expressed  in 
human  tissues,  have  as  yet  been  identified.  However,  coexpression  of  increasing 
amounts  of  hERP2  and  hERpS  with  ERa  resulted  in  inhibition  of  LYl  17018-ERa 
transcriptional  activity  but  not  LY117018-ERP1  activity  from  the  TGF-P3  promoter. 
In  contrast  to  an  ERE  containing  promoter  is  the  observation  that  the  wild-type 
hERpi  did  not  significantly  inhibit  the  transcriptional  activity  of  the  LYl  17018-ERa 
complex  at  the  TGF-p3  promoter.  At  this  promoter  the  differences  in  the  hERp 
isoform  activity  on  LYl  17018-ERa  complexes,  was  correlated  to  the  ability  to  recruit 
coregulatory  factors.  Significant  effects  of  hERP2  on  hERa  were  seen  under 
conditions  of  equimolar  expression,  as  determined  by  Western  blot  analysis  of  the 
similarly  tagged  proteins,  but  hERpSwas  less  active  than  hERp2  and  this  is 
consistent  with  a  reduced  efficiency  of  dimerization.  There  appears  to  be  differential 
expression  of  hERP  isoforms  at  least  at  the  RNA  level  in  different  human  tissues,  as 
well  as  altered  relative  expression  during  breast  tumorigenesis  (9,  36),  and  altered 
levels  of  hERPcx  (hERp2)  as  well  as  hERpi  during  prostate  cancer  progression  (8). 
Therefore,  it  is  possible  that  the  differential  activities  of  hERP  isoforms  on  some 
genes  may  have  both  physiologically  and  pathophysiologically  importance. 
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In  conclusion  we  have  characterized  some  potential  functions  of  several 
commonly  expressed  hERP  isoforms.  Generally,  the  ligand  binding  wild-type  hER^l 
has  transcriptional  activity  alone  on  both  ‘classical’  and  ‘non-classical’  estrogen 
responsive  promoters,  although  it  is  less  efficient  than  ERa.  Furthermore,  the 
hERP  family  of  receptors  generally  negatively  modulate  ERa  transcriptional  activity 
when  coexpressed,  at  ‘classical’  as  well  as  ‘non-classical’  ER  responsive  promoters. 
However,  promoter  specific  differential  activity  of  the  various  hERP  isoforms  was 
found,  in  particular  between  the  wild-type  hERpi  and  its  C-terminally  truncated 
variants  hERp2  and  hERps.  The  possibility  that  there  is  differential  expression  of  the 
hERP  isoforms  suggests  that  they  may  have  a  role  in  differentially  modulating 
estrogen  action. 
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Legends  to  Figures. 

Figure  1. 

A.  cDNA  structure  of  human  ERpl,  human  ERP2  and  human  ERPS. 

B.  Genomic  structure  of  the  human  estrogen  receptor  P  locus  on  chromosome  14. 
Human  ERP  cDNA,  ESTs  and  published  promoter  sequences  (accession  numbers 
in  the  text)  where  aligned  with  sequences  from  two  genomic  clones  of  human 
chromosome  14  (AL162756/CNS01RHJ  and  AF215937).  Exons  1C  to  IH  are 
found  only  in  one  cDNA,  accession  number  AB006589). 


Figure  2. 

Predicted  proteins  of  the  human  ERP  isoform  cDNAs. 

Figure  3. 

Determination  of  estradiol  binding  to  hERpi  (A),  hERp2(B),  and  hERp5(C). 
Increasing  amounts  of  ^H-estradiol-17p  (0.01-100  nM)  were  incubated  with  a 
constant  amount  of  each  in  vitro  transcribed/translated  hERP  protein.  Two  separate 
experiments  were  performed.  Total  binding  is  shown  by  the  closed  circles  and  non¬ 
specific  binding  is  shown  by  the  empty  circles. 

Figure  4. 


Sensitivity  of  human  Estrogen  Receptor  Isoforms  to  Protease  Digestion. 
Radiolabelled  ER  was  made  in  vitro  as  described  in  the  Material  and  Methods 
section,  and  digested  with  increasing  levels  of  trypsin,  with  and  without  ligand.  The 
products  were  visualized  by  autoradiography  after  SDS-PAGE.  Resistant  bands  are 
shown  by  arrows.  The  asterisk  shows  agonist  induced  resistant  32.5  kDa  band. 

Figure  5. 

Determination  of  the  ability  of  hERa,  hERpi,  hERp2  and  hERps  to  bind  to  DNA. 
Autoradiograph  of  an  electrophoretic  mobility  gel  shift  analysis  of  in  vitro 
transcribed/translated  hERa,  hERpi,  hERp2  and  hERps  proteins  binding  to  a  35  mer 
double  stranded  ERE  oligonucleotide  containing  the  vitellogenin  A2  ERE  sequence. 
Free  ERE  and  the  shifted  complexes  are  indicated.  The  presence  of  the  appropriate 
ER  isoform  in  the  shifted  complex  was  determined  by  the  ability  of  a  specific 
antibody  (H222  for  ERa;  PAl  for  the  ERP  proteins)  to  super-shift  the  complex,  data 
not  shown.  Specificity  of  the  complexes  was  determined  by  the  ability  of  a  200  fold 
excess  of  the  unlabelled  ERE  (specific  competitor)  to  compete  for  the  shifted 
complex  and  non-specific  interactions  were  determined  using  a  200  fold  excess  of 
unlabelled  nonspecific  33  mer  oligonucleotide  (nonspecific  competitor). 
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Figure  6. 

Transcriptional  activity  of  hERa,  hERpl,  hERP2  and  hERps  on  ERE  regulated  CAT 
reporter  gene.  Increasing  amounts  of  HER  isoform  expression  vectors  were 
cotransfected  with  an  ERE-CAT  reporter  gene  into  Cos-1  cells.  The  transfected  cells 
were  treated  for  48  hours  with  or  without  estradiol- 17p.  Results  show  the  mean  CAT 
activity  after  correction  for  P-galactosidase  activity  (transfection  efficiency)  +  SEM 
of  5  independent  experiments  for  hERa,  hERpi  and  hERP2,  and  3  independent 
experiments  for  hERps.  See  text  for  statistical  analysis. 


Figure  7. 

Western  blot  analysis  of  Cos-1  cell  extracts  48  hours  after  transfection  of  the 
indicated  amounts  of  tagged  -  hERa,  -  hERpi  and  -  hERP2  expression  vectors,  as 
detailed  in  Materials  and  Methods.  Lane  1  is  extract  of  Cos-1  cells  transfected  with 
empty  expression  vector  alone;  lane  2  is  tagged  -  hERa  expression  vector;  lane  3  is 
tagged  -  hERa  (50ng)  +  tagged  -  hERpl  (50  ng)  expression  vectors;  lane  4  is  tagged 
-  hERa  (50  ng)  +  tagged  -  hERpl  (150  ng)  expression  vector;  lane  5  is  tagged  - 
hERa  (50ng)  +  tagged  -  hERP2  (50  ng)  expression  vectors;  lane  6  is  tagged  -  hERa 
(50ng)  +  tagged  -  hERp2  (450  ng)  expression  vectors.  The  tagged  proteins  were 
visualized  with  anti-Xpress  antibody  as  described  in  Materials  and  Methods. 


Figure  8. 

Effect  of  increasing  amounts  of  coexpressed  hERP  isoforms  on  the  ability  of 
hERa  (50  ng)  to  activate  transcription  from  an  ERE  (vitellogenin  A2)  regulated  CAT 
reporter  gene  in  the  presence  and  absence  of  ligand  following  co-transfection  into 
Cos-1  cells.  The  results  show  the  mean  CAT  activity  after  correction  for 
P-galactosidase  activity  (transfection  efficiency)  ±  SEM  of  3  independent 
experiments.  A.  Effect  of  hERpl  and  hERp2  on  hERa;B.  Effect  of  hERp5  on 
hERa.  See  text  for  statistical  analysis. 


Figure  9. 

Transcriptional  activity  of  hERa,  hERpl  and  hERP2  on  a  novel  raloxifene  responsive 
element  regulated  reporter  gene,  TGF-P3-CAT.  Increasing  amounts  of  hER  isoform 
expression  vectors  were  cotransfected  with  a  TGF-P3-CAT  reporter  gene  into  Cos-1 
cells.  The  transfected  cells  were  treated  for  48  hours  with  or  without  or  the  raloxifene 
analog  LYl  17018.  The  results  show  the  mean  CAT  activity  after  correction  for 
P-galactosidase  activity  (transfection  efficiency)  +  SEM  of  3  independent 
experiments.  See  text  for  statistical  analysis. 


Figure  10. 


Effect  of  increasing  amounts  of  coexpressed  hERP  isoforms  on  the  ability  of 
hERa  (50  ng)  to  activate  transcription  from  a  TGF^3-CAT  reporter  gene  in  the 
presence  and  absence  of  the  raloxifene  analog  LY 11 70 18  following  co-transfection 
into  Cos-1  ceils.  The  results  show  the  mean  CAT  activity  after  correction  for 
p-galactosidase  activity  (transfection  efficiency)  ±  SEM  of  3  independent 
experiments.  A.  Effect  of  hERpi  and  hERP2  on  hERa;  B.  Effect  of  hERp5  on 
hERa.  See  text  for  statistical  analysis. 


characterized  promoter  (AF 191 544) 


Figure  2 


-Estradiol-17p  Binding  (fmol)  H  -Estradiol-17p  Binding  (pmol)  H^-Estradiol-17p  Binding  (fmol) 


hERpi 


Figure  3 


Estrogen  Receptor  Alpha 

Ugand"(100  nM)  -  E2 

Trypsin  (ug)  q  o2  10  50  0  02  10  50 

4^1  kDa  ^  -*■  ^ 

16.5  kDa ->  -M 


Estrogen  Receptor  Alpha 

_ 4^0H-TAM  LY 

0  0.2  1.0  5.0  0  0.2  1.0  5.0 


I  Resistant 
jMfci  <  Bands 


Estrogen  Receptor  Beta1 

Ligand  (100  nM)  --  E2 


47.5  kDa  ^ 

32.5  kDa 


25.0  kDa 


16.5  kDa 


ww 


Estrogen  Receptor  Betal 
4-OH-TAM  Ly" 


0  0.2  1.0  5.0  Q  0.2  1.0  5.0 


^  Resistant 
Bands 


Estrogen  Receptor  Beta2 

Ligand  (100  nM)  -  E2 

Trypsin  (ug)  q  0.2  1.0  .5  0  n  nz  in  «;n 
kga  X 

62  kDa  ^ 


47.5  kDa  -> 

32.5  kDa  -> 
25.0  kDa  -> 


16.5  kDa 


V’-' 


Estrogen  Receptor  Beta2 
4-OH-TAM  LY 


_0  0.2  1.0  5.0  0  0.2  1.0  5.0 


^  Resistant 
<  Bands 


Estrogen  Receptor  Beta5 

Ligand  (100  nM)  -  E2 

5.9  D0|  1.0  S.0 


47.5  kDa  g 

32.5  kDa  ^ 
25.0  kDa  ->  ^ 


'igM 

gg 


16.5  kDa 


Estrogen  Receptor  Betas 
4-OH-TAM  LY 

_0_  0.2  1.0  5.0  0  0.2  1.0  fi  n 

liB  Jli  Wl  I 


■'fSJ  |g|E 


^  Resistant 
Bands 


Figure  7. 


1  2  3  4  5  6 


Tag-ERa 

Tag-ERp 


Figure  8A. 


O 

O 


o 

o 


+ 

1 

o 

V) 

Tf 

+ 

1 

o 

+ 

450 

1 

+ 

o 

tn 

1 

+ 

50 

1 

+ 

1 

1 

I  I 


I  I 

I  I 


o 

o 

o 

ro 


ieO-d/^»!A!»3VlV3 


1 

+ 

1 

o 

If) 

1 

+ 

1 

O 

o 

If) 

1 

+ 

If) 

Tf 

1 

1 

+ 

150 

1 

1 

+ 

50 

1 

1 

+ 

1 

1 

B  ©f»  WD 
O  0  £2 

»/)  W  sS 

q  C2D.  oa 

^ 

M  ^  'tq 


200 


CAT  Activity/p-Gal 


QC 

O 

o 


1000 


Figure  lOA 


Figure  lOB 


ERa  (50ng) 
ERP5  (ng) 


APPENDIX  4 
in  press,  Brit.  J.  Cancer 


APPENDIX  4 


Relationship  of  co-regulator  and  estrogen  receptor  isoform  expression  to  de  novo  tamoxifen 

resistance  in  human  breast  cancer.^ 

Leigh  C  Murphy^,  Etienne  Leygue,  Yulian  Niu,  Linda  Snell,  S-M  Ho 

md  Peter  H  Watson. 

Affiliations  of  authors;  Manitoba  Institute  of  Cell  Biology  (LCM),  Department  of  Biochemistry  and 
Medical  Genetics  (E.L.,  Y.N,,  L.C.M.)  and  Department  of  Pathology  (L.S.,  P.H.W.),  University  of  Manitoba, 

Faculty  of  Medicine,  Winnipeg,  Manitoba,  Canada,  R3EOW3;  SMH . 

Running  title:  ER  isoform  and  coregulator  expression  in  primary  breast  tumors 

Key  words;  SRA,  AIB 1 ,  REA,  co-regulators,  ER  isoforms,  human  breast  cancer,  tamoxifen. 


Footnotes: 

1  This  work  was  supported  by  grants  from  the  Canadian  Breast  Cancer  Research  Initiative  (CBCRI),  the 
Canadian  Institutes  for  Health  Research  (CIHR)  and  the  U.S.  Army  Medical  Research  and  Materiel 
Command  (USAMRMC).  P.  H.  W.  is  a  recipient  of  a  CIHR  Scientist  Award  and  a  US  Army  Acadanic 
Award.  E.L.  has  a  USAMRMC  Career  Development  Award. 

2  To  whom  requests  for  reprints  should  be  addressed,  at  Manitoba  Institute  of  Cell  Biology,  CancerCare 
Manitoba,  675  McDermot  Ave,  Winnipeg,  MB.  R3E  OV9,  Canada.  Phone:  (204)  787-4071;  Fax;  (204)  787 
2 1 90;  E-mail :  lcmurph@cc.umanitoba.ca 

2  The  abbreviations  used  are:  ER,  estrogen  receptor;  PR,  progesterone  receptor;  SRA,  steroid  receptor  RNA 
activator;  AIB  1 ,  amplified  in  breast  cancer  1 ;  REA,  repressor  of  estrogen  receptor  activity. 


Abstract 


This  study  addresses  the  hypothesis  that  altered  expression  of  estrogen  receptor-beta  and/or  altered  relative 
expression  of  coactivators  and  corepressors  of  estrogen  receptors  are  associated  with  and  may  be  mechanisms 
oide  novo  tamoxifen  resistance  in  estrogen  recq)tor  positive  breast  tumors.  All  cases  were  ER+,  node 
negative,  primary  breast  tumors  from  patients  who  later  had  no  disease  progression  (tamoxifen  sensitive)  or 
whose  disease  progressed  while  on  tamoxifen  (tamoxifen  resistant).  Using  an  antibody  to  estrogen  receptor- 
beta  that  detects  multiple  forms  of  this  protein  (total)  but  not  an  antibody  that  detects  only  wild-type  estrogen 
receptor-beta  1,  it  was  found  that  high  total  ER  beta  protein  expressors  were  more  frequently  observed  in 
tamoxifen  sensitive  tumors  than  resistant  tumors  (Fishers  exact  test,  P=0.046).  However,  no  significant 

differences  in  the  relative  expression  of  ERp2,  ERp5  and  wild-type  ERpl  RNA  in  the  tamoxifen  saisitive  and 

resistant  groups  were  found.  As  well,  when  the  relative  expression  of  two  known  coactivators  SRA  and  ABl 
RNA  to  the  known  corepressor,  REA  RNA,  was  examined,  no  significant  differences  between  the  tamoxifen 

sensitive  and  resistant  groups  was  found.  Altogether,  there  is  little  evidence  for  altered  coregulators 
expression  in  breast  tumors  that  are  de  novo  tamoxifen  resistant.  However,  our  data  provide  preliminary 
evidence  that  the  expression  of  ERp  protein  isoforms  may  differ  in  primary  tumors  of  breast  cancer  patients 
who  prove  to  have  differential  sensitivity  to  tamoxifen  therapy. 

Introduction 

The  ability  of  anti-estrogens  such  as  tamoxifen  to  compete  with  estrogens  for  binding  to  ER  and  to  antagonize 
their  mitogenic  action  provides  the  basic  rationale  for  endocrine  therapy  and  prevention  (for  a  review  see  (1)) 
in  breast  cancer.  Adjuvant  tamoxifen  post-operative  therapy  reduces  the  nianber  of  recurrences  and  prolongs 
survival  in  women  whose  primary  tumors  are  ER  positive  (2).  However,  even  though  ER  level  is  considered  a 
marker  for  predicting  the  likelihood  of  responding  to  adjuvant  hormonal  therapies,  some  patients,  whose 
primary  tumors  are  ER  positive  do  not  respond  to  tamoxifen  treatment.  Such  apparent  de  novo  tamoxifen 
resistance  does  not  depend  upon  the  level  of  ER  within  the  primary  tumor.  As  well  many  of  those  patients 
whose  disease  initially  responds  to  tamoxifen,  progress  under  tamoxifen  having  acquired  resistance  to 
tamoxifen  and  this  occurs  despite  continued  expression  of  ER.  Thus  suggesting  other  components  of  the 
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estrogen  signaling  pathway  may  be  altered.  Recent  observations  using  laboratory  models  (3-6)  have 
demonstrated  that  altered  levels  of  ER  isoforms  and/or  alteration  of  expression  of  co-activators  and  co¬ 
repressors  can  result  in  altered  estrogen  and  antiestrogen  activity  in  target  cells,  suggesting  the  hypothesis  that 
altered  levels  of  ER  isoforms  and/or  co-regulators  in  vivo  could  be  a  mechanism  of  tamoxifen  resistance. 
Previously  we  have  demonstrated  that  the  relative  expression  of  ERo/ERp  as  well  as  the  relative  expression  of 
some  ER  coactivators  to  corepressors  is  sigmficantly  altered  during  breast  tumorigenesis  in  vivo  (7,  8). 
Furthermore,  since  these  alterations  parallel  the  marked  changes  in  estrogen  action  that  accompany  breast 
tumorigenesis,  they  may  have  a  role  in  this  process.  To  explore  the  hypothesis  that  such  changes  could 
underlie  de  novo  tamoxifen  resistance  in  vivo  the  expression  of  ER  isoforms,  two  known  co-activators 
(steroid  receptor  RNA  activator,  SRA,  (4)  and  amplified  in  breast  cancer-1,  AIBl  (9))  and  one  co-repressor 
(repressor  of  estrogen  receptor  activity,  REA  (5))  of  ER  activity  have  been  investigated  in  primary  breast 
tumors  fi'om  node  negative  patients  whose  tumors  were  ER  positive  and  that  subsequently  responded  or  bad 
disease  progression  while  on  adjuvant  tamoxifen  therapy. 

Materials  and  methods 
Human  breast  tumors. 

All  breast  tumor  cases  used  for  this  study  were  selected  fi-om  the  NCIC-Manitoba  Breast  Tumor  Bank 
(Winmpeg,  Manitoba,  Canada).  As  previously  described  (10),  tissues  are  accrued  to  the  Bank  fi'om  cases  at 
multiple  centers  within  Manitoba,  rapidly  collected  and  processed  to  create  matched  formalin-fixed-embedded 
and  fi-ozen  tissue  blocks  for  each  case  with  the  mirror  image  surfaces  oriented  by  colored  inks.  The  histology 
of  every  sample  in  the  Bank  is  uniformly  interpreted  by  a  pathologist  in  Hematoxylin/Eosin  (H&E)  stained 
sections  fi'om  the  face  of  the  paraffin  tissue  block.  This  information  is  available  in  a  computerized  database 
along  with  relevant  pathological  and  clinical  information  and  was  used  as  a  guide  for  selection  of  specific 
paraffin  and  firozen  blocks  from  cases  for  this  study.  For  each  case,  interpretations  included  an  estimate  of  the 
cellular  composition  (including  the  percentage  of  invasive  epithelial  tumor  cells  and  stroma),  tumor  type  and 
tumor  grade  (Nottingham  score).  Steroid  receptor  status  was  determined  for  all  cases  by  ligand  binding  assay 
performed  on  an  adjacent  portion  of  tumor  tissue.  Tumors  with  estrogen  receptor  levels  above  3  finol/mg  of 


total  protein  were  considered  ER  positive. 

To  identify  cases  diat  responded  divergently  to  tamoxifen,  review  of  approximately  1000  consecutive 
cases  was  undertaken  to  identify  cases  that  were  ER  positive,  node  negative  and  that  had  been  treated  with 
adjuvant  tamoxifen  following  surgery  +/-  local  radiation.  From  these  the  first  cohort  of  27  cases  was  selected 
to  include  a  subset  (n=13)  that  had  shown  progression  of  disease  (either  died  or  alive  with  recurrent  disease, 
referred  to  as  tamoxifen  resistant  cases)  and  a  similar  control  subset  (n=14)  specifically  selected  to  comprise 
cases  with  similar  lengths  of  follow-up  (resistant  34  versus  sensitive  39  months),  ER  status,  tumor  grade  and 
tumor  histology,  but  that  had  shown  no  progression  of  disease  (referred  to  as  tamoxifen  sensitive  cases).  The 
tumor  characteristics  were:  1)  ‘Tam  Sensitive’  group  median  ER  was  60.5  finol/mg  protein  (range  6  —146 
finol/mg  protein),  median  PR  was  32  finol/mg  protein  (range  8-216  finol/mg  protein);  median  grade  was  5 
(range  4-8);  median  age  at  biopsy  was  69  years  (range  35-87  years);  median  follow-up  time  was  56  months 
(range  18-  79);  2)  Tam  Resistant’  group  median  ER  was  57  finol/mg  protein  (range  4-136  finol/mg 
protein);  median  PR  was  14  finol/mg  protein  (range  4  —  288  finol/mg  protein);  median  grade  was  6  (range  4- 
9);  median  age  at  biopsy  was  67  years  (range  49-83  years);  median  follow-up  time  was  56  months  (range  9  - 
85). 

For  the  RNA  studies,  fi'ozen  tissue  corresponding  to  the  blocks  for  several  of  the  first  cohort  of  older 
cases  used  above,  were  not  available.  Therefore  a  second  study  cohort  was  selected  that  had  fi-ozen  tissue 
available.  The  relevant  patient^tumor  characteristics  were  similar  to  the  above  cohort,  although  the  follow-up 
time  was  shorter:  1)  ‘Tam  Sensitive’  group  (n=16)  median  ER  was  37.5  finol/mg  protein  (range  4.4  -146 
finol/mg  protein),  median  PR  was  44  finol/mg  protein  (range  13.1  —216  finol/mg  protein);  median  grade  was 
6  (range  4-9);  median  age  at  biopsy  was  72  years  (range  47-87  years);  median  follow-up  time  was  39  months 
(range  13-  76);  2)  Tam  Resistant’  group  (n=16)  median  ER  was  21.5  finol/mg  protein  (range  5.6  -  107 
finol/mg  protein);  median  PR  was  14.3  finol/mg  protein  (range  7.8  -  288  finol/mg  protein);  median  grade  was 
6  (range  4-9);  median  age  at  biopsy  was  71  years  (range  60-89  years);  median  follow-up  time  was  34  months 
(range  9  -  63). 

Immunohistochemistry.  Inununohistochemistry  was  performed  on  serial  5  pm  sections  from  a 
representative,  formalin  fixed  paraffin  embedded  archival  tissue  block  from  each  tumor.  Immunohistochemical 
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staining  for  ERp  was  performed  using  two  different  primary  antibodies.  IgYERB503  (a  gift  from  Dr  Jan-Ake 
Gustafson)  detects  total  ERp  isoforms  (11, 12)  and  GC17  (a  gift  from  Dr  Shuk-Mei  Ho)  detects  only  the  wild- 
type  ERp  (13).  Antibodies  were  applied  using  an  automated  tissue  immunostainer  (Discovery  module,  Ventana 
Medical  Systems,  Phoenix,  AZ,  U.S.A.),  DAB  immunohistochemistry  kit  and  bulk  reagents  that  were  supplied 
by  the  manufacturer.  Briefly,  the  Discovery  staining  protocol  was  set  to  “Standard  Cell  Conditioning” 
procedure,  followed  by  12  hours  incubation  with  primary  antibody  and  32  minutes  incubation  with  secondary 
antibody.  Concentrations  of  primary  antibodies  initially  applied  to  the  Ventana  instrument  were  1 :200  for 

IgYERB503  and  1:50  for  GC17,  which  translates  into  final  concentrations  of  1:600  and  1:150  after  a  1:3 
dilution  with  buffer  dispensed  onto  the  slide  with  the  primary  antibody.  Levels  of  nuclear  ERp  expression  were 
scored  semi-quantitatively,  under  the  light  microscope.  Scores  were  obtained  by  estimating  average  signal 
intensity  (on  a  scale  of  0  to  300)  and  the  proportion  of  epithelial  cells  showing  a  positive  signal  (0,  none;  0.1, 
less  than  one  tenth;  0.5,  less  than  one  half;  1.0  greater  than  one  halQ.  The  intensity  and  proportion  scores  were 
then  multiplied  to  give  an  overall  H-score.  Cases  with  a  score  lower  than  or  equal  to  100  were  considered 
negative  or  weakly  positive,  whereas  tumors  with  scores  higher  than  100  were  classified  as  positive  for  ERp 
expression  (14). 

RNA  Extraction  and  RT-PCR  conditions.  Total  RNA  was  extracted  from  20  pm  frozen  tissue  sections 
(20  sections  per  tumor)  using  Tiizol™  reagent  (Life  Technologies,  NY)  according  to  the  manufacturer’s 
instructions  and  quantified  spectrophotometrically.  One  pg  of  total  RNA  was  reverse  transcribed  in  a  final 
volume  of  25  pi  as  previously  described  (15). 

Primers  and  PCR  conditions. 

Coregulators:  The  primers  used  were:  SRAcoreU  primer  (5’-  AGGAACGCGGCTGGAACGA  -3’;  sense; 
positions  35-53,  Genbank  accession  number  AF092038)  and  SRAcoreL  primer  (5’- 
AGTCTGGGGAACCGAGGAT  -3’;  antisense;  position  696-678,  Genbank  accession  number  AF092038); 
AIBl-U  primer  (5’- ATACTTGCTGGATGGTGGACT  -  3’;  sense;  positions  110-130,  Genbank  accession 
number  AF012108)  and  AIBl-L  primer  (5’- TCCTTGCTCTTTTATTTGACG  -3’;  antisense;  positions  458 
-  438,  Genbank  accession  number  AFO 1 2 1 08); 
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REA-U  primer  (5’-  CGAAAAATCTCCTCCCCTACA  -  3’;  sense;  positions  385-405,  Genbank  accession 
number  AF150962)  and  REA-L  primer  (5’-  CCTGCTTTGCTTTTTCTACCA  -3’;  antisense;  positions  781- 
76 1 ,  Genbank  accession  number  AF 1 50962). 

Radioactive  PCR  amplifications  for  SRA  were  performed  and  PCR  products  anal3^ed  as  previously  described 
(16)  with  minor  modifications.  Briefly,  1  jil  of  reverse  transcription  mixture  was  amplified  in  a  final  volume  of 

1 5  fil,  in  the  presence  of  1 .5  pCi  of  (a-^^p)  dCTP  (3000  Ci/mmol),  4  ng/pl  of  each  primer  and  0.3  unit  of  Taq 
DNA  polymerase  (Gibco  BRL,  Grand  Island,  NY).  For  SRA  each  PCR  consisted  of  30  cycles  (30  sec  at 
94°C,  30  sec  at  60°C  and  30  sec  at  72°C).  PCR  products  were  then  separated  on  6%  polyacrylamide  gels 
containing  7M  urea.  Following  electrophoresis,  the  gels  were  dried  and  exposed  2  hours  to  a  Molecular 
Imager™-FX  Imaging  screen  (Bio-Rad,  Hercules,  CA). 

PCR  amplifications  for  AIBl  and  REA  were  performed  and  PCR  products  analyzed  as  previously  described 
(15)  witii  minor  modifications.  Briefly,  1  pi  of  reverse  transcription  mixture  was  amplified  in  a  final  volume  of 
20  pi,  in  the  presence  of  4  ng/pl  of  each  primer  and  0.3  unit  of  Taq  DNA  polymerase  (Gibco  BRL,  Grand 
Island,  NY).  For  AIBl,  each  PCR  consisted  of  30  cycles  (30  sec  at  94°C,  30  sec  at  55°C  and  30  sec  at  72°C). 
For  REA  each  PCR  consisted  of  30  cycles  (30  sec  at  94°C,  30  sec  at  57°C  and  30  sec  at  72°C).  PCR  products 
were  then  separated  on  agarose  gels  stained  with  ethidium  bromide  as  previously  described  (15). 

Primers  for  ER  isoforms: 

ERa-U  primer  (5'-TGTGCAATGACTATGCTTCA-3';  sense;  located  in  ERa  792-811)  and  ERa-L  primer 
(5'-GCTCTTCCTCCTGTmTA-3';  antisense;  located  in  ERa  940-922).  Nucleotide  positions  given 
correspond  to  published  sequences  of  the  human  ERa  cDNA  (17).  PCR  amplifications  were  performed  and 
PCR  products  analyzed  as  previously  described  with  minor  modifications  (18).  Briefly,  1  pi  of  reverse 

transcription  mixture  was  amplified  in  a  final  volume  of  15  pi,  in  the  presence  of  1  pCi  (a-^^p)  dCTP  (3000 
Ci/mmol),  2  ng/pl  of  ERa-U/ERa-L  and  0.3  unit  of  Taq  DNA  polymerase  (Gibco  BRL,  Grand  Island,  NY). 
Each  PCR  consisted  of  30  cycles  (30  sec  at  94°C,  30  sec  at  60°C  and  30  sec  at  72°C). 

A  previously  validated  triple  primer  assay  was  used  to  determine  the  relative  expression  of  ERpi  and  its 
variant  isoforms  ERp2  and  ERp5  (19).  Briefly,  1  pi  of  reverse  transcription  mixture  was  amplified  in  a  final 


volume  of  15  ^il,  in  the  presence  of  1  jnCi  of  (a-32p)  dCTP  (3000  Ci/mmol),  4  ng/pl  of  each  primer  (ERplU, 
ERpiL  and  ERp2L)  and  0.3  unit  of  Taq  DNA  polymerase  (Gibco  BRL,  Grand  Island,  NY). 

All  ER  PCRs  consisted  of  30  cycles  (30  sec  at  94°C,  30  sec  at  60°C,  and  30  sec  at  72°C).  PCR  products  were 
then  separated  on  6%  polyacrylamide  gels  containing  7M  urea.  Following  electrophoresis,  the  gels  were  dried 
and  autoradiographed.  Three  independent  PCRs  were  performed. 

Quantification  of  SRA  and  ER  RNA  expression.  Exposed  screens  were  scanned  using  a  Molecular 
ImagerTM-FX  (Bio-Rad,  Hercules,  CA)  and  the  intensity  of  the  signal  corresponding  to  SRA  or  the 
appropriate  ER  isoform  PCR  fragments  was  measured  using  Quantity  One™  software  (Bio-Rad,  Hercules, 
CA).  Three  independent  PCRs  were  performed.  In  order  to  control  for  variations  between  experiments,  a 
value  of  1  was  arbitrarily  assigned  to  the  signal  of  one  particular  tumor  measured  in  each  set  of  PCR 
experiments  (always  the  same  tumor)  and  all  signals  were  expressed  relative  to  this  signal.  Levels  of  SRA 
was  expressed  relative  to  REA  (SRA/REA),  AIBl  (SRA/AIBl)  or  ERa  (SRA/ERa)  in  each  individual  tumor 
sample.  Levels  of  ERp  isoforms  were  expressed  relative  to  other  ERp  isoforms  shown  imder  statistical 
analysis  and  as  previously  described  (19). 

Quantification  of  the  relative  expression  of  the  deleted  SRA  variant  RNA.  It  has  previously  been 
shown  that  the  co-amplification  of  a  wild-type  and  a  deleted  variant  SRA  (SRA-Del)  cDNA  resulted  in  the 
amplification  of  two  PCR  products,  the  relative  signal  intensity  of  which  provided  a  reliable  measurement  of 
the  relative  expression  of  the  deleted  variant  (16).  For  each  sample,  the  signal  corresponding  to  the  SRA-Del 
was  measured  using  Quantity  One™  software  (Bio-Rad,  Hercules,  CA)  and  expressed  as  a  percentage  of  the 
corresponding  core  SRA  signal.  For  each  case,  3  independent  assays  were  performed  and  the  mean 
determined. 

Quantification  of  REA  and  AIBl  RNA  expression.  Following  analysis  of  PCR  products  on  prestained 
agarose  gels,  signals  were  quantified  by  scanning  using  MultiAnalyst™  (Bio-Rad,  Hercules,  CA).  At  least, 
three  independent  PCRs  were  performed.  A  value  of  1  was  arbitrarily  assigned  to  the  REA  or  AIBl  signal  of 
the  one  particular  tumor  and  is  the  same  tumor  as  described  above  and  all  signals  were  expressed  relative  to 


this  signal.  Levels  of  AIBl  were  expressed  relative  to  REA  (ABl/REA)  and  ERa  (AIBl/ERa),  and  levels  of 
REA  were  expressed  relative  to  ERa  (REA/ERa ). 

Statistical  analysis.  Differences  between  tamoxifen  sensitive  and  tamoxifen  resistant  cases  were  tested 
using  the  Mann- Whitney  rank  sum  test,  two-tailed.  Potential  differences  in  expression  between  the  two 
groups  with  respect  to  each  ERp  isoform  RNA  relative  to  other  ERp  isoform  RNA  expression  (e.g.  log 
ERp  1/total  ERp;log  ERp2/total  ER|3;log  ERp5/total  ERp,  as  previously  described  (19))  ,  and  the  relative 
expressions  of  co-regulators  (i.e.  logAIBl/REA;  logSRA/REA;  logSRA/AIBl;  logAIBl/ERa; 
logSRA/ERa ;  logREA/ERa )  were  determined. 

Tumors  were  classified  as  low  (scores  between  0  and  100)  and  high  (150  to  300)  ERp  expressors,  and 
differences  between  tamoxifen  sensitive  and  tamoxifen  resistant  cases  were  tested  using  Fisher’s  exact  test. 
Correlation  between  ERp  protein  expression  (H-score)  and  tumor  characteristics  was  tested  by  calculation  of 
the  Spearman  coefficient  r. 

Results. 

Expression  of  ERp  protein  in  primary  human  breast  tumors  from  patients  who  later  progressed  on 
tamoxifen  treatment  or  showed  no  progression  on  tamoxifen  treatment.  ERp  protein  was  determined 
immunohistochemically  on  adjacent  sections  firom  each  tumor,  using  two  different  antibodies.  GC-17  is  an 
antibody  recognizing  an  epitope  in  the  C-terminus  of  full-length  ERpl  (13).  Normal  breast  tissue  was  used  as 
a  positive  control  and  is  shown  in  Figure  1  A.  Examples  of  staining  in  hxmian  breast  tumor  sections  are  shown 
in  Figure  IB  to  D.  Some  tumor  sections  showed  no  (Figure  IB,  wild-type  ERp  score  =  0)  or  low  (Figure  1C, 

wild-type  ERp  score  =  100),  while  others  showed  strong  wild-type  ERp  signals  (Figure  ID,  wild-type  ERp 
score  =  300).  Tumors  were  classified  as  low  (scores  between  0  and  100)  and  high  (150  to  300)  wild-type 
ERp  protein  expressors,  and  differences  between  between  tamoxifen  sensitive  and  resistant  tumors 
determined  by  Fishers  exact  test.  No  significant  differences  were  found. 

IgYERbeta503  is  an  antibody  that  recognizes  all  known  putative  ERp  protein  isoforms  (11,  12)  and 
which  we  refer  to  as  total  ERp  protein.  Normal  breast  tissue  was  used  as  a  positive  control  and  is  shown  in 
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Figure  2A.  Examples  of  staining  with  this  antibody  in  human  breast  tumor  sections  are  shown  in  Figure  2  B  to 
D.  Some  sections  showed  no  (Figure  2B,  total  ERp  score  =  0)  or  low  (Figure  2C,  total  ERp  score  =  100)  total 
ERp  expression  whereas  others  had  strong  total  ERp  protein  signal  (Figure  2D,  total  ERp  score  =  300). 
Tumors  were  classified  as  low  and  high  total  ERp  protein  expressors,  and  there  was  a  statistically  significant 
difference  in  high  total  ERp  protein  between  the  Tamoxifen  sensitive  and  resistant  groups  (Fishers  exact  test, 
P=0.046).  High  total  ERp  protein  expressors  were  more  fi-equently  observed  in  tamoxifen  sensitive  tumors 
than  resistant  tumors. 

Correlation  between  ERp  protein  expression  and  tumor  characteristics  was  tested  by  calculation  of  the 
Spearman  coefficient.  A  positive  correlation  between  ERpl  (GC17)  protein  and  PR  levels  (Spearman  r=0.44, 
P=0.022)  was  foimd  when  each  was  examined  as  continuous  variables.  When  tumors  were  divided  into  PR+ 
(>  10  finol/mg  protein)  or  PR-  (<  10  finol/mg  protein)  groups  was  used  there  was  a  significantly  higher  level 
of  ERpl  (GC17)  protein  in  PR+  tumors  compared  to  PR-  tumors  (Mann- Whitney  test,  P  =  0.0268;  median 
for  PR+  tumors  =  55,  range  5  to  150  and  median  for  PR-  tumors  =10,  range  0  to  75).  As  well,  there  was  also 
a  significantly  higher  level  of  total  ERp  (IgY503)  protein  in  PR+  tumors  compared  to  PR-  tumors  (Mann- 
Whitney  test,  P  =  0.0085;  median  for  PR+  tumors  =  125,  range  25  to  270  and  median  for  PR-  tumors  =  50 
range  0  to  100). 

Relative  expression  of  ERp  isoform  RNA  in  primary  human  breast  tumors  from  patients  who  later 
progressed  on  tamoxifen  treatment  or  showed  no  progression  on  adjuvant  tamoxifen.  To  determine 
if  the  differences  described  above  in  ERp  protein  expression  were  correlated  with  differences  in  ERp  variant 
isoform  RNA  expression,  we  compared  the  relative  expression  of  ERp  variant  RNA  to  wild-type  ERp  RNA  in 
the  tamoxifen  sensitive  and  resistant  groups.  Unfortunately,  frozen  tissue  samples  corresponding  to  some  of 
the  paraffin  blocks  from  patients  in  the  cohort  used  for  immunohistochemistry  were  not  available.  The  tumor 
bank  was  queried  and  additional  cases  selected  as  described  in  Materials  and  Methods.  Using  previously 
validated  assays  (7,  19)  the  relative  expression  of  ERp2,  ERp5  and  wild-type  ERpl  RNA  in  the  tamoxifen 
sensitive  and  resistant  groups  was  not  significantly  different. 
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Relative  expression  of  coregulators  in  primary  human  breast  tumors  from  patients  who  later 
progressed  on  tamoxifen  treatment  or  showed  no  progression  on  tamoxifen  treatment.  To  address  the 
hypothesis  that  altered  relative  expression  of  steroid  receptor  coactivators  and  corepressors  could  underly 
altered  tamoxifen  sensitivity  in  human  breast  tumors,  and  since  we  previously  showed  that  the  relative 
expression  of  two  coactivators  (SRA  and  AIBl)  to  a  corepressor  (REA)  is  altered  in  ER+  breast  tumors 
compared  their  adjacent  normal  breast  tissue,  we  chose  these  coregulators  to  study.  They  were  measured  by 
RT-PCR  in  the  above  tumor  cohorts.  SRA,  AIBl,  and  REA  mRNAs  were  detectable  in  most  samples,  even 
though  their  level  of  expression  differed  from  one  sample  to  another.  Consistent  with  our  previous  results 
(16),  an  additional  fragment,  migrating  at  an  apparent  size  of  459  bp  was  also  observed  in  most  tumors  when 
using  SRA  specific  primers.  Sequencing  analysis  revealed  that  this  band  corresponded  to  a  variant  form  of 
SRA  (referred  to  as  SRA-Del)  deleted  in  203  bp  between  positions  155  and  357  (positions  givai  correspond 
to  Genbank  accession  number  AF092038).  There  were  no  significant  differences  between  the  tamoxifen 
sensitive  and  the  de  novo  tamoxifen  resistant  breast  cancers  in  the  relative  expression  of  any  of  the 
coactivators  to  corepressor  RNA,  or  in  the  relative  expression  of  SRA/AIBl  RNA,  or  in  expression  of  any  of 
these  coregulator  RNAs  relative  to  ERa  or  total  ERp  RNA  expression.  As  well,  there  was  no  significant 
difference  in  the  relative  expression  of  variant  SRA/wild-type  SRA  between  the  groups  either. 

Tumor  Characteristics.  No  statistically  significant  differences  were  found  between  the  tamoxifen  sensitive  and 
tamoxifen  resistant  cohorts  in  any  of  the  tumor  characteristics  described  in  the  Materials  and  methods  section 
except  for  PR.  PR  levels  were  statistically  significantly  different  (P=0.044)  between  the  two  groups  using  a 
Mann- Whitney  rank  sum  test  (two  sided).  PR  levels  were  higher  (median  PR  was  32  fmol/mg  protein;  range  8  - 
216  finol/mg  protein)  in  the  tamoxifen  ‘sensitive’  group  compared  to  the  tamoxifen  ^de  novo  resistant’  group 
(median  PR  was  14  finol/mg  protein;  range  4  -  288  finol/mg  protein).  This  was  a  consistent  finding  in  both 
selected  cohorts  (that  used  for  immunohistochemistiy  and  that  used  for  the  RNA  study),  and  provides  strong 
support  for  differences  in  estrogen  signaling  pathways  in  these  two  groups  since  PR  is  a  marker  of  ER  signal 
transduction  (20, 21). 

Discussion 


We  and  others  have  shown  that  the  relative  expression  of  ERa  and  ERp  is  significantly  altered  during 
breast  tumorigenesis  (7, 22),  and  a  similar  mechanism  has  been  proposed  to  underlie  tamoxifen  resistance  in 
breast  cancers  (23).  The  current  study  shows  no  significant  differences  in  expression  of  wild-lype 
ERp  (ERpl)  between  tamoxifen  sensitive  and  resistant  tumors.  Interestingly,  in  this  small  cohort  of  tumors 
when  total  ERp  expression  was  examined,  there  were  significantly  more  high  total  ERp  expressors  in  the 
tamoxifen  ‘sensitive’  compared  to  the  ‘resistant’  group.  The  data  suggest  the  possibility  that  increased  and 
altered  ERp  isoform  protein  expression  may  have  a  role  in  de  novo  tamoxifen  resistance,  or  at  least  together 
with  other  parameters  may  provide  better  markers  of  endocrine  sensitivity.  The  increased  expression  of  ERp 
proteins  in  the  tamoxifen  sensitive  group  is  also  consistent  with  recently  published  data  where  patients  with 
ERp  positive  tumors  (determined  using  an  antibody  to  an  N-terminal  epitope  of  the  ERp  protein,  and  defined 
as  nuclear  staining  in  >  10%  of  cancer  cells)  had  a  significantly  better  overall  survival  than  patients  with  ERp 
negative  tumors  while  receiving  adjuvant  tamoxifen  therapy  (24).  Both  these  latter  data  and  those  presented 
currently  in  this  manuscript  are  in  contrast  to  data  showing  increased  ERp  RNA  expression  in  tamoxifen 
resistant  tumors  versus  tamoxifen  sensitive  tumors  previously  published  (25).  Together  these  studies  suggest 
that  the  ERp  status  and  the  nature  of  ERp  isoforms  together  with  ERa  status  in  human  breast  cancers  may  be 
important  biomarkers  on  endocrine  sensitivity,  and  warrants  further  study,  in  larger,  clinical  trials  cohorts.  The 
association  of  increased  ERp  isoform  expression  with  tamoxifen  sensitivity,  suggests  a  possible  mechanistic 
role,  and  one  possible  mechanism  may  be  suggested  by  several  publications  which  have  shown  that 
ERp  isoforms  have  a  modulatory  effect  on  ERa,  both  in  normal  tissues  (26)  as  well  as  in  cell  culture  models 
(6,27). 

Interestingly,  if  one  examines  previous  studies  in  which  ERp  protein  expression  was  investigated  in 
breast  tumors  immunohistochermcally,  marked  differences  in  correlations  of  ERp  protein  with  various  clinical 
and  pathological  parameters  have  been  documented  (24, 28-30).  Much  of  the  difference  seen,  can  be  related  to 
the  antibody  used  and  often  if  only  the  wild-type  ERp  was  detected  or  if  total  (wild-type  plus  variant) 
ERp  isoforms  were  likely  detected.  In  our  hands  measurement  of  wild-type  ERp  protein  expression  was 
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positively  correlated  with  PR  levels  when  measured  as  continuous  variables  as  well  as  when  a  cut-off  value  for 
PR  positivity  was  used.  Total  ERp  protein  expression  was  correlated  with  PR  levels  only  when  a  cut-off  values 
for  PR  positivity  was  used.  This  is  a  reasonably  consistent  finding  in  the  literature  usually  when  only  wild- 
type  ERp  is  measmed  (24,  28,  30,  31),  although  recently  the  same  relationship  has  been  observed  using  an 
antibody  recogmzing  an  N-terminal  epitope  and  therefore  total  ERp  protein  expression  (30),  as  well  as  in  the 
current  study.  Indeed  most  of  these  studies  also  find  that  wild-type  ERp  expression  is  correlated  with  both 
ER  and  PR.  The  current  study  did  not  find  a  correlation  with  ER  levels,  which  is  likely  due  to  the  selection 
criteria  for  our  tumor  cohorts,  retrospectively,  to  be  node  negative  and  ER+  tumors  and  then  to  have  either  no 
disease  recurrence  while  on  tamoxifen  treatment  or  to  have  died  or  had  disease  recurrence/progression  while 
on  tamoxifen  treatment.  Indeed  the  tumors  in  the  two  cohorts  had  no  significant  differences  in  level  of  ER 
expression  as  determined  by  ligand  binding.  The  data  in  this  current  study  are  in  contrast  with  our  previous 
observations  (32)  tiiat  total  ERp  RNA  expression  is  inversely  correlated  with  PR  expression  and  progestins 
downregulate  total  ERp  RNA  levels  (32).  However,  our  previous  data  ware  measured  in  tumors  that  ranged 

widely  in  ER/PR  status  and  included  ER-/PR+  and  ER-/PR-  tumors  (32),  whereas  all  the  tumors  in  the 
current  study  are  ER+,  the  majority  are  PR+,  as  well  and  none  are  ER-/PR-. 

The  potential  differences  between  tamoxifen  sensitive  and  resistant  groups  with  respect  to  ERp-like 
proteins,  was  not  correlated  with  differences  in  the  relative  expression  of  wild  t3Ape  ERp  and  two  known 
variants  ERp2  and  ERp5  at  the  RNA  level  between  the  tamoxifen  ‘sensitive’  versus  the  tamoxifen  ‘resistant’ 
groups,  however.  This  may  due  to  differential  regulation  of  protein  versus  RNA  level  or  the  likelihood  that 
there  are  other  potential  ERp  isoforms  (known  and  unknown)  expressed  in  breast  tissues  in  addition  to 
ERpl,  ERp2  and  ERp5  (33, 34),  whose  cognate  proteins  would  be  detected  by  the  antibody  but  not  measure  in 
the  triple  primer  RT-PCR  assay. 

Another  mechanism  for  differential  tamoxifen  sensitivity  in  ER+  breast  tumors  could  be  altered 
coregulator  expression.  Since  we  previously  had  shown  differential  expression  of  some  coregulators  during 
breast  tumorigenesis  in  ER+  tumors,  a  semi-quantitative  reverse-transcription  polymerase  chain  reaction  (RT- 


PCR)  assay  was  used  to  detennine  the  relative  expression  of  ER  co-regulators  SRA,  AIBl  and  REA  in  the 
primary  breast  tumors  of  patients  who,  following  tamoxifen  adjuvant  therapy,  relapsed  or  remained  disease 
free.  To  our  knowledge,  no  study  has  previously  compared  the  expression  of  these  co-regulators  (SRA,  AIBl 
and  REA)  in  such  cohorts.  No  significant  differences  were  observed  in  the  ratios  of  any  of  the  coactivators 
/corepressors  or  any  of  the  ratios  of  these  coregulators  to  ERaRNA  levels  between  the  two  breast  tumor 
groups,  from  patients  who  were  later  formd  to  be  disease  free  (sensitive)  or  have  disease  progression 
(resistant)  while  on  adjuvant  tamoxifen  treatment.  These  data  suggest  that  altered  relative  expression  of  these 
coregulators  is  unlikely  to  be  a  marker  of  tamoxifen  sensitivity  in  ER+,  node  negative,  primary  breast  tumors, 
and  unlikely  to  have  a  functional  role  in  de  novo  tamoxifen  resistance.  Although  SRA  is  functional  as  an 
RNA  molecule,  REA  and  AIBl  are  functional  as  proteins.  Furthermore,  other  factors  can  affect  protein 
activity  for  example  phosphorylation  in  the  case  of  AIBl  (35)  or  sequestration  by  other  proteins  such  as 
prothymosin-alpha  in  the  case  of  REA  (36).  Our  studies  do  not  exclude  differences  at  the  protein  and/or 
activity  levels  of  REA  and  AIBl  being  involved  in  de  novo  tamoxifen  resistance,  nor  do  they  exclude  altered 
expression  of  these  factors  having  a  role  in  acquired  tamoxifen  resistance.  The  possible  involvement  of 
altered  coregulators  in  acquired  resistance  has  support  from  previously  published  data  where  it  was 
demonstrated  that  decreased  levels  of  the  nuclear  receptor  co-repressor  (N-CoR)  complexes,  that  play  a  role  in 
the  mediation  of  the  inhibitory  effects  of  tamoxifen,  correlate  with  the  acquisition  of  tamoxifen  resistance  in  a 
mouse  model  system  of  acquired  tamoxifen  resistance  for  human  breast  cancer  (37).  One  should  also  note 
that  Bems  et  al.  (38)  recently  reported  that  a  high  expression  of  SRC-1  co-activator  indicated  a  favorable 
response  to  tamoxifen  of  patients  with  recurrent  breast  cancer,  and  more  recently  Chan  et  al.(39)  reported  no 
significant  differences  compared  to  control  groups  in  expression  of  the  coactivators  TIF-1,  SUG-1,  RIP  140 
and  the  corepressor  SMRT  RNAs  in  ER+  human  breast  tumors  from  patients  that  were  de  novo  tamoxifen 
resistant  or  in  recurrent  tumors  obtained  from  patients  who  acquired  tamoxifen  resistance.  Altogether,  there  is 
little  evidence  for  altered  coregulators  expression  in  breast  tumors  that  are  de  novo  tamoxifen  resistant. 
However,  our  data  provide  preliminary  evidence  that  the  expression  of  ERp  protein  isoforms  may  differ  in 
primary  tumors  of  breast  cancer  patients  who  prove  to  have  differential  sensitivity  to  tamoxifen  therapy.  As 


well  oxir  data  support  distinct  differences  in  the  ER  signaling  pathways  between  these  two  groups  of  patients 
since  the  expression  of  a  known  estrogen  responsive  gene  PR  is  significantly  different  between  the  two 
groups,  the  precise  mechanisms  underlying  these  differences  remain  to  be  elucidated. 


Legends  to  Figures. 

Figure  1.  Examples  of  immunohistochemistty  using  the  GC-17  antibody  which  will  only  recognized 
the  wild-type  ERpl:  A,  normal  human  breast  tissue;  B,  ERpl  negative  human  breast  tumor,  H- 
score  =  0;  C,  ERpl  low  expressing  human  breast  tumor,  H-score  =  100;  D,  ERpl  high  expressing 
human  breast  tumor,  H-score  =150. 

Figure  2.  Examples  of  immxmohistochemistry  using  the  IgYERpSOS  antibody  which  recognizes 
most  ERp  isoforms;  A,  normal  human  breast  tissue;  B,  ERp  negative  human  breast  tumor,  H-score 
=  25;  C,  ERp  low  expressing  human  breast  tumor,  H-score  =  100;  D,  ERp  high  expressing  human 
breast  tumor,  H-score  =  225. 
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Abstract 

Recently,  Lanz  et  al.  (Cell  97,  17-27,  1999)  identified  and  characterized  a  new  transcriptional  co-activator,  called 
steroid  receptor  RNA  activator  (SRA),  which  was  found  to  fimction  as  an  RNA  transcript  and  to  be  specific  for 
steroid  receptors.  Lanz  et  al.  reported  that  none  of  the  different  SRA  transcript  isoforms  isolated  encoded 
detectable  levels  of  SRA  protein  following  in  vitro  translation  experiments.  We  have  identified  three  new  SRA 
isoforms  encoding  in  vitro  stable  SRA  proteins.  One  of  these  isoforms  contains  a  point  mutation  followed  by  a 
codon  insertion  that  is  found  in  both  normal  and  tumor  breast  tissue  as  well  as  in  some  breast  tumor  cell  lines. 
This  insertion  is  present  at  the  DNA  level  and  is  likely  to  represent  a  genetic  polymorphism.  We  hypothesize 
that  wild-type  and  inserted  SRA  mRNA  and  protein  could  play  a  role  in  the  mechanisms  underlying  breast 
tumorigenesis  and  tumor  progression. 


Key  words:  SRA,  steroid  receptor  co-activator,  normal  and  tumor  human  breast  tissue,  polymorphism, 
tumorigenesis,  PCR. 

Introduction 

Endogenous  hormones  such  as  estrogen  and  progesterone,  known  to  regulate  the  growth  and  the  development  of 
female  reproductive  organs,  also  play  key  roles  in  other  target  tissues  such  as  bone,  central  nervous  system  and 
cardiovascular  systems.  Estrogen  and  progesterone  actions  are  mainly  mediated  through  specific  receptors  that 
belong  to  the  steroid/thyroid/retinoic  acid  receptor  super-family  and  act  as  ligand-dependent  transcription  factors 
(1-3).  The  transcription  of  target  genes  by  hormone-liganded  receptors  depends  upon  interactions  between  these 
receptors  and  several  members  of  a  complex  co-regulator  population  (see  2  and  references  herein).  To  an  already 
long  list  of  nuclear  receptor  co-activators  (3),  Lanz  et  al.  (4)  recently  added  the  steroid  receptor  RNA  activator 
(SRA).2  SRA  differs  from  other  co-activators  in  two  main  features.  Firstly,  SRA  transcripts  do  not  appear  to  be 
translated  and  therefore,  this  co-activator  acts  as  an  RNA  and  not  as  a  protein.  Indeed,  Lanz  et  al.  were 
unsuccessful  when  trying  to  generate  in  vitro  or  in  vivo  stable  SRA  encoded  protein  and  showed  that  SRA  exists 
in  a  ribonucleoprotein  complex  that  contains  SRC-1  and  is  recruited  by  steroid  receptors  (4).  Secondly,  SRA 
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appears  to  be  actually  specific  for  steroid  receptors.  Indeed,  most  of  the  receptor-interacting  factors  interact  with 
and  co-activate  both  class  I  and  class  II  nuclear  recqitors  (4). 

Our  recent  data  suggest  that  the  expression  of  SRA  and  AIBl  (another  steroid  receptor  co-activator)  is  higher  in 
estrogen  receptor  (ER)  positive  breast  tumors  than  in  matched  adjacent  normal  breast  tissues  (5  and  unpublished 
data).  As  well,  we  showed  that  SRA  expression  correlated  with  estrogen  recq)tor  and  progesterone  receptor  levels 
in  particular  tumor  subgroups  (6).  A  variant  form  of  SRA,  containing  a  deletion  of  203  bases  wittun  the  SRA 
core  sequence,  was  also  observed  in  breast  tumor  tissues.  The  relative  expression  of  this  new  SRA  isoform 
correlated  with  tumor  grade  (6).  Altogether,  these  data  suggested  that  changes  in  the  expression  of  SRA  related 
molecules  occur  during  breast  tumorigenesis  and  breast  tumor  progression.  We  therefore  speculated  that 
changes  in  SRA  expression  could  be  involved  in  the  mechanisms  underlying  these  phenomena. 

SRA  sequence  (Genbank  AF092038)  is  fully  contained  in  a  genomic  sequence  of  chromosome  5  (Genbank 
AC005214),  within  five  separated  exon-like  regions  (FIG.l).  Within  this  latter  clone  we  have  found,  using  Gene 
finder  (http://dot.imgen.bcm.tmc.edu:9331/gene-finder/gfhtml),  the  sequence  of  a  hypothetical  mRNA  which 
could  encode  a  236  amino-acid  protein,  that  corresponds  exactly  to  the  published  SRA  sequence,  except  for  an 
additional  37  nucleotides  in  the  5'  region  (FIG.l).  This  5’  region  now  contains  an  AUG  codon  encoding  the  first 
methionine  of  a  novel,  putative  236  amino-acid  SRA  protein  in  contrast  to  the  previously  predicted  162  amino- 
acids  which  was  unstable.  We  have  investigated  the  expression  of  this  hypothetical  mRNA  in  normal  and  tumor 
breast  tissues. 

Materials  and  methods 

Human  breast  tissues  and  cell  lines. 

14  cases  were  selected  fi’om  the  NCIC-Manitoba  Breast  Tumor  Bank  (Winnipeg,  Manitoba,  Canada).  For  each 
case,  matched  adjacent  normal  and  tumor  fi'ozen  tissue  blocks  were  available.  The  quality  of  each  block  was 
determined  by  the  histo-pathological  assessment  of  sections  from  adjacent  mirror  image  paraffin  embedded 
tissue  blocks,  as  previously  described  (7).  The  presence  of  normal  ducts  and  lobules  as  well  as  the  absence  of 
any  atypical  lesion  were  confirmed  in  all  normal  tissue  specimens.  Breast  epithelial  cell  lines  (BT-20,  MDA- 
MB-468,  MDA-MB-231,  MCFIOATI,  MCFIOAI,  MCF10AT3B,  ZR-75,  T47D,  T5,  MCF7,  HBLIOO)  were 
grown,  harvested  and  cell  pellets  stored  at  -70°C,  as  previously  described  (8).  Total  RNA  was  extracted  from  20 
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(im  frozen  tissue  sections  (15  and  5  sections  for  normal  and  tumor  breast  tissue,  respectively)  or  cell  pellets 
using  Trizol™  reagent  (Gibco  BRL,  Grand  Island,  NY)  according  to  the  manufacturer's  instructions. 

Primers  and  RT-PCR  conditions. 

Detection  of  the  hypothetical  coding  SRA  RNA\  primers  used  consisted  of  SRAUl  primer  (5'- 
TCCTTTGGTGCCTTGTGAC-3';  sense;  position  36132-36149,  Genbank  accession  number  AC005214),  and 
SRALl  primer  (5'-AGTCTGGGGAACCGAGGAT-3';  antisense;  position  43128-43110,  Genbank  accession 
number  AC0052 14). 

Pilot  study:  200  ng  of  total  RNA  was  reverse  transcribed  in  a  final  volume  of  25  |li1  using  ThermoScript^''^  RT- 
PCR  system  (Gibco  BRL,  Grand  Island,  NY)  and  oligo-dT  according  to  the  manufacturer's  instructions.  PCR 
amplifications  were  performed  as  previously  described  (10).  Briefly,  1  pi  of  RT  mixture  was  amplified  in  a 
final  volume  of  30  pi,  in  the  presence  of  20  mM  Tris-HCl  (pH  8.4),  50  mM  KCl,  2  mM  MgCl2,  0.2  mM 
dATP,  0.2  mM  dTTP,  0.2  mM  dGTP,  0.2  mM  dCTP,  4  ng/pl  of  each  primer  and  1  unit  of  Taq  DNA 
polymerase  (GIBCO-BRL).  Each  PCR  consisted  of  a  2  minute  pre-incubation  step  at  94°C  followed  by  30 
cycles  of  amplification  (30  sec  at  94°C,  2  minutes  at  72°C  and  30  seconds  at  55  °C)  using  a  Thermocycler 
(Perkin  Elmer).  PCR  products  were  then  separated  on  2%  agarose  gels  before  staining  with  ethidium  bromide 
(15  pg/ml).  PCR  products  were  sub-cloned  using  TOPO  TA  cloning  ®  kit  for  sequencing  (Invitrogen, 
Carlsbad,  CA)  according  to  the  manufacturer's  instruction  and  sequenced  as  previously  described  (10). 
Optimized  conditions:  With  further  experimentation,  we  found  that  the  above  conditions  w^e  not  optimal:  we 
used  the  PCR  Optimizer™  kit  (Invitrogen,  Carlsbad,  CA)  to  set-up  optimal  RT-PCR  amplification  conditions. 
One  pg  of  total  RNA  was  reverse  transcribed  in  a  final  volume  of  25  pi  using  Moloney  Murine  Leukemia 
Virus  (M-MLV)  reverse  transcriptase  and  random  hexamers  as  previously  described  (10).  One  pi  of  RT 
mixture  was  amplified  in  a  final  volume  of  30  pi,  in  the  presence  of  60  mM  Tris-HCl  (pH  8.5),  15  mM 
[NH4]'S04,  1.5  mM  MgCl2,  0.2  mM  dATP,  0.2  mM  dTTP,  0.2  mM  dGTP,  0.2  mM  dCTP,  4  ng/pl  of  each 
primer  and  1  unit  of  Taq  DNA  polymerase  (GIBCO-BRL).  Each  PCR  consisted  of  a  5  minutes  pre-incubation 
step  at  94°C  followed  by  30  cycles  of  amplification  (30  sec  at  94°C,  30  seconds  at  55  °C,  and  2  minutes  at 
72°C). 
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Detection  of  SRA  isoform  3  isoform:  primers  used  consisted  of  SRAU2  primer  (5- 
GGGCCTCCACCTCCTTCAAGTA-3';  sense;  position  41664-41685,  Genbank  accession  number 
AC005214),  SRAL2  primer  (5'-GCAGTCTTCCAATGCCTG-3';  antisense;  position  41813-41796,  Genbank 
accession  number  AC005214)  and  SRAL-GTCG  primer  (5'-CACATCCTCCATCAGTCG-3’;  antisense; 
position  41780-41767,  Genbank  accession  number  AC005214  plus  GTCG  sequence).  Radioactive  PCR 
amplifications  were  performed  and  PCR  products  analyzed  as  previously  described  (11)  with  minor 
modifications.  Briefly,  1  pi  of  reverse  transcription  mixture  was  amplified  in  a  final  volume  of  15  pi,  in  the 

presence  of  1.5  p.Ci  of  (a-^^P)  dCTP  (3000  Ci/mmol),  4  ng/pl  of  each  primer  (SRAU2/SRAL2  or  SRAU2/ 
SRAL-GTCG  primers  )  and  0.3  unit  of  Taq  DNA  polymerase  (Gibco  BRL,  Grand  Island,  NY).  Each  PCR 
consisted  of  a  2  minutes  pre-incubation  step  at  94°C  followed  by  30  cycles  (30  sec  at  60®C,  30  sec  at  72°C  and 
30  sec  at  94°C).  PCR  products  were  then  separated  on  6%  polyacrylamide  gels  containing  7M  urea.  Following 
electrophoresis,  the  gels  were  dried  and  exposed  30  minutes  to  a  Molecular  Imager™-FX  Imaging  screen  (Bio- 
Rad,  Hercules,  CA).  Exposed  screen  was  then  scanned  using  a  Molecular  Imager™-FX  (Bio-Rad,  Hercules, 
CA).  As  positive  control,  we  amplified  10  ng  of  sequenced  plasmids  previously  shown  to  correspond  to 
SRAcodl,  SRAcod2  and  SRAcodl  inserted. 

In  vitro  synthesis  of  SRA  protein 

In  vitro  transcription/translation  reactions  were  performed  using  the  TnT  coupled  Reticulocyte  Lysate  System 
(Promega,  Madison,  El)  according  to  the  manufacturer's  instructions  and  using  SRAcodl,  SRAcod2  and 
SRAcodl  inserted  corresponding  plasmids.  The  programmed/im-programmed  lysates  were  then  subjected  to 
SDS-PAGE  and  [^^S]-methionine  labelled  protein  bands  visualized  after  exposing  overnight  to  a  Molecular 
Imager™-FX  Imaging  screen  (Bio-Rad,  Hercules,  CA)  subsequently  scanned  using  a  Molecular  Imager™-FX 
(Bio-Rad,  Hercules,  CA). 

Results 

Detection  of  hypothetical  coding  SRA  isoforms  in  normal  and  tumor  breast  tissues. 

In  order  to  determine  if  indeed  our  hypothetical  coding  SRA  mRNA  was  expressed  in  vivo,  primers  were 
designed  corresponding  to  sequences  upstream  of  the  putative  first  AUG  codon  (SRAUl,  contained  in  the 
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genomic  sequence  of  chromosome  5,  Genbank  AC005214)  and  dowstream  of  the  putative  stop  codon 
(SRALl).  See  Materials  and  Methods  and  FIG.  1  for  the  exact  primer  positions. 

Initially,  total  RNA  was  extracted  from  two  normal  and  one  tumor  breast  tissues.  Reverse  transcription  was 
performed  using  Thermoscript  and  PCR  amplification  performed  as  described  in  Materials  and  Methods  {pilot 
study).  Using  these  primers,  we  obtained  the  predicted  920  bp  PCR  product  in  both  normal  and  tumor 
components  (data  not  shown).  Cloning  and  sequencing  of  this  fragment  revealed  that  it  essentially  corresponded 
to  the  hypothetical  SRA  coding  sequence.  Interestingly,  three  different  SRA  cDNAs  were  identified:  SRA 
isoform  1  (Genbank  AF293024),  SRA  isoform  2  (Genbank  AF293025)  and  SRA  isoform  3  (Genbank 
AF293026).  These  all  migrated  with  a  similar  apparent  size  of  920  bp  (see  FIG.  2).  SRA  isoform  1  contained 
the  full  coding  hypothetical  SRA,  whereas  SRA  isoform  2  contained  two  point  mutations  at  positions  338 
(C->T)  and  348  (A^C),  and  SRA  isoform  3  contained  a  point  mutation  followed  by  an  insertion  of  three 
nucleotides  at  position  520  (G->CGAC).  All  these  sequences  contained  a  potential  open  reading  frame  able  to 
encode  two  236  aa  and  one  237  aa  protein  from  SRA  isoform  1, 2  and  3  respectively  (FIG.3). 

In  vitro  translation  of  three  new  SRA  isoforms:  SRA  isoforms  1,  2  and  3. 

Previously  cloned  SRA  cDNAs,  in  which  only  the  third  ATG  codon  (encoding  Met3,  FIG.2)  was  present,  were 
xmable  to  support  detectable  protein  synthesis  in  vitro  (4).  To  determine  if  the  three  new  SRA  cDNAs  isolated 
in  our  laboratory  could  be  translated  in  vitro,  vectors  containing  SRA  isoform  1, 2  and  3  sequences  downstream 
of  a  T7  polymerase  promoter  were  used  in  a  TnT  coupled  Reticulocyte  Lysate  system  as  described  in  Materials 
and  Methods.  Results  are  shown  FIG.  4.  The  three  different  SRA  isoforms  encoded  stable  SRA  proteins  under 
these  conditions.  It  should  however  be  noted  that  two  bands  were  observed  for  each  construction,  consistent 
with  the  use  of  two  different  initiating  methionine  codons  (FIG.2,  Metl  and  Met2).  Moreover,  the  observed 
molecular  masses  (31  and  30  kDa  for  SRA  isoforms  1-2  and  3,  respectively)  were  slightly  higher  that  those 
predicted  (25.7  and  25.8  kDa  for  SRA  isoform  1-2  and  3,  respectively).  Possible  reasons  for  such  the 
discrepancy  between  observed  and  theoretical  molecular  weight  are  under  investigation. 

Detection  of  hypothetically  coding  SRA  isoforms  in  breast  epithelial  cell  lines. 

To  determine  if  the  SRA  putatively  coding  sequences  (that  will  be  referred  to  as  hypothetical  coding  SRAs) 
were  expressed  in  breast  cancer  cells,  total  RNA  was  extracted  from  different  cancer  cell  lines  and  amplified  as 
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described  above.  It  should  be  noted  that  the  conditions  used  in  the  pilot  study  failed  to  give  reproducible  results, 
even  on  the  control  plasmid  DNA  preparations  (data  not  shown).  Optimized  conditions  were  therefore 
determined  (see  Materials  and  Methods),  and  used  to  obtain  the  results  presented  in  the  following  sections.  The 
predicted  PCR  product  corresponding  to  the  hypothetical  coding  SRAs  was  observed  in  all  cell  lines  (FIG.5A), 
confirming  that  these  transcripts  are  indeed  expressed  in  tumorigenic  and  non  tumorigenic  human  breast 
epithelial  cells. 

Blast  searches  of  human  EST  databases  (http://www.ncbi.nlm.nih.gov/blast/)  revealed  that  the  specific  sequence 
of  SRA  isoform  3  (i.e  G-^CGAC  at  position  520)  had  already  been  cloned  by  others  (Genbank  accession 
number:  AW954396,  AW957456,  AW630779,  AA305793,  AA410852,  AA353911).  This  confirmed  our  own 
data  resulting  firom  independent  RT-PCRs,  cloning  and  sequencing  and  underlined  that  the  SRA  isoform  3  was 
not  the  result  of  a  technical  artifact.  Therefore,  it  was  of  interest  to  investigate  the  expression  of  this  isoform  in 
breast  cancer  cells.  PCR  primers  were  designed  to  specifically  amplify  a  firagment  overlapping  the  putatively 
inserted  region  (SRAU2,  SRAL2)  or  to  specifically  anneal  with  the  inserted  sequence  (SRAL-GTCG)  (see 
FIG.2  for  primer  positions).  PCR  products  amplified  using  the  former  set  of  primers  (SRAU2,  SRAL2)  were 
expected  to  migrate  at  an  apparent  size  of  150  bp  and  153  bp  for  SRA  isoforms  1-2  and  3  cDNA,  respectively. 
Using  the  latter  set  of  primers  (SRAU2,  SRAL-GTGC),  a  PCR  product  117  bp  long  was  expected  only  in 
samples  expressing  SRA  isoform  3  mRNA.  Results  obtained  using  SRAU2/SRAL2  and  SRAU2/SRAL- 
GTGC  are  shown  FIG.5B  and  FIG.5C,  respectively.  Interestingly,  some  cell  lines  expressed  only  SRA  isoform 
3  specific  fragment  (MDA-MB-231,  MCF-7)  whereas  others  expressed  both  SRA  isoforms  1-2  and  3  (T47D, 
T5). 

DNA  extracted  from  these  cell  lines  was  also  amplified  using  SRAU2  and  SRAL2.  As  shown  in  FIG.5D  a 
perfect  correspondence  existed  between  SRA  isoform  expression  and  DNA  sequence,  suggesting  the  possible 
existence  of  a  genetic  polymorphism  and  therefore  of  at  least  two  alleles  of  the  SRA  gene. 

Detection  of  similar  SRA  isoforms  in  matched  normal  and  tumor  breast  tissues. 

Total  RNA  was  extracted  from  matched  normal  and  tumor  breast  tissues  from  twelve  different  patients,  reverse 
transcribed  and  amplified  using  the  different  pairs  of  primers  mentioned  earlier.  Results  corresponding  to  4 
different  patients  are  shown  FIG. 6.  Hypothetical  coding  SRAs  were  detectable  in  matched  normal  and  tumor 
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components  (FIG.6A).  Interestingly,  even  though  the  cohort  studied  was  too  small  to  look  for  statistical 
significance,  it  was  noticed  that  the  expression  of  the  hypothetical  coding  SRAs  was  higher  in  tumor  than  in 
matched  normal  breast  tissues  (data  not  shown).  Using  primers  overlapping  the  inserted  region  of  SRA  isoform 
3  (FIG.6B)  or  specifically  amplifying  this  isoform  (FIG.6C),  it  appeared  that,  as  observed  for  breast  epithelial 
cells  grown  in  culture,  some  patients  expressed  either  one  (patients  2-4)  or  two  (patient  1)  isoforms.  Further,  the 
pattern  of  expression  (either  one  or  two  isoforms  expressed)  was  conserved  between  matched  normal  and  tumor 
components. 

Discussion 

Three  new  SRA  isoforms  have  been  identified,  that  are  able  to  encode  stable  SRA  proteins  in  vitro.  Interestingly, 
using  “classical”  RT-PCR  conditions  (Cf.  pilot  study)  it  was  not  possible  to  reproducibly  amplify  the  PCR 
fragments  corresponding  to  these  isoforms.  The  reasons  for  this  lack  of  reproducibility  are  unknown  but  RNA 
or  cDNA  secondary  structures  might  be  involved.  Indeed,  changing  the  set  of  primers  during  the  PCR  was  not 
sufficient  to  overcome  this  problem  (our  unpublished  data)  suggesting  inherent  properties  of  these  particular 
nucleotide  sequences.  The  technical  difficulties  we  have  encountered  associated  with  the  reverse-transcription 
and/or  amplification  of  these  isoforms  together  with  the  possible  tissue  specificity  of  their  expression  could 
underlie  the  lack  of  detection  of  these  isoforms  in  previous  studies. 

Indeed,  database  searches  revealed  that  although  many  partial  SRA-like  sequences  isolated  firom  various  normal 
and  tumor  tissues  have  been  entered  in  the  human  EST  sequence  database,  only  one  appears  to  correspond  to  a 
full  length  hypothetical  coding  SRA.  This  particular  EST  (Genbank  R51952),  1500  bp  long,  was  isolated  fi’om 
normal  breast  tissue,  but  had  only  been  sequenced  in  to  its  first  460  5’  nucleotides.  This  partial  sequence 
contained  exactly  the  5’  terminal  extremity  of  the  hypothetical  SRA  mRNA  we  had  predicted  fi'om  tiie  genomic 
clone.  We  have  fully  sequenced  this  EST  and  found  that  it  indeed  corresponded  to  the  SRA  isoform  2  which  we 
had  isolated  firom  human  breast  tissue.  Also  of  interest  is  the  fact  fliat  the  Northern  blot  analysis  performed  by 
Lanz  et  al.  has  shown  an  apparently  less  abundant  SRA  transcript  of  similar  length  ~1500  bases,  in  addition  to 
the  0.7-0.8  transcript,  in  some  normal  tissues  (e.g.  muscle),  but  not  all  (4).  Such  data  suggest  that  the  expression 
of  the  hypothetical  coding  SRA  isoforms  may  be  tissue  specific. 
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Our  data  are  the  first  to  show  that  naturally  expressed  SRA  RNAs  can  generate  a  cDNA  from  which  a  detectable 
SRA  protein  is  generated  in  vitro.  Indeed,  Lanz  et  al.,  who  reported  the  identification  of  three  different  SRA 
isoforms,  concluded  that  their  cloned  SRA  cDNA  sequences  (i.e.  predicted  to  encode  a  162  a.a.  protein)  did  not 
encode  a  detectable  translation  product  (4).  However,  when  engineered  to  form  a  fusion  product  e.g.  with  GST, 
Gal4  or  HSV-thymidine  kinase  initiation  sequences  attached  to  the  N-terminal  region,  the  appropriate  fusion 
SRA-like  product  was  detected.  This  suggests  the  possibility  that  the  extra  74  N-terminal  amino  acids  predicted 
in  the  hypothetical  coding  SRA  mRNAs  could  provide  sequences  which  would  result  in  a  translation  product 
detectable  in  vivo. 

We  observed  that  one  or  two  SRA  isoforms  can  be  expressed  in  breast  epithelial  cells  as  well  as  in  normal  and 
tumor  breast  tissues.  It  should  be  noted  that  our  PCR  assays  (using  SRAU2/SRAL2  and  SRAU2/SRAL- 
GTCG)  do  not  establish  whether  SRA  isoform  1  and/or  2  are  expressed.  Indeed,  these  assays  only  establish 
whether  or  not  SRA  isoform  3  is  expressed,  alone  or  together  with  one  of  the  other  isoforms  (SRA  isoform  1  or 
2).  Interestingly,  we  also  foimd  that  in  breast  cancer  cells  the  pattern  of  expression  was  directly  related  to  the 
detection  of  these  particular  isoforms  within  the  genomic  DNA.  Moreover,  when  normal  breast  tissue  expressed 
specific  isoforms,  the  adjacent  tumor  tissue  also  expressed  the  same  isoforms.  Altogether,  these  data  suggest  the 
existence  of  a  genetic  polymorphism  within  the  SRA  gene.  Moreover,  some  patients  appeared  to  be  homozygotes 
or  heterozygotes  for  this  gene.  Whether  or  not  this  polymorphism  contributes  to  a  predisposition  to  breast 
tumorigenesis  or  to  development  of  breast  tumors  with  a  particular  phenotype  (i.e  steroid  receptor  status, 
resistance  to  tamoxifen,  grade  etc...)  is  under  investigation. 

In  our  previous  studies,  we  showed  that  the  expression  of  SRA  (Genbank  AF092038)  correlated  with  estrogen 
receptor  and  progesterone  receptor  levels  in  particular  tumor  subgroups.  We  also  determined  that  the  expression 
of  SRA  was  higher  in  tumor  that  in  adjacent  normal  human  breast  tissues  (5  and  unpublished  observation).  It  is 
therefore  reasonable  to  speculate  that  the  different  SRA  isoforms,  at  the  RNA  and/or  at  the  protein  levels,  could 
be  involved  in  the  mechanisms  underlying  breast  tumorigenesis  and  tumor  progression.  The  role  of  the  different 
SRA  isoforms  and  their  putative  protein  products  on  ER  and  PR  signaling  pathway  remains  to  be  determined. 

In  conclusion,  we  have  identified  several  SRA  mRNA  isoforms,  which  are  able  to  encode  stable  proteins  in  vitro 
and  are  expressed  in  normal  and  neoplastic  breast  tissues.  We  therefore  speculate  that  SRA  proteins,  encoded 
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by  these  mRNA  isoforms,  are  present  in  vivo.  We  also  hypothesize  that  these  isoforms  and  their  translation 
products  could  be  involved  in  the  mechanisms  underlying  breast  tumorigenesis  and  tumor  progression. 


References 

1  Evans,  R.  M.  (1988).  The  steroid  and  thyroid  receptor  super-family.  Science  240, 889-895. 

2  Shibata,  H.,  Spencer,  T.  E.,  Onate,  S.  A.,  Jenster,  G.,  Tsai,  S.  Y,  Tsai,  M.  J.,  and  O'Malley,  B.  W.  (1997). 
Role  of  co-activators  and  co-repressors  in  the  mechanism  of  steroid/thyroid  receptor  action.  Recent  Prog. 
Horm.  Res.  52, 141-164. 

3  McKenna,  N.  J.,  Lanz,  R.  B.,  and  O'Malley,  B.  W.  (1999).  Nuclear  recq)tor  coregulators:  cellular  and 
molecular  biology.  Endocr.  Rev.  20, 321-344. 

4  Lanz,  R.  B.,  McKenna,  N.  J.,  Onate,  S.  A.,  Albrecht,  U.,  Wong,  J.,  Tsai,  S.  Y.,  Tsai,  M-J.,  and  O'Malley,  B. 
W.  (1999).  A  steroid  receptor  coactivator,  SRA,  functions  as  an  RNA  and  is  present  in  an  SRC-1  complex. 
Cell  97,  17-27. 

5  Murphy  L.  C.  and  Watson,  P.H.  (2000).  Estrogen  receptors  (ER)  in  human  breast  tumorigenesis.  Breast 
Cancer  Research:  Reasons  for  Hope.  Abstract  529. 

6  Leygue,  E.,  Dotzlaw,  H.,  Watson,  P.  H.,  and  Murphy,  L.  C.  (1999).  Expression  of  the  steroid  receptor 
RNA  activator  (SRA)  in  human  breast  tumors.  Cancer  Res,  59, 4190-4193. 

7  Hiller,  T.,  Snell,  L.,  and  Watson,  P.  (1996).  Microdissection/RT-PCR  analysis  of  gene  expression. 
Biotechniques  21,  38-44. 

8  Leygue,  E.,  Dotzlaw,  H.,  Watson,  P.  H.,  and  Murphy,  L.  C.  (1999).  Expression  of  estrogen  receptor  beta  1, 
beta  2  and  beta  5  mRNAs  in  human  breast  tissue.  Cancer  Res.  59,  1 175-1 179. 

9  Dotzlaw,  H.,  Leygue,  E.,  Watson,  P.  H.,  and  Murphy,  L.  C.  (1999).  Estrogen  receptor-_  mRNA  expression 
in  human  breast  tumor  biopsies:  relationship  to  steroid  receptor  status  and  regulation  by  progestins.  Cancer 
Res.  59, 529-532. 

10  Leygue,  E.,  Watson,  P.  H.,  and  Murphy,  L.  C.  (1999).  Estrogen  receptor  variants  in  normal  human 
mammary  tissue.  J.  Natl.  Cancer  Inst.  88, 284-290. 

1 1  Leygue,  E.,  Dotzlaw,  H.,  Watson,  P.  H.,  and  Murphy,  L.  C.  (1998).  Altered  estrogen  receptor  alpha  and  beta 
expression  during  human  breast  tumorigenesis.  Cancer  Res.  58, 3 1 97-320 1 . 


12 


^  et  al. 

Figures  legend 

FIG.  1  Schematic  representation  of  SRA  mRNA  and  gene  structure.  Human  SRA  isoforms  identified  to 
date  differ  in  their  5'  and  3'  terminal  regions  (white  circles)  but  present  an  identical  nucleotide  sequence  (dark 
gray  boxes)  in  between.  The  only  SRA  sequence  entered  in  Genbank  (AF092038)  is  fully  contained  within 
10,000  bp  of  Bac  5  genomic  sequence  (AC005214).  Numbers  indicate  positions  relative  to  this  latter  sequence. 
Hyp  coding  SRA:  hypothetical  coding  SRA  sequence  found  using  Gene  Finder.  Position  of  SRAUl  and 
SRALl  primers  are  indicated  by  arrows. 

FIG.  2  Alignment  of  SRA  isoforms  and  primer  positions.  Previously  (4)  cloned  SRA  sequence 
(AF092038)  is  aligned  with  new  SRA  isoforms:  SRA  1  (AF293024),  SRA  2  (AF293025)  and  SRA  3 
(AF293026).  Differences  between  sequences  are  circled.  Positions  of  two  new  putative  starting  ATG  codons  are 
indicated  (Met  1  and  Met  2),  together  with  that  of  the  initially  predicted  starting  codon  (Met  3,  reference  4)  and 
the  common  stop  codon.  The  positions  of  SRAUl,  SRAU2,  SRAL-GTCG,  SRAL2  and  SRALl  primers  are 
also  depicted. 

FIG.  3  Alignment  of  proteins  putatively  encoded  by  SRA  isoforms.  SRA,  SRA  1,  SRA  2,  and  SRA  3: 
Proteins  putatively  encoded  by  open  reading  frames  contained  in  AF092038,  AF293024,  AF293025  and 
AF293026,  respectively.  Positions  of  two  putative  first  methionine  (Met  1  and  Met  2),  together  with  initially 
predicted  first  methione  (Met  3,  reference  4)  are  indicated.  Differences  between  sequences  are  circled. 

FIG.4  In  vitro  translation  of  SRA  isoform  1,  2  and  3  cDNAs.  In  vitro  transcription/translation  reactions 
were  performed  using  SRA  isoform  1,  2  and  3  corresponding  plasmids  and  programmed  /  un-programmed 
(UPL)  lysates  analyzed  as  described  in  the  Materials  and  Methods  section.  Size  in  kDa,  corresponding  to  the 
prestained  marker,  broad  range  (Premixed  format,  BioLabs),  are  shown  on  the  left. 

FIG.5  Detection  of  SRA  isoforms  in  breast  epithelial  cell  lines.  Total  RNA  was  extracted  fi'om  breast 
epithelial  cell  lines  (BT-20,  MDA-MB-468,  MDA-MB-231,  MCFIOATI,  MCFIOAI,  MCF10AT3B,  HBL- 
100,  ZR-75,  T47D,  T5,  MCF-7),  reverse  transcribed,  PCR  amplified  using  SRAUl /SRALl  (A), 
SRAU2/SRAL2  (B)  or  SRAU2/SRAL-GTCG  (C)  primers  and  PCR  products  were  separated  on  agarose  (A, 


C)  or  aciylamide  (B)  gels,  as  described  in  the  Materials  and  Methods  section.  (D)  genomic  DNA 
corresponding  to  these  cell  lines  was  amplified  using  SRAU2/SRAL2  and  PCR  products  separated  on 
acrylamide  gel  as  described  in  the  Materials  and  Methods  section. 

FIG.6.  Detection  of  SRA  isoforms  in  matched  normal  and  tumor  breast  tissues.  Total  RNA  was 
extracted  fi'om  matched  normal  (N)  and  tumor  (T)  breast  tissues  from  different  patients  (1-4),  reverse 
transcribed,  amplified  using  SRAUl/SRALl  (A),  SRAU2/SRAL2  (B)  or  SRAU2/SRAL-GTCG  (C)  primers 
and  PCR  products  were  separated  on  agarose  (A)  or  acrylamide  (B,  C)  gels,  as  described  in  the  Materials  and 
Methods  section.  SRA  isol,  SRA  iso3,  and  SRA  isol/iso3:  positive  controls  amplified  in  parallel  consisting  of 
50  pg  of  plasmid  containing  SRA  isoform  1,  SRA  isoform  3  cDNA,  and  a  mix  of  these  plasmids  respectively. 


SR  A  I 
SR  A  2 
SR  A  3 


TCCTTTGCTCCCTTGTGAC 
ITCCTTTGGTGCCTTCTGAC 
Itcctttggtgccttg  T  G  a  C  I 


SRAUl 


CAGGGCCCTGATGGTTCATTAGATGCAGCCTTCGAGTCTTA 

CAGGGCCCTGATGGTTCATTAGATGGAGCCTTCGAGTCTTA 

CAGGGCCCTGATGGTTCATTACATGGAGCCTTCGAGTCTTA 


SRA  /  GCGAGTTGCCGCAGGGTCCCCACAGCCGCTCCCGACGGTTGTCAACCAGCATCCATCCTC 
SRA  2  CCGAGTTGCCGCAGGGTCCCCACAGCCGCTCCCGACGGTTGTGAACCAGCATCCATCCTC 

SRA  3  gggacttgccgcagggtccccacagcgcctcccgacggttgtgaaccagcatccatcctc 

180 

SRA  /  cacggattccggcaacccccctggccctggacgtgtctcaactggcccgcgtcaggggcc 
SRA  2  cacggattccggcaacccgcctggccctggacgtgtctcaactgccccgcgtgaggggcc 
SRA  3  cacggattccggcaacccgcctggccctggacgtgtctcaactggcccgcgtgaggggcc 


SRA 
SRA  I 
SRA  2 
SRA  3 


SRA 
SRA  I 
SRA  2 
SRA  3 


GCCCCGGAA 

A  T  G 

ACGCGCTGCCCCGCTGGCCAAGCGGAAGTGGAG 

A  T  G 

GCCCCGGAA 

A  T  G 

ACGCGCTGCCCCGCTGGCCAAGCGGAAGTGGAG 

A  T  G 

GCCCCGGAA 

A  T  G 

ACGCGCTGCCCCGCTGGCCAAGCGGAAGTGGAG 

A  T  G 

Metl 


Met  2 


240 

gcggagctgtac 

gcggagctgtac 

gcggagctgtac 

Igcggagctgtac 

300 


gtgaagccgggcaacaaggaacgcggctggaacgacccgccgcagttctcatacgggctg 

gtgaagccgggcaacaaggaacgcggctggaacgacccgccgcagttctcatacgggctg 

gtgaagccgggcaacaaggaacgcggctggaacgacccgccgcagttctcatacgggctg 

gtgaagccgggcaacaaggaacgcggctggaacgacccgccgcagttctcatacgggctg 


SRA 
SRA  / 
SRA  2 
SRA  3 


338  348  360 

cagacccaggccggcggacccaggcgctcgctgcttaccaagagcgtagccgcaccccag 
cagacccaggccggcggacccaggcgctcgctgcttaccaagacggtagccgcaccccag 
cagacccaggccggcggacccaggcgctcgctgctt  a(T)c  a  a  g  a  g  g  g  t©g  ccgcaccccag 
cagacccagccccgcggacccaggcgctcgctgcttaccaagagcgtagccgcaccccag 


SRA 
SRA  J 
SRA  2 
SRA  3 


SRA 
SRA  / 
SRA  2 
SRA  3 


SRA 
SRA/ 
SRA  2 
SRA  3 


SRA 
SRA  1 
SRA  2 
SRA  3 

SRA 
SRA  / 
SRA  2 
SRA  3 


gatggatcccccagagtccccgcatcagagacttctcctgggcctccccca 

gatggatcccccagagtccccgcatcagagacttctcctgggcctccccca 

gatggatcccccagagtccccgcatcagagacttctcctggccctccccca 

gatggatcccccagagtccccgcatcagagacttctcctgggcctccccca 


A  T  G 

G  G  G  C  C  T 

A  T  G 

G  G  G  C  C  T 

A  T  G 

G  G  G  C  C  T 

A  T  G 

G  G  G  C  C  T 

ccacctccttcaagta 

CCACCTCCTTCAAGTA 

ccacctccttcaagta 

CCACCTCCTTCAAGTA 


A  G  G  C  T 
A  G  G  C  T 
A  G  G  C  T 
A  G  G  C  T 


SRAU2 


Meta 

CCCAGGTCCCCACCTGTGGGGAGTGGTCCTGCCT 

CCCAGGTCCCCACCTGTGGGGAGTGGTCCTGCCT 

CCCAGGTCCCCACCTGTGGGGAGTGGTCCTGCCT 

CCCAGGTCCCCACCTGTGGGGAGTGGTCCTGCCT 


GTGGAGCCCACAAGTTTCCCA 

GTGGAGCCCACAAGTTTCCCA 

GTGGAGCCCACAAGTTTCCCA 

GTGGAGCCCACAAGTTTCCCA 


520 

GTCGAGTCTGAGGCTG- 
GTCGAGTCTGAGGCTG- 
CTCGAGTCTGAGGC  T  G  ♦ 
GTCGAGTCTGAGGC  TIS©C 


TGATGGAGGATGTGC 
TGATCGAGGATGTGC 
TGATCGAGCATGTG C 


CCTTTGGAA 

CCTTTGGAA 

CCTTTGGAA 

CCTTTGGAA 


CAGGCATTGGAAGACTGC 

CAGGCATTGGAAGACTGC 

CAGGCATTGGAAGACTGC 

CAGGCATTGGAAGACTGC 


SRAL2 


)tgatgcaggatgtg)c 

SRAL-GTCG 

cgtggccacacaaggaaccaggtatgtg 

CGTCGCCACACAAGGAACCAGGTATGTG 

cgtggccacacaaggaaccaggtatgtg 

CGTGGCCACACAAGGAACCAGGTATGTG 


atcagccgacgcctggcactgctgcaggaacagtgggctggaggaaacttgtcaa 

ATCAGCCGACGCCTGCCACTGCTGCAGGAACAGTGGGCTGGAGGAAAGTTGTCAA 

atcagccgacgcctggcactgctgcaggaacagtgggctggaggaaagttgtcaa 

ATCAGCCGACGCCTGGCACTGCTGCACGAACAGTGGGCTGGAGGAAAGTTGTCAA 


480 

C  T  G  G  C 
C  T  G  G  C 
C  T  G  G  C 
C  T  G  G  C 

540 

T  G  A  G  A 
T  G  A  G  A 
T  G  A  G  A 
T  G  A  G  A 

600 

A  T  G  A  C 
A  T  G  A  C 
A  T  G  A  C 
A  T  G  A  C 

660 

T  A  C  C  T 
T  A  C  C  T 
T  A  C  C  T 
T  A  C  C  T 


720 

SRA  GTAAAGAAGAGAATGGCTCTACTGGTGCAAGAGCTTTCAAGCCACCGGTGGGACGCAGCA 
SRA  /  GTAAAGAAGAGAATGGCTCTACTGGTGCAAGAGCTTTCAAGCCACCGGTGGGACGCAGCA 
S/iA  2  GTAAAGAAGAGAATGGCTCTACTGGTGCAAGAGCTTTCAACCCACCGGTGGGACGCAGCA 
SRA  3  GTAAAGAAGAGAATGGCTCTACTCGTGCAAGAGCTTTCAAGCCACCCGTGGGACGCAGCA 

780 

SRA  GATGACATCCACCGCTCCCTCATGCTTGACCATGTGACTGAGGTCAGTCACTGGATGGTA 
SRA  I  GATGACATCCACCGCTCCCTCATGGTTGACCATGTGACTGAGCTCAGTCAGTGGATGGTA 
SRA  2  GATGACATCCACCGCTCCCTCATGGTTCACCATGTCACTGACGTCAGTCAGTGGATGGTA 
SRA  3  GATGACATCCACCGCTCCCTCATGGTTGACCATGTGACTGAGGTCAGTCAGTGGATGCTA 

K40 

SRA  CGAGTTAAAAGATTAATTGCAGAAAAGAGGAGTCTGTTTTCACACCAGGCAGCCAATGAA 
SRA  I  GGAGTTAAAAGATTAATTGCAGAAAAGAGGAGTCTGTTTTCAGACGAGGCAGCCAATGAA 
SRA  2  GGAGTTAAAAGATTAATTGCAGAAAAGAGGAGTCTGTTTTCACAGCAGGCAGCCAATGAA 
SRA  3  GGAGTTAAAACATTAATTGCAGAAAAGAGGAGTCTGTTTTCAGAGCAGGCAGCCAATGAA 
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GAGAAATCTGCAGCCACAGCTGAGAAGAACCATACCATACCACGCTTCCAGCAGGCTTCA 
SRA  /  CAGAAATCTGCAGCCACAGCTGAGAAGAACCATACCATACCAGGCTTCCAGCAGGCTTCA 
SRA  2  GAGAAATCTGCAGCCACAGCTCAGAAGAACCATACCATACCACGCTTCCAGCAGGCTTCA 
SRA  3  GAGAAAT  CT  GCAGCC  ACACCT  GAGAAGAACCATACCATACCAGGCTTCCAGCAGGCTTCA 
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